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ON THE INTERPRETATION OF SPIRAL STRUCTURE 
IN THE NEBULAE* 


BERTIL LINDBLAD 


ABSTRACT 


The development of spiral structure in the nebulae is discussed in adherence to 
principles expressed in earlier papers but with a strenghthening of some of the argu- 
ments and renewed comparison with observed conditions. 

The series of types of increasing flattening among the extragalactic nebulae is inter- 
preted as a series of increasing angular momentum when this quantity is considered in 
relation to the mass and mean density of the system. An internal condensation leading 
to an increase of flattening may be the formation of small solid particles in gaseous 
material and the formation of the stars themselves in an original extended medium. 

The instability of matter at the edge of a rotating stellar system of great flattening 

is discussed, and the mechanism by which a spiral arm may be formed in a steady 
process of disturbed motion is described, with reference to the properties of the char- 
acteristic diagram of the system. The reasons for the natural selection of a type of dis- 
turbance giving two spiral arms, extending from two opposite points which are fixed in 
space, are explained, and the possible analogy with the sectorial harmonic wave which 
_ makes the Maclaurin ellipsoid unstable at a certain flattening is mentioned. 
Previous work had indicated that the full development of spiral structure is likely 
to give an approximately logarithmic shape of the arms, in which the particles of the arm 
describe roughly circular orbits of uniform angular velocity around the center. It is conclud- 
ed here that such a formation is possible. Still, the process of segregation leading to such 
a formation is perhaps of a more complicated nature than indicated by the schematic 
theoretical picture. 

The characteristics of the ‘‘amorphous”’ central system of a nebula, its successive 
formation, the physical character of the individual stars—which appear to be dwarf 
stars of late type, in contrast to the more massive white stars in the spiral arms— 
and the angular velocity of rotation are discussed. The normal case is probably the 
one where the nuclear regions rotate with a higher angular speed than the spiral struc- 
ture. An exception is perhaps found in the barred spirals, where the ‘‘bar’’ is possibly to 
be considered as a remnant of a kind of interior spiral structure and has, presumably, 
a low angular speed of rotation. 


* This paper was presented at the symposium on galactic and extragalactic struc- 
ture, held in connection with the dedication of the McDonald Observatory on May 5-8, 


1939. 
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The verification of the theoretical results by observational data is discussed in some 
detail. H. W. Babcock’s results concerning the rotation of the Andromeda nebula are 
in good agreement with the main theoretical result concerning the possibility of a 
uniform rotation of the spiral structure. 

The question of orientation of the nebular plane is discussed for the Andromeda 
nebula, for M 81, and for twelve typical nebulae which show different degrees of de- 
velopment of the central mass. The formation of characteristic lanes on the near and 
the farther sides of the nebular disk is discussed. The relative arrangement of bright 
and dark matter in nebulae of small central mass, the polarization of light in the dark 
clouds as observed by Ohman, and the different illumination of dark clouds on the near 
and farther side, respectively, are factors of great importance in this connection. The 
evidence seems to be in favor of the theoretical result that the rotation is in the direc- 
tion in which the spiral arms wind when followed outward from the center. In the case 
of “‘spindle’’ nebulae, where the apparent direction in which the spiral arms wind can- 
not be determined directly, the intensity distribution along the major axis may indicate 
the side of concave and convex arms, respectively, as seen by the observer. 

A method of deriving the direction of rotation independently of the dark lanes in 
a tilted nebula is afforded by the phenomenon of a spiral arm breaking up into a fan- 
shaped structure of secondary branches. In the absence of detailed observations of the 
velocities of rotation for the individual branches, the assumption of a uniform value of 
the area constant has been made. The arms of the nebula M ror then give a direction 
of rotation in agreement with the theory. 


The main principles of the theory of the development of spiral 
structure in the nebulae which will be discussed below have been 
presented in previous papers.’ A review of the general line of thought 
of the theory will be given here, with references to previous papers 
for some of the details, but in a development which to a considerable 
extent brings new points of view into the discussion. The important 
question of the observational test of some of the fundamental con- 
clusions will be approached as far as appears to be possible at pres- 
ent. 

A fundamental fact which is of the greatest importance for an in- 
terpretation of the phenomenon of spiral structure is the appearance 
of the spiral type of nebula as a natural continuation of a series of 
elliptical types of increasing degree of flattening. This feature has 
been recognized very fully by Jeans? in his theory of the development 
of spiral structure. It lies very near at hand, from a general theo- 
retical point of view, to assume that this series of types corresponds 

* “Cosmogonic Consequences of a Theory of the Stellar System,” Ark. mat., astr. och 
fys., 19A, No. 35, Upsala Medd., No. 13, 1926; ‘“‘On the Spiral Orbits in the Equatorial 
Plane of a Spheroidal Disk,”’ K. Sv. Vetenskapsakademiens Handlingar, 4, No. 7, U psala 
Medd., No. 31, 1927; ‘‘Contributions to the Theory of Stellar Systems,” Stockholm Ann., 
12, No. 4, 1936; ‘“‘On the Theory of Spiral Structure in the Nebulae,” Zs. f. Ap., 15, 124, 
Stockholm Medd., No. 40, 1938. 

2 Problems of Cosmogony and Stellar Dynamics, Cambridge University Press, 1917; 
Astronomy and Cosmogony, Cambridge University Press, 1928. 
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to an increase of the angular momentum, considered in a certain 
relation to the mass and the mean density of the system. In the case 
of the series of homogeneous Maclaurin ellipsoids, if H is the angular 
momentum, M the mass, and p the density, the quantity 


(1) 


is a function only of the flattening of the system. We may assume 
that the increase of flattening in the case of stellar systems corre- 
sponds to an increase of an analogous quantity, where p is a repre- 
sentative mean density. It is possible that one and the same mass, 
with a given value of the angular momentum H, may go through 
successive forms of equilibrium on account of an internal process of 
condensation. We may have in mind here a condensation of inter- 
stellar gas into small, solid particles and even the formation of the 
stars themselves in an original, extended medium. We must con- 
sider in the first place, however, that in all probability the systems 
were ‘“‘born’’ with a wide range in angular momentum and that the 
variation of type is thus to a great extent dependent on original con- 
ditions. 

When treating the dynamics of the nebulae, we may consider the 
mass motions of matter in a way which is in several respects analo- 
gous to the developments of hydrodynamics and the dynamics of 
gases. A fundamental difference, however, lies in the circumstance 
that in stellar systems the motions of the individual particles must 
be considered as free motions under the gravitational forces of the 
system. The assumption made in hydrodynamics that a small 
closed surface, supposed to move with the fluid, will always inclose 
the same matter cannot be true here, even approximately, for any 
great length of time. 

The special clue to an understanding of the development of spiral 
structure which forms the foundation of the present theory starts 
from the properties of individual motions close to the edge of a 
flattened system. Moreover, this theory has immediate connection 
with the theory of galactic rotation. We know that in our own galac- 
tic system there is, according to this theory, a concentration of mat- 
ter toward the circular motions. When examining the circular mo- 
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tion at the edge of a system of great flattening toward the plane of 
rotation, which we assume to be an equatorial plane of symmetry, 
we find that at a certain degree of flattening the circular motion at 
the edge becomes unstable, and by an infinitely slight disturbance 
the particle will depart in an asymptotic orbit which at first winds in 
a close spiral around the system but at greater distances departs 
rapidly from the system. After attaining a certain maximum dis- 
tance, the orbit returns symmetrically toward the original circular 
motion. In the case of a spheroidal, homogeneous system the insta- 
bility begins at the ratio 0.55 between the polar and equatorial axes 
of the spheroid, corresponding to an eccentricity of the meridian sec- 
tion amounting to 0.834. With increasing flattening, the maximum 
distance of the orbit increases; for an infinite degree of flattening, the 
maximum distance extends to infinity. If the system has symmetry 
of rotation, as well as an equatorial plane of symmetry, the compo- 
nent of force at right angles to the equatorial plane close to this plane 
may be written 


OV 


where G is the constant of gravitation, and z is the elevation above 
the equatorial plane. The condition for the appearance of the insta- 
bility then is 


<p-—p, (3) 


where w, is the angular velocity of the circular orbit and p is the 
density of matter at the point considered. This relation may be 
written 


we 


(4) 


Dild 


The first term may be assumed to decrease continuously with in- 
creasing degree of flattening of the system, in analogy with a corre- 
sponding quantity in the case of the equilibrium figures of liquid 
spheroids. The second term depends more on local conditions. A 
low value of this term is likely to occur if there is a rapid increase of 
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density with distance in the equatorial plane in the neighborhood of 
the point considered. This may be fulfilled in particular if there is a 
strong radial density gradient in the equatorial plane close to the 
point in question. 

The instability may be fully described by the so-called ‘‘character- 
istic diagram,” which has the energy /, and the area velocity J, of a 
particle as axes. If II, 0, and Z, are the velocity components in the 
radial, tangential, and vertical directions, respectively, and if, fur- 
ther, V is the potential and R the distance from the axis of rotation, 
we have 


,L=P+@+ 27 — 2V 

The region in the diagram which represents possible combinations 
of J, and J, at a given point in the system is bounded by a certain 
parabola, and the region of possible combinations of J, and J, for the 
system as a whole is defined by a certain curve, the ‘“‘characteristic 
envelope-curve,”’ which ordinarily represents the envelope of the 
characteristic parabolae for successive values of R in the equatorial 
plane. The point of contact between a characteristic parabola and 
the envelope represents the circular motion in this plane for the 
value of R in question. For homogeneous spheroidal systems the 
envelope-curves for the interior part of the system are straight lines, 
as shown in Figure 1. In a more general case, which is illustrated in 
Figure 2, the angle ¢ between the tangent of the envelope-curve and 
the /, axis obeys the relation 


tan g = 24, (6) 


where w, is the angular velocity of the circular orbit for the value of 
7, in question. In the case of an instability at the edge, the envelope- 
curve of the outer regions recedes “backward” from a turning-point 
P, to another turning-point P,, so that in this interval J, and J, de- 
crease with increasing Rk, whereas for the regions of stable circular 
orbits, 7, and /, increase with R. In this case the /,, 7, diagram may 
be divided into three parts,’ two situated below the singular part 


3 Cf. Bergstrand festskrift, p. 23, 1938. 
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P,P, of the envelope-curve and one above this part of the curve. 
Points in the last-mentioned region represent particles which under 
some circumstances may depart from the central system in wide 
orbits. 

We must now consider that the real significance of the develop- 
ment of the singular course of the envelope-curve will be that a great 
part of the matter situated close to the edge of the system moves in a 
state close to instability, though still in stable orbits. It should be 
expressly remarked that a point in the diagram above P,P, does 


I, 


Fic. 1 Fic. 2 
Fic. 1.—The characteristic diagram for homogeneous spheroids 


Fic. 2.—The characteristic diagram for a general system of rotational symmetry and 
great flattening. 


not necessarily, at least for a considerable time, correspond to an 
orbit of ejection. If the particle has a finite motion in s, it may still 
remain entirely within a region of stable orbits, passing through the 
equatorial plane with a finite value of the velocity component Z. 
The successive development of an instability of the kind in question, 
however, would lead to an ejection of a certain amount of matter 
across the circle where the instability begins. In the first instance 
this would tend to smooth the radial density gradient at this circle 
and thereby to counteract the development of the instability by 
making p/p increase in the region of instability. A further increase 
of the flattening, however, leading to a further increase of p and 
thereby to a decrease of the first term in equation (4) and a restora- 
tion of the instability, would not mean an immediate ejection of any 
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| 
| i le 

I, fe 
E, / i 

E, 


SPIRAL STRUCTURE IN NEBULAE 7 


considerable amount of new material, as the matter likely to be 
ejected at the circle of instability has already been supposed to have 
crossed this circle. A further development of the instability by fur- 
ther increase of the flattening of the system would make all the eject- 
ed matter recede to greater distances. For this reason and by the 
increased density of the central system due to the increase of flatten- 
ing and the general contraction, we may expect in the first instance 
the formation of a central lens with a fairly sharp edge and at the 
same time around the central lens the appearance of a faint ring- 
shaped structure formed by the ejected matter. According to ob- 
servations by Hubble,‘ this seems to agree well with what is actually 
observed in those nebulae which represent the transition from the 
elliptical to the spiral type. , 

The ejection of maiter will not, in general, be perfectly evenly dis- 
tributed around the edge, and on account of this fact the assembled 
mass of the ejected matter will exert a complicated disturbing action 
on the central system. By the law of areas in the motion of the 
ejected matter this will have a general tendency to recede along the 
edge relative to the matter following the circular motion here, in a 
direction opposite to the direction of the circular motion. According 
to Hubble, it seems as if the dark matter in the outer parts of the 
central lens would be most sensitive to the disturbing forces, which 
might possibly be attributed to an increased radiation pressure on 
such matter helping to upset the balance by reducing the net central 
force. In an advancing disintegration we are most interested in a 
type of ejection which is of a steady character, in which the ejection 
continues steadily by the tidal action of the previously ejected mat- 
ter. A possible disturbance of this type is shown in Figure 3. In this 
case the point of ejection is supposed to be fixed in space. The rela- 
tive displacement of the ejected matter from its original position at 
the edge to the positions in the “‘arm”’ causes a tangential disturbing 
force at 7. We may now construct the characteristic envelope-curve 
along the path OTS and compare it with the envelope-curve along 
the radius under undisturbed conditions (£, in Fig. 3). The curve 
Eors will be depressed relative to E, in about the way indicated. In 


4Communicated by Dr. Hubble in a lecture at the dedication of the McDonald 
Observatory in May, 1939. 
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particular, the increased depression of the “‘receding”’ part of the en- 
velope-curve will cause a considerable amount of matter arriving 
at T with representative points above this part of the curve to follow 
the arm TS in orbits of ejection. The diagram explains the way in 


Fic. 3.—The formation of a spiral arm 


which the edge of the central system is successively ‘“‘cut off” to form 
the arm. 

We may ask if other assumptions regarding the speed of the point 
of ejection T are possible. If we introduce a rotating co-ordinate 
system moving with the angular speed w, in order to follow the mat- 
ter at the edge, we have assumed above that the point of ejection 
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moves in the rotating co-ordinate system with the relative angular 
speed w in the direction against the rotation. It is clear that the 
speed actually selected will depend not only on the kind of tidal 
force which is developed as a consequence of this speed but also on 
the resonance of the disturbing action with the free oscillation of the 
particles. If the system has approximately uniform angular velocity 
of rotation w, the speed of internal density variations in a co-ordinate 
system following the rotation will be approximately 2.’ If the dis- 
turbance is in resonance with this oscillation, its intensity at the 
edge, which we may measure by the intensity of the tangential force 
in the direction of rotation, may be developed in a series 


(As cos [2wt + + B, sin [2wt + n6]) . (7) 


n 


For disturbances proceeding against the velocity of rotation, m is 
positive. The value 2 = 1 gives a disturbance with only one maxi- 
mum along the circumference, proceeding with the speed 2w against 
the velocity of rotation. It does not seem that an arm produced by 
such a disturbance would, in general, give a sufficient disturbing 
force.° The case n = 2, however, corresponds to a disturbance with 
two opposite maxima, proceeding with the relative speed w against 
the velocity of rotation, i.e., fixed in space. The contraction of the 
matter in the arm due to the law of areas in the motion from the 
fixed point of ejection will cause a considerable disturbance at this 
point, as explained above in reference to Figure 3. For higher values 
of the disturbance becomes more complicated, as the number of 
maxima around the circumference increases. It seems, therefore, 
probable that the tidal force will be most effective in the case n = 2, 
in which case we should expect the appearance of two diametrically 
opposite spiral arms, in agreement with observation. 

If an ejection process of the kind just described is present at a 
certain moment, it may either decrease or increase with time. In 
the former case the system is stable for a disturbance of this type; in 
the latter, unstable. It is evident that the instability will be favored 


5 Cf. Zs. f. Ap., 15, 132, 1938. 
6 An exception may be a nebula with one main arm like NGC 4725 on PI. II. 
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by a great flattening of the system, which gives asymptotic orbits of 
wide extension. This may be seen also by reference to the sectorial 
harmonic waves in the case of a rotating homogeneous ellipsoid. The 
proceeding disturbance n = 2 in our case has a certain analogy to 
the two-dimensional, sectorial harmonic wave of two opposite maxi- 
ma, which makes the Maclaurin ellipsoid ‘ordinarily’ unstable at 
the eccentricity e = 0.9529 of the meridian section.’ The spiral 
arms may perhaps be considered as the advancing front of such a 
wave of steadily increasing amplitude, proceeding with higher angu- 
lar speed than the wave as a whole, which in the case of the liquid has 
the speed w/2 against the direction of rotation. We may then expect 
the spiral phenomenon to take place at a flattening corresponding to 
about the eccentricity e = 0.95 of the meridian section. This corre- 
sponds to the ratio of the axes in this section equal to 0.3 and seems 
to agree well with the limiting flattening of the elliptical nebulae 
where these merge into the spiral type. 

In the case of nebulae of fairly thin arms, like M 81, we should ex- 
pect the arms to follow the asymptotic orbits from the circular mo- 
tion at the edge. It is actually found that the arms of M 81 corre- 
spond very well to the asymptotic orbits for an eccentricity e = 
0.965 of the meridian. When the disintegration of the central sys- 
tem becomes more complete, however, the spiral arm no longer rep- 
resents the orbit of an individual particle. On account of the succes- 
sive decrease of the central system, the particles must instead be ex- 
pected to describe approximately circular orbits entirely outside of 
the remaining nucleus. Under these circumstances I have shown‘ 
that a roughly logarithmic form of the arms may be expected and 
that the speed of rotation of such a spiral structure may be approxi- 
mately uniform, if the original system was one of uniform angular 
speed. 

In the asymptotic orbit extending from the edge of a flattened 
system, the ratio of the centrifugal force at the point of maximum 
extension to the centrifugal force at the edge of the central system 
will be f, if f is the ratio of the radii vectors for the two positions. 
The ratio of the central forces will only slightly exceed the value 


7 Bryan, Phil. Trans. R. Soc. London, 180, 187, 1889. 
8 Zs. f. Ap., 15, 133, 1938. 
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f-3, which makes the orbit at the maximum distance adhere rather 
closely to a circular orbit. The slight excess of the ratio over f~3 
makes the orbit return toward the center. In the case of a develop- 
ment of heavy spiral structure, however, the contraction of the re- 
maining central system and the decrease of its mass due to the de- 
parture of matter in the spiral arms will mean a decrease of the cen- 
tral force at the maximum distance of the orbit, and the ejected mat- 
ter will, therefore, no longer return to the central system. The bal- 
ance of centrifugal and central 
force in approximately circular 
motion may thus be restored in 
the spiral structure. 

The logarithmic shape of the 
arms, 


r= Ce, (8) 


follows if we assume the dis- 
tances from the center in the 
spiral structure to be propor- 
tional to the distances from Since 

the center in the original sys- 

tem and the radial thickness of the arm to be proportional to 
these distances, taking into account also the contraction of the 
arm in its own direction at larger distances from the center, which 
follows by the law of areas in the motion. If the original system has 
roughly uniform angular velocity of rotation, this will evidently be 
the case also for the spiral structure under these assumptions. That 
the central force in the spiral structure will be roughly proportional 
to the distance from the center may be seen in the following way. In 
the schematic picture of the two spiral arms in Figure 4, we consider 
the radial component of force at two consecutive intersections P,, 
and P,,, of an arm with the vertical line through the center. We 
denote the force at the lower point with F,,, at the upper point with 
F,,:. We further denote with F;,, the force at the upper point when 
we prolong the spiral arms as shown by the dotted curves. Assuming 
the radial thickness of the arm to be proportional to the distance from 
the center and the thickness at right angles to the plane of the figure 
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to be small and one and the same throughout the structure, we 
should have 
ntt = F,, (9) 


because the difference between the two cases is one of scale only. 
The mass of similarly situated parts in the two cases increases as the 
square of the scale factor and thus balances everywhere the effect 
of the scale factor on the inverse square of the distances. In order 
to obtain F,,,,, we must subtract from F;,,, the gravitational force 
due to the dotted parts of the arms. We propose to calculate first 
this force on a point P situated on the vertical line at the distance a 
from the center. We take the positive x» axis upward in the figure 
and count the angle @ from this direction. We denote the distance 
from the center for a point on the spiral arm with 7, and the distance 
from the point P with p. The mass dm of a part of the spiral arm 
corresponding to an interval d#@, when we assume the radial thick- 
ness of the arm to be proportional to r, may be expressed 


dm = wurde, (10) 


where yp is a certain constant. If the spiral structure expands, pu will 
keep its value unchanged on account of the constancy of r’dé@ in the 
motion. The force G, in the x direction exerted on P by the upper 
dotted part of the spiral arm will be 


2 
G, -{f w(x — a)d0, (11) 


ps 


where 


p = (r? + a? — 2ar cos r= x= 6. 
We may develop in powers of a/r and assume a/r < 0.4, so that 
a\? a 
(2) —1. 
r r 


In this case we know that the rest term in the development of p™ 
will be of the order of the first term omitted. In the present case the 
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terms of even order in a have no influence on the final result, and 
the third-order terms may be shown to have only small influence. 
We are, therefore, entitled to consider, as an approximation, a de- 
velopment keeping only the terms of the first order, and we get 


a 3a 
G, (cos ce + cos 0) dé. 
In the third term, which contains cos? 6 as a factor, we may with very 
good approximation take together numerically equal @ of opposite 
sign and for the factor e~® take the geometrical mean, which is equal 
to unity. Then by performing the integration, we obtain 


a 
G,; = 2p + (e~e(*/2) ete(n/2) + 


The force G, due to the lower dotted part of the spiral arm is ob- 
tained by changing the sign of G,, after changing the sign of a. Add- 
ing together, we get 


2ua 


[Bem — (ete(/2) — = ka. (12) 


G,+ G, = 


Thus, in the first approximation the combined gravitational force of 
the dotted parts of the spiral arms in Figure 4 is directed outward 
along the radius and is proportional to the distance from the center. 
We have, then, in particular for the points on the spiral arms con- 
sidered above 


att = Fat, + 
or, by equation (9), 
Pais — Fy = — bays = — — (13) 
where the constants k and x obey the relation 


k 
(14) 
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We may now repeat the process for two similarly situated points 
which we denote by P,,_, and P,, and obtain, in an analogous way, 


Assuming the two spiral arms to wind indefinitely inward to an in- 
finitely small distance from the center, we may put 


| 
| 


— Fa = — Tn) 
Thanks to the small but finite thickness at right angles to the spiral 
plane, we may assume F,,_; to vanish with r,_; for very large values 
of 7. Equations (15) are still valid without perceptible changes of x. 
Adding together equations (15), we get as the general expression for 
the radial force 


F=—kxr (16) 


for points along the arms. That the result will hold approximately 
also for other sections through the center than the one just consid- 
ered can easily be verified. 

We may now calculate, on the same assumptions, the change of 
the radial force F, on one and the same particle of an arm by an ex- 
pansion of the spiral structure in the ratio f, during a motion obeying 
the law of areas. The coefficient c will increase in the ratio f? to c.f’, 
say, if c, is an initial value. The coefficient C, when 6 is reckoned as 
before, will obey the relation 


The coefficient & will change as follows: 
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The force F, = — xr, when r increases in the ratio f, is finally found 
to be 
2 


— (e%f2(4/2) — (2/2) —I 


(19) 


For observed spiral arms, which according to Danver correspond in 
the mean nearly to cr = 1, the rate of change of the central force 
according to formula (19) lies between f~? and f~3, but for slightly 
larger values of c the rate of change becomes much accelerated and 
more rapid than f-$. According to equation (19), F, may even pass 
zero and become positive for very great values of ¢,f?. Though this 
cannot be true for the outer regions in higher approximations, it 
seems that we may draw the conclusion that the spiral structure 
would become unstable for very large values of ¢.f?. 

Though we must keep in mind that the numerical results will of 
course be subject to modifications if the schematic assumptions re- 
garding the mass distribution are changed, the considerations above 
show the harmony between, on the one hand, a logarithmic shape of 
the arms and a uniform speed of rotation deduced from the process of 
disintegration and, on the other hand, the gravitational forces to be 
expected along the arms in a fully developed spiral structure. The 
assumption of uniform thickness of the arms in the direction at right 
angles to the spiral plane may be justified by the circumstance that 
the mass detached to form the spiral arms probably constitutes a 
higher percentage of the entire mass in the outer regions of the sys- 
tem than in the inner regions. On the other hand, for this reason we 
must take into account a remaining central mass not included in the 
spiral structure. In the case of very well-developed spiral structure, 
however, we assume this remaining central mass to be small compared 
with the mass of the arms, which seems to be a reasonable assump- 
tion from the actual appearance of such nebulae. 

The adjustment assumed above between the centrifugal force and 
the central force at the maximum distance in the asymptotic orbit 
is, of course, a somewhat schematic assumption. An adjustment to 
balance in circular motion may mean considerable shifts in the radius 
vector. Furthermore, it is not to be excluded that the final motion at 
various points may differ perceptibly from the circular motions and 
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may include pulsating motions in the radial direction. The value of 
the coefficient c in equation (8), when determined from observed 
spirals, agrees, however, with the assumption of an adjustment into 
circular orbits at distances from the center which, at least to the 
order of magnitude, agree with the maximum distances in asymp- 
totic orbits. We have, approximately, 


c= ab’, (20) 


where a is the ratio between the radial thickness of the arm and the 
original distance from the center and b is the ratio between the dis- 
tance from the center in the spiral structure and the original dis- 
tance from the center. If vis the angle between the radius vector and 
the tangent of the arm, we have 


(21) 


According to Danver,’ v is close to 74°, which gives c close to 3. 

We may also look upon this problem from the point of view of the 
virial theorem. If a system with constant energy settles over from 
one equilibrium form to another, the potential energy must be the 
same in the two cases. The increase of potential energy due to the 
expansion in radius vector will be partly compensated by the de- 
crease of potential energy due to the contraction of the arms in their 
own direction, which is a consequence of the law of areas, but to some 
extent the compensation will also be due to a contraction of the re- 
maining central mass. Such a contraction must follow when a con- 
siderable part of the outer regions departs to greater distances, which 
must give an increase of the central force on the remaining matter. 
The angular speed of the central mass will increase by the contrac- 
tion, which will, in a way, favor the continued development of spiral 
structure in this mass, but the angular momentum of the central 
mass will, nevertheless, decrease successively, and at last the cen- 
tral system remains only as an “amorphous” elliptical mass. This 
mass may have a high angular speed of rotation, but in spite of this, 
by a high degree of concentration toward the nucleus it will possess 


9 Vierteljahrschrift, 72, 340, 1937. 
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too low an angular momentum to take on a very much flattened 
figure of equilibrium. We thus get the typical picture of the nebula 
as a spiral structure surrounding an apparently “amorphous” cen- 
tral system of low degree of flattening and of strong concentration 
toward a dense nucleus. The characteristic difference in color be- 
tween the spiral arms and the nuclear regions can then be explained. 
The nuclear regions contain stars which have formed outside of the 
equatorial plane, often with smaller original velocities of rotation 
than the stars which follow circular orbits in this plane at the same 
distance from the center. The former stars, which have a tendency 
to condense toward the center, are mostly dwarf stars of advanced 
color. The abundance of matter with low internal velocity in the 
neighborhood of the equatorial plane, on the other hand, favors the 
development of massive stars of high surface brightness in these re- 
gions. These heavy stars thus belong to the rotating “‘subsystem” 
which contributes most of the material to the spiral structure. 

We have considered heretofore a uniform angular motion of a 
fully developed spiral structure as a quality which may follow from 
the geometrical and dynamical properties of a process of dissolution 
of an original closed central system, when the process takes place 
from the edge inward. We must consider, however, that the picture 
thus given of the dissolution into spiral structure is possibly a con- 
siderable simplification of the real process as it actually occurs in 
nature. It is probable that in many cases the dissolution begins at 
various points in the original rotating mass, at different distances 
from the center where instabilities of motion originate. The spiral 
arms may develop to more and more distinct features and catch in 
their whorls the matter of greatest angular momentum and of low 
internal velocity dispersion, by a process similar only in its broad 
lines to the one which we have described. It is possible that in this 
way an exchange of momentum over different intervals of distance 
from the center may take place along the spiral arms in such a way 
that in a regular, fully developed spiral structure the uniformity of 
angular speed becomes well realized. 

This does not mean, however, that one and the same angular 
speed prevails throughout the entire system. We have seen that the 
central regions where the spiral structure is not fully developed are 
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likely to possess a higher angular speed than the spiral structure. 
This should apply especially to the elliptical mass immediately sur- 
rounding the nucleus, which in several nebulae is known to have a 
high and uniform angular speed. The region of transition between 
the spiral structure and the central mass may be expected to show 
complicated conditions. In the Andromeda nebula this region shows 
traces of spiral structure, especially in the dark markings, which is 
probably due to direct tidal effects of the surrounding main spiral 
arms. A complicated variation of the observed angular speed with 
distance from the center may be expected here, as has actually been 
observed at the Lick Observatory by H. W. Babcock.*’ 

In the case of the “barred” spirals the conditions in the nuclear 
regions may deviate still more strongly from normal conditions. 
With “normal conditions” we mean here the case when the nuclear 
regions show an increased angular speed compared with the spiral 
structure. It is possible that the “bar” is a remnant of a kind of 
inner spiral structure, more developed than in the case of the An- 
dromeda nebula, and perhaps caused by, but largely independent of, 
the outer spiral structure. The degeneration of this inner spiral 
structure into a straight bar, presumably by a process of extension, 
would suggest that the angular speed of rotation of the bar structure 
is low, perhaps equal to or even smaller than the rotation of the outer 
structure. 

We shall at last proceed to consider the verification by observa- 
tion of some of the vital points in the theory developed above. The 
conclusion that a uniform angular speed of rotation may prevail 
in a fully developed spiral structure agrees well with Babcock’s ob- 
servations of the variation of the radial velocity along the large axis 
of the Andromeda nebula. In particular, the very high velocities ob- 
served by him for a couple of gaseous nebulae situated outside of the 
main spiral structure are of extremely great interest. The fact that 
these nebulae seem to belong to secondary branches of the main 
spiral arms is important in this connection. If they have the same 
area constant in their motion around the center as the neighboring 
parts of the main arm, we should expect the linear velocities of the 
extreme parts of the main arm to be still higher, which would actual- 


10 Pub. A.S.P., 50, 174, 1938. 
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ly improve the linearity of the variation of the velocity with the dis- 
tance from the center in Babcock’s diagram. The variation in the 
region of transition between the nuclear regions and the main spiral 
structure has been found" to stand in intimate connection with the 
distribution of light in the region in question. If we multiply the in- 
tensity J with the distance R from the center, this product has a 
minimum where Babcock finds the minimum velocity. It is possible 
that this fact may have something to do with the diffuse spiral 
structure in this region and that the matter left between the nuclear 
region and the diffuse spiral arm has a low velocity of rotation. It is 
also possible that we are concerned here with a higher level in the 
nebula. The minimum of /R may be due partly to interior absorp- 
tion in the nebula, which prevents us from penetrating here to the 
more rapidly rotating equatorial layers of matter. In this connection 
it may be added that the importance of the absorption of light in the 
interior of stellar systems for the interpretation of spectrographic 
observations of internal motions has recently been emphasized by 
I. Holmberg.’? He shows how a strong internal absorption in a 
nebula seen edgewise may cause an apparently uniform rotation. In 
the case of the nuclear regions of the Andromeda nebula it does not 
seem probable, however, that the linear variation of velocity may be 
caused to any great extent in this way. The decreasing flattening of 
the isophotes surrounding the strong nucleus (Fig. 5) suggests that 
the nucleus is fairly transparent for the angle of incidence under 
which we see the nebula. On the whole, the measured velocities are, 
therefore, probably fairly representative for the motion of matter 
close to the equatorial plane of the nebula. 

The present theory favors the opinion that the direction of rota- 
tion of the system agrees with the direction in which the spiral arms 
wind around the center when followed outward from the central 
mass. In order to test this conclusion, we need, in the first place, 
observations of spectrographic rotations for typical spiral nebulae 
for which the true inclination toward the line of sight can be in- 
ferred with some degree of probability. An independent method 
based on the formation of secondary branches of spiral arms will be 


™ Ark, mat., astr. och fys., 27A, No. 2; Stockholm Medd., No. 45, 1939. 
M.N., 99, 650, 1939. 
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Fic. 5.—The photographic isophotes of the Andromeda nebula according to Lind- } 
blad and Ramberg. 
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mentioned below. When trying to solve the question of orientation 
of tilted nebulae from the appearance of the absorption lanes due to 
light absorption in dark matter, it is reasonable to begin with the 
great Andromeda nebula. A detailed photometric analysis of the 
nebula is in progress at the Stockholm Observatory, and a prelimi- 
nary picture of the run of the isophotes, drawn as straight lines be- 
tween different radii for which the intensity distribution has been 
measured, is given in Figure 5. The isophotes of individual dark 
clouds are not indicated in the figure, except in the southern (lower) 
part, where some isophotes which have been drawn as broken lines 
surround an intensity minimum due to certain widespread dark 
clouds. On the north-preceding side, however, when going outward 
from the center, the influence of the clouds on the general intensity 
distribution is especially visible in a sudden crowding of the iso- 
photes along certain lanes. From the general appearance of the ab- 
sorption, which, on the whole, proceeds from intense dark lanes step- 
wise in the direction away from the nucleus, I have drawn the con- 
clusion that the part of the nebula which shows heavy obscuration is 
in all probability the one farther from us."? The obscuring clouds 
have probably considerable extension at right angles to the equa- 
torial plane, and they may also be asymmetrically distributed in the 
angular co-ordinate around the center. An investigation by Ohman" 
on the polarization of light in the Andromeda nebula has strength- 
ened the conclusion as to the orientation. He has found that the 
dark clouds near the minor axis on the north-preceding side of the 
nebula are slightly deeper and more marked in the polarized com- 
ponent, where the plane of vibration is perpendicular to the plane 
passing through the center of the nebula and the observer, than in 
the other component. This would mean a negative polarization of 
the light scattered in the clouds, which it would be possible to explain 
if the angle nucleus-cloud-observer is far from 180°. In this connec- 
tion, the possibility of a different degree of illumination of the clouds 
on the near and the far side of the nucleus, respectively, on account 
of a different intensity of scattered light in the direction of the ob- 

"3 Ark. mat., astr. och fys., 24A, No. 21; Stockholm Medd., No. 14, 1934; Ann. d’ap., 1, 
173, 1938. 
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server, should be mentioned. As we must have a considerable flux 
of radiation from the central region outward, the angle center-cloud- 
observer should be of great importance for the intensity of diffused 
radiation if the diameters of the particles are of the order of the 
wave length of light. In general, for such particles the intensity of 
scattered light is much greater in a direction close to that of the in- 
cident beam than in directions forming an angle exceeding go° with 
the direction of the incident light. Under such conditions we must 
expect clouds at angles close to 180°—i.e., clouds on the near side — 
to be more strongly illuminated, and hence less visible by contrast, 
than clouds for which the angle is small. If, in accordance with these 
indications, the north-preceding side of the nebula is the farther one 
from us, then the direction of rotation found by Slipher and by Pease’ 
will be the direction in which the spiral arms wind outward, in ac- 
cordance with our theoretical rule. 

The difficulty of applying a criterion derived directly from the 
“spindle’’ nebulae is illustrated by M 81. Plate I gives two repro- 
ductions, with the same orientation, of a photograph of this nebula 
taken at the Lick Observatory, one showing only the nuclear region, 
which is ordinarily much overexposed in prints of the nebula, the 
other showing also the surrounding regions. In the first-mentioned 
picture there is a characteristic lane below and to the left of the nu- 
clear region; in the second picture the lanes are stronger in the upper 
right-hand part. If the inner lane close to the nucleus is situated on 
the near side of the nebula, in analogy with the lanes of the spindles, 
then the orientation of the nebula is such that the rotation and the 
direction of the arms agree with the present theory, according to a 
determination of the direction of rotation by Babcock. Certain 
other details in the distribution of bright and dark matter in the 
nebula, especially the appearance and position of certain neatly de- 
fined dark lanes along the spiral arms, are strongly in favor of this 
conclusion. In such a case, however, the stronger appearance of the 
lanes in the right-hand part of Plate I on the farther side of the 
nebula, if it is not due to an actual dissymmetry of the distribution 
of dark matter in the angular co-ordinate, might be caused by an 
actual elevation of the dark matter in these regions of the central 


1s Proc. Nat. Acad., 4, 21; Mt. W. Comm., No. 51, 1918. 
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PLATE II 


DIFFERENT TyPEs OF TILTED NEBULAE 


df 


a-c, NGC 4504, 3003, 7817; d-f, NGC 5746, 1337, 6503; g-i, NGC 4725, 4192, 972; 
j-l, NGC 2841, 3675, 5078. Lick Observatory photographs, except NGC 4725, which 
is a Mount Wilson photograph. 
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system over the most intensely radiating layers close to the equa- 
torial plane or may be an effect of a different illumination of the dark 
matter on the near and the far side of the central system. A combi- 
nation of all these effects, in analogy with what has been said con- 
cerning a corresponding effect in the Andromeda nebula, is perhaps 
the most likely explanation. 

An attempt to form a system of tilted nebulae of varying types is 
given here in Plate II. The nebulae are all orientated in the pic- 
ture in such a way that the edge which is supposed to be the near one 
is turned downward (right or left), the farther edge upward. The 
three successive columns are intended to show different degrees of 
the development of central mass. The left-hand column shows early 
types with large amorphous central parts, the right-hand column 
shows nebulae of very small central mass, and the column in the 
middle shows an intermediate type. 

The first nebula in the middle column shows unmistakably a lane 
on the near side analogous to that of the two real spindles in the first 
column. The central mass is small, and it is evident that the appar- 
ent disk has its maximum intensity close to the lane on the near side 
of the nucleus. 

The second and third nebulae in the middle column are obviously 
related in type to the uppermost one. The third nebula in this col- 
umn shows a very broad dark lane along the lower edge of the cen- 
tral system, though the faintly luminous lower border of the lane does 
not show well in the reproduction. A similar very large lane is faintly 
visible in the third nebula of the first column, in which the central 
mass is very much larger and where the spiral arm extends from a 
ring-shaped formation. 

The third column contains nebulae of very small central mass. 
In the nebula at the top of this column the nucleus is, in a character- 
istic way, displaced toward the lower edge of the apparent ellipse 
which is formed by the inner spiral arms, and the true orientation of 
the object is scarcely doubtful. The analogy between this nebula 
and the one immediately below it is very great, though the inner spi- 
ral structure is not so well developed in the last-mentioned object. 
That the inner structure in this case is not an amorphous central 
mass analogous to that of the early types in the first column is 
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shown by its richness of detail, and the real nuclear region is small. 
The nebula is classed by Hubble as Sc. The true orientation of the 
two lower nebulae in this column is rather obvious. On the whole, 
the sequence formed by the nebulae in the last column and the first 
nebula in the second column, taken together, is rather strong. The 
fact that in all these nebulae the light-intensity in the apparent el- 
lipse formed by the inner spiral arms is a maximum at a lane similar 
to that of the spindles and presumably situated on the near side, 
will be most readily understood if we assume the degree of flattening 
of the luminous matter in the central region to be so great that a con- 


Fic. 6.—Schematic vertical section of luminous and dark matter in a nebula of 
small central mass. Intersections with density surfaces of luminous matter are full 
drawn, intersections with dark matter are dotted curves. The oblique lines indicate 
the line of sight. 


siderable amount of dark matter actually may be considered as lying 
above the luminous matter. By an effect somewhat analogous to the 
darkening of the limb on the solar disk we then get a darkening to- 
ward the farther limb of the nebular disk. 

We now proceed to the last nebulae of the two first columns. By 
analogy with the neighboring nebulae their pr »bable orientation has 
been derived. The nebula in the middle column shows the darkening 
toward the farther limb extremely strongly developed. The general 
explanation is analogous to that just given and is further illustrated 
by Figure 6, where the luminous matter is schema‘ically represented 
by full-drawn curves, dark matter by dotted curves. Though the 
dark matter may have a very strong condensation toward the equa- 
torial plane, it is very probable that widespread and thin clouds ex- 
tend in an irregular way to a considerable height over this plane, 
forming, more or less, an envelope around the central regions. The 
effect of a longer path through this dark matter on the farther side 
may then produce a very effective darkening of the far half of the 
nebula. It is further conceivable, however, that the greater intensity 
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of the near side may be due partly to the diffusion of light emitted 
from the intense central region in the surrounding dark clouds. If the 
scattering of light by a particle of the dark matter is greater in the 
direction of the incident radiation than in the opposite direction, the 
dark clouds on the near side of the nebula are likely to appear illumi- 
nated by scattered radiation to a higher degree than the dark clouds 
on the farther side. Accurate measurements of the polarization of 
light may throw light on the extent to which such an effect may con- 
tribute to the apparent asymmetry of the light-distribution. 

A spectrographic determination of the rotation has been made 
by N. U. Mayall at the Lick Observatory for the second nebula in 
the last column of Plate Il, NG 6503. With the orientation in 
space suggested here the rotation is in the direction of the spiral 
whorls when proceeding outwaid, in agreement with the present 
theory. The nebula M 7331, for which the spectrographic rotation 
has been measured by H. W. Babcock, resembles, at least in its cen- 
tral regions, the nebula NGC 3675, which is the last nebula of the 
second column in Plate II. If NGC 7331 conforms to the ideas 
expressed here concerning the orientation, the observed rotation is 
in agreement with the present theory for this nebula as well. A 
photometric study of NGC 6503 and 7331 in various colors based 
on plates taken with the Crossley reflector of the Lick Observatory 
is at present in progress, and it is hoped that the results will help to 
interpret further the distribution of light in nebulae of these types. 

In the case of nebulaé which are so much inclined that we cannot 
distinguish the spiral arms sufficiently to determine which way they 
wind around the nucleus, we can perhaps infer this direction indi- 
rectly by the following way of reasoning. If we go in the direction 
of the major extension of the nebula to both sides from the nucleus, 
the dark matter may be expected to have the strongest obscuring 
effect on the concave side of the spiral arms, as seen by us. This 
follows if we assume that the obscuring matter of some elevation 
above the equatorial plane and close to the arms is more abundant 
inside the arms than outside of them. The half of the nebula which 
is relatively faint is then likely to be the one which turns the con- 
cave side of the arms in our direction. This criterion, applied to 
NGC 4594—for which the spectrographic rotation has been deter- 
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mined by Slipher'® and by Pease'?—appears to give a result in ac- 
cordance with the present theory. 

A phenomenon of a quite different kind, which may contribute to 
throw light on the dynamics of the nebulae, is the way in which a 
spiral arm often breaks up into a fan-shaped structure of “second- 
ary’’ branches. This phenomenon is very well shown, for instance. in 
M tor and M 33. Neighboring parts of the secondary arms are like- 
ly to possess about the same values of the area constant in the motion 
around the center. In such a case we should expect the branches ex- 
tending to the greatest distances from the center to recede relative 
to the branches of smaller distance from the center in a direction op- 
posite to the direction of rotation. This criterion, applied, for in- 
stance, to the fan-shaped spiral arm of M 1o1, evidently shows the 
direction of rotation to be the direction of winding of the arms, when 
we go outward. Measurements of radial velocity for different points 
of such fan-shaped spiral arms in nebulae of some inclination against 
the line of sight might contribute to test the assumption just made 
concerning the value of the area constant. 

STOCKHOLM OBSERVATORY 
SALTSJOBADEN, SWEDEN 
October 1939 
© Pop. Astr. 23, 24, 1915. 
17 Proc. Nat. Acad., 2, 517; Mt. W. Comm., No. 32, 1916. 


THE OPACITY OF GAS MIXTURES IN 
STELLAR INTERIORS 


PHILIP M. MORSE 


ABSTRACT 

The opacity of mixtures of hydrogen, helium, oxygen, sodium, magnesium, silicon, 
potassium, calcium, and iron is computed for any proportions of hydrogen and helium, 
and for four different mixtures of the heavier elements, for temperatures from 4X 104 
to 109 degrees centigrade, and for densities from 107? to 104 gm/cm3. Absorption cross- 
sections for Coulomb-field wave functions are used; but these are corrected for the 
screening of filled shells by using effective nuclear charges and quantum numbers, and 
for the effect of pressure ionization (excluded volume) by the use of the solutions of the 
Thomas-Fermi equation for electronic distribution in dense materials. A table of cor- 
rection factors for Kramers’ opacity formula (guillotine factors) is given for the four 
mixtures, for the ranges of temperature and density covered. Instructions for determin- 
ing the value of the opacity, for a given temperature and density and for intermediate 
chemical mixtures, are given in the Summary. There are comments on the dependence 
of the opacity on the chemical constitution of the mixture, and on the effect of these 
results on the calculation of hydrogen concentrations in stars. 

With increasing knowledge about the processes of energy genera- 
tion in stellar interiors it has become important to have accurate and 
extensive tables of the opacity of stellar material. B. Strémgren'’ has 
made calculations for a particular mixture of heavy elements (the 
Russell mixture) and for a limited range of temperature and density. 
It is the purpose of this paper to extend and improve Strémgren’s 
results, covering a wider range of temperature and density, using 
several different mixtures of heavy elements, and utilizing improved 
quantum mechanical techniques to obtain a more accurate result. 

The three processes that are important in opacity are the photo- 
electric absorption (bound-free transitions), collision absorption 
(free-free transitions), and scattering. Stobbe, Gaunt, Menzel and 
Pekeris,? and others have computed the probability of the absorp- 
tion processes for hydrogenic wave functions. Their results can be 


expressed in terms of effective cross-sections for absorption of light 


*Zs. f. Ap., 4, 118, 1932. See also Chandrasekhar, An Introduction to the Study of 
Stellar Structure, chap. vii, University of Chicago Press, 1939. In the present paper the 
notation of Chandrasekhar is followed as closely as possible. 

2M. Stobbe, Ann. d. Phys., 7, 661, 1930; J. A. Gaunt, Phil. Trans., A, 229, 163, 
1930; D. H. Menzel and C. H. Pekeris, M.N., 96, 77, 1035. 
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of frequency w atomic units’ from a filled shell of total quantum 
number or from an energy band d(k?) of free states (k? in atomic 
energy units), for an atom of nuclear charge Ze, 


3V3 me (Ry) gi(w) bem 


24¢? 
dQ. (w) = g(k)d(k?) cm?. 


(1) 


The factor (e?/mc) = 8.401073 cm?/sec, (Ry) is the Rydberg con- 
stant 3.29X10'5 sec"; and g is the so-called Gaunt factor, which 
varies slowly with w and u or k and has been shown "? to have the 
order of magnitude of unity. The energy W, is the absorption limit 
in atomic units; its value will be given in equations (8) and (12). 

To obtain the absorption coefficient per gram of material one 
must multiply these cross-sections by the probability of an electron’s 
being in the mth shell or in the range d(k?) (only the Fermi factor 
needs to be used, since the statistical weight of the state is already 
included in the Q’s), by the probability that the final state is unoc- 
cupied, and by the number of atoms Z per gram of material; and, 
finally, one must sum the products over all initial states and over all 
values of Z. 

At high temperatures this absorption opacity becomes negligible, 
compared with the opacity due to scattering. In these cases the 
Klein-Nishina modification of the Thomson scattering formula can 
be used, and the scattering cross-section per electron can be written 


= 6.57 X 10°*N(w) cm’, (2) 


3 Many of the intermediate formulae in the calculations are expressed in simpler 
form in atomic units: 


Symbol Quantity One Atomic Unit Equals: 
Energy =13.53 ev =2.152 ergs 
Frequency = 3.29 X 10'S sec~! = Rydberg constant 


(wave length of light of frequency w atomic 
units equals g00/w angstrom units) 


| Seer Density 11.2 gm/cm (one hydrogen atom per unit atomic 
volume corresponds to unit density in atomic 
units) 

Temperature | =157,000° C (i.e., the corresponding 


Boltzmann factor is unity) 
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where values of the function N are given in Table 1. The function 
has the following limiting behavior: 


—— (1000 > w) 
8, 
(3) 


— (105 < w) 
w 11,400 


The scattering coefficient is the product of Q, and the total number 
of electrons per gram. The total opacity is the sum of the absorp- 
tion and scattering coefficients and is a function of the light fre- 
quency w, of the chemical composition, and of the temperature. 


THE EXCLUDED VOLUME 


Equations (1) have been derived for hydrogenic wave functions 
and for exceedingly small gas densities. In stellar interiors each 
electronic orbit is distorted by the presence of other electrons in or- 
bits about the same nucleus, and also by the presence of other atoms 
near by. Equations (1) should be corrected for both of these effects. 
The correction due to the presence of other electrons about the same 
nucleus can be made by using the effective charge and the effective 
quantum number, instead of actual Z and n, in the equations. This 
will be discussed later. 

The effect of the presence of neighboring atoms can be computed 
by the help of the quantum theory of metals. According to this 
theory, the atomic energy states are modified as the atoms are 
squeezed closer and closer together to form a dense material. The 
bound states have sharp, discrete levels in an isolated atom; but, as 
the density is increased, the levels spread out into bands. The band 
width is greater the higher the level, so that, above a certain band, 
all the bands will overlap, and the upper states will be more like free 
states than bound ones. In other words, for any given density there 
will be a lower limit of energy, above which the states, originally 
bound, are now virtually free states and below which the states are 
still tightly bound with very little spreading-out of the energy bands. 
This depression of the lower limit of the free electron states increases 
with increasing density, and for high-enough densities the bound 
states will all be absorbed by the free states. The effect has been 
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TABLE 1* 


VALUES OF THE SCATTERING FUNCTION, THE NORMALIZING CONSTANT 
FOR FERMI STATISTICS, AND THE STROMGREN FUNCTION S 


logw N(w) logn log A Fo S(u) u S(u) 
(Eq. [18]) (Eq. [6]) (Eq. [16]) (Eq. [15]) 
2.0 0.995 2.0 2.0015 0.9981 ©.0 0.00000 7.5 94.0856 
2.2  .ggl 2 2.2024 .9969 0.5 0.00029 7.6 97.5219 
2.4 .985 4 2.4038 .9952 1.0 0.01071 7.7 100.3346 
2.6 .976 6 2.6060 9931 1.5 0.08954 103.3223 
2.8 .963 8 2.8096 9902 2.0 0.39955 7.9 105.8825 
2.5 1.2434 
3:0 045 f:0 . .9857 8.0 108.6119 
.936 .9825 3.0 3.0480 8.1 
2.2 .922 2 .9785 3.2 4.1483 8.2 113.9608 
3-3 -904 3 1.3303 .9733 3:4 5.5041 8.3 116.5707 
3.4 4 1.4382 .9670 3.6 7.1469 8.4 119.1587 
3-5  .856 5 1.5480 .9592 3.8 9.0057 8.5 121.6939 
3.6 .826 6 1.6605 .9499 4.0 11.3676 8.6 124.1825 
2.9 29092 7 1.7760 .9385 4.2 13.9808 8.7 126.6251 
236 2955 8 1.8956 .9245 4.4 16.9409 8.8 129.0210 
2:9 0.0201 .9082 4.6 20.2476 8.9 131.3691 
4.8 23.8940 
4:0 672 0.0 0.1511 .8880 9.0 133.6682 
1 0.2898  .8640 5.0 27.8664 Q.2 138.1127 
4.2 .583 2 0.4380 .8370 5.1 29.9725 9.4 142.3489 
4:3 3. 0.5982 8061 5.2 32.1510 9.6 146.3742 
4.4  .498 4 0.7730 .7726 5-3 34.4014 9.8 150.1856 
4-5 -457 5 0.9654 .7367 5.4 36.7225 10.0 153.7856 
4.6 .417 6 1.1804  .6993 5.5 39.1128 10.2 157.1759 
4.97. 2378 7 1.4215 .6601 5.6 41.5698 10.4 160.3624 
4.8  .340 8 1.6959 .6187 5-7 44.0895 10.6 163.3490 
4:9 «803 Q 2.0090 5751 5.8 46.6656 10.8 166.1404 
5-9 49.2017 11.0 168.7424 
5.0 .268 1.0 2.3699 5281 
me 2.7816 .4746 6.0 51.9633 174.4525 
5.2 . 205 2 3.2672 .4123 6.1 54.6755 12.0 179.1535 
3 3.8191 6.2 57.4249 12.5 182.9749 
4 4.4628 .2772 6.3 60.2082 13.0 186.0382 
5 5-2193 2213 6.4 63.0213 13.5 188.4766 
5.6 .115 6 6.0944  .1758 6.5 65.8637 14.0 190.3846 
5-7 -0©99 43307 6.6 68.7225 15.0 193.0267 
5.8 .085 8 8.3010 .1109 6.7 71.5978 16.0 194.5800 
5.9 .073 9g 9.6872 o881 6.8 74.4853 18.0 195.9583 
6.9 77.3810 20.0 196.3685 
6.0 .063 2.0 I1.299 .0700 
6.2 .046 233576 .0556 7.0 80.2800 22.0 196.4813 
6.4 .033 .2 15.304 .0442 7.1 83.1773 24.0 196.5103 
6.6 023 -3 “27.055 0351 7.2 86.0698 20.0 196.5174 
6.8 -4 20.892 0279 7.3 88.9546 30.0 196.5194 
7.0 0.009 -§ 24.364 0.0222 7.4 91.8283 © 196.5194 


* The symbol “‘log’’ characterizes the logarithm to the base 10, whereas ‘‘In’’ in text characterizes the 
logarithm to the base e. 
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called ‘‘pressure ionization.”” The amount of the depression can be 
computed approximately by means of the Thomas-Fermi equation.‘ 

The usual Thomas-Fermi solution gives the effective potential 
energy of an electron in an isolated atom when the atomic electrons 
are considered to form a degenerate gas. If the uppermost filled 
levels are placed at zero energy, then the potential is approximately 
equal to —(Ze?/r) for small values of r (the distance from the nu- 
cleus) and approaches zero asymptotically for large r. This potential 
is valid for an isolated atom. Slater and Krutter’ have computed the 
Thomas-Fermi potential when the atomic electrons are confined in- 
side a sphere of radius a,, corresponding to the case of an atom in a 
compressed gas. In this case the potential is again equal to — (Ze?/r) 
for small 7; but it rises to a maximum value, less than zero, at r = a. 
The upper limit of the filled levels is above zero, since ‘he squeezing 
has increased the average kinetic energy of the electrons. It turns 
out® that the position of this potential maximum coincides fairly 
closely with the position of the lower limit of the free states men- 
tioned above. To a fairly good approximation the atomic states 
with levels below the maximum can be considered to be bound states, 
unchanged by the presence of the near-by atoms, whereas the levels 
that were above this limit in the isolated atom are effectively merged 
in the continuous states for free electrons. Figure 1 shows the poten- 
tial energy-curves for a given density for three elements and shows 
the lower limit of the continuum and the bound levels below this 
limit. 

The problem of joining the effective potential-curves about each 
nucleus to the curves for the next nuclei becomes somewhat compli- 
cated when the gas is a mixture of different kinds of atoms. One 
must choose a radius for each nucleus such that the energy difference 
between the potential maxima at the edge of each atomic sphere and 
the top of the filled levels is the same for every nucleus. When this 
is true, the top of the filled levels is the same throughout the gas, and 
the potential energy in the interspaces between atoms is nearly con- 
stant. Slater and Krutter’s tables® enable one to compute the ratio 

4L. H. Thomas, Proe. Camb. Phil. Soc., 23, 542, 1927; E. Fermi, Zs. f. Phys., 48, 


73, 1928. 
5 Phys. Rev., 47, 559, 1935- ® Rev. Mod. Phys., 6, 210, 1934. 
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of the atomic radius a, for atom Z in terms of the radius a, for hydro- 
gen. These ratios are approximately proportional to the cube root of 
the nuclear charge; (a./a,) © (Z/Z’)"/3, for most of the values of a, 
that are important in the present problem. It is possible, therefore, 
to find an approximate value for a, as a function of density for any 
mixture of elements. 

Let U, be’ the proportion by weight of the element of nuclear 
charge Ze and atomic weight M.. Then the number of such nuclei 


Continuum 


Fic. 1.—Effective electronic potential as a function of distance from three nuclei, as 
computed by the Thomas-Fermi equation, for p = 100 gm/cm3 and for X = 0.16. 
The dashed lines above the bottom of the continuum show the asymptotic value of the 
potential of the corresponding isolated atom. The dashed lines below the bottom of the 
continuum show the positions of bound states not merged in the continuum. 


per unit atomic volume is (U.p,/M.), where p, is the density of the 
gas in atomic units.’ The volume occupied by such an atom is 
(47a3/3); and, since a3 is approximately equal to a}Z, the total volume 
occupied by all the atoms in the unit volume must be 


7 In later discussions in this paper we shall find it more useful to let the proportions 
by weight of hydrogen and helium in the total mixture be UV; = X and U, = Y, respec- 
tively, and the proportions of the heavier elements be Uz = (1 — X — YV)cz, so that cz 
is the proportion by weight of atom Z in the heavy mixture before it is diluted with 
hydrogen and helium. The use of X and Y follows the notation of Chandrasekhar (n. 1, 


above). 


| 
Hf 
1 t LLLLS /, 
Electron 
Potential ase 
Fe oO 
470 }Pa U.Z 
3 <M." 


OPACITY OF GAS MIXTURES 33 


which must equal unity. For hydrogen the ratio (Z/M.) = 1, and 
for all other elements (Z/M,) 2 3. Therefore, the summation re- 


duces’ to 3(1+.X), since = U, = 1. The approximate value of a, 
is thus 


where (p,/2) (1+X) is the number of electrons per unit atomic vol- 
ume. 


Fc. 2.--Depression of continuum £; in atomic units as a function of the number of 
electrons per unit atomic volume, pa(1+X), for three elements. 


Once a, is known for a given p, and_X, it is possible to go back to 
the tables of Slater and Krutter to obtain more accurate values for 
a. for the various nuclei in the mixture, and also to calculate E., the 
energy difference (in atomic units) between the potential maximum 
for the computed a, and the asymptotic value of the potential for 
the same atom when isolated. The energy £, is, therefore, the 
amount by which the lower limit of the continuous free states has 
been depressed for this nucleus by the squeezing-together of the 
atoms. The values of £, in atomic units for three different elements 
are given as functions of p,(1+X) in Figure 2. 

These calculations have been made for zero temperature, for a 
completely degenerate gas; and a correction should be made for the 
fact that the gas under consideration is hot. It turns out, however, 
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that whenever absorption by the bound states is important this cor- 
rection has a negligible effect on the computed opacity, so that the 
correction has been neglected. 


THE FERMI DISTRIBUTION 


Since the electrons follow Fermi-Dirac statistics, the probability 
of an electron’s being in a given state of energy W is 


P(W) = +4 (5) 


where W and 7, are both in atomic units and where the constant A 
is determined by the requirement that the sum of the P’s over all 
states per unit volume equals the number of electrons per unit vol- 
ume (p,/2) (1+X). The zero of energy is taken to be the lower 
limit of the free states. At the temperatures and densities occurring 
in stellar interiors the great majority of the electrons are in free 
states, either because the temperature has forced them out of the 
bound states or because pressure ionization has squeezed them out. 
Therefore, A is determined primarily by the distribution of free 
states. 

The number of free states per unit volume having momentum in 
atomic units between k and k + dk is (k?/1?)dk. The energy of these 
states is k?. Therefore, an approximate equation for A is 


= 2p.(1 + X) (: = . (6) 


A table of values of log A as a function of log 7 is given* in Table 1. 
The solutions for A for limiting ranges of are 


exp[(1.21)*/3(1 — 10) 


8 The quantity 7 corresponds to NV-/G(T7) in the notation of Chandrasekhar (n. 1, 
above). 
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The values of A for intermediate values of 7 must be obtained from 
Table 1. 

Once values of A are obtained, it is possible to correct formulae (1) 
for the effect of other electrons inside the atom. For a given density 
and temperature one computes for each atom the number of shells 
which are filled (i.e., those for which P is close to unity) and applies 
the screening constants given by Slater? to obtain an effective 
charge Z, to use instead of Z in equations (1). Sometimes these cor- 
rections alter the number of shells shown to be filled, and the calcula- 
tions must be remade until a self-consistent result is obtained. The 
values for W, in equations (1) are then given by 


Wir = (2) — E, (see also eq. [12]) , (8) 


where E, is the depression of the lower limit of the free states given 
in Figure 2. Bethe’® has shown that formulae (1) still hold when 
screening and pressure ionization occur, as long as one considers 
states with W,,, > o to be the only bound states. 

After all these corrections have been made, one can form a general 
picture of the relevant phenomena that occur in the gas as the pres- 
sure and temperature are varied. Figure 3 shows a general sketch 
of these occurrences in the p, T plane (cgs units). To the right, for 
high temperatures, the scattering opacity predominates; and to the 
left, the absorption does. As we reduce the temperature, the various 
shells fill up, first the K shell of iron (the heaviest atom considered), 
and so on, until at about 500,000° C so many shells are filled that 
equations (1) begin to be inadequate even with the corrections made 
for screening given above. Consequently, no results are given in this 
paper for temperatures below 400,000 °. 

As the density is increased, one shell after another is merged in 
the continuum of free states, until finally even the K shell of iron dis- 
appears. By this time the free states have become degenerate and 
electronic conductivity begins to short-circuit the opacity, as is 


9 Phys. Rev., 36, 57, 1930. 
© Handbuch der Physik (2d ed.), 24, Part I, 478, Berlin: Springer, 1933. 
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shown by Marshak." Consequently, no results are given in this 
paper for densities larger than 104 gm/cc. 


THE ABSORPTION AND SCATTERING OPACITIES 
We are now in a position to calculate the opacity as a function of 
frequency. The opacity due to scattering may be obtained from 


Continuum Degenerate 


hell 
4/ 


KS! 
11, 


But 
Scattering” 
= Jer 
= = reporderant, 
~ 
oT oS 
o~ ~~ 
FAS 


=, 
-2 
5 7 & 9 
log T 


Fic. 3.—Sketch of p, T plane, showing phenomena important in opacity 


equation (2), by multiplying by the number of electrons per unit 
mass of material, 3.03 X 107} per gram, giving 


ks(w) = 0.20(1 + ¥)N(w) cm?/gm. (9) 


The opacity due to free-free transitions is somewhat more com- 
plicated. We multiply dQ; by the number of Z atoms per gram and 
sum over Z. Since there is a strong dependence of Y on Z, we may 
neglect the effect for hydrogen and helium and sum over only the 
heavier elements. There results a factor T = Dd c(Z*/ M.), where 
c. is the proportion by weight’ of atom Z in the mixture of heavy ele- 
ments, and I is the average value of (Z?/M.) over this mixture. I 


may be called the “composition factor” of the mixture. 
It is also necessary to multiply by P(X’) to give the probability of 


" Ap. J.,in press. The writer is indebted to Dr. Marshak for the opportunity to 


read his paper in advance of publication. 
| 
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the electron’s being in initial state k?, and by [1 — P(k?+w)] to give 
the probability of the final state being vacant. Integration over k 
gives the absorption for free-free transitions, 


oo ) 
= 4.76 X 10°(1 — X — 


[(1 + (1 + dx 


= 4.76 X 10%(1 — X — (1 — (10) 


x in| cm /gm 
where x = k?/T,, y = w/T., x = InA, and g is the average value 
of the Gaunt factor. Values of A can be obtained from Table r. 
To obtain the total absorption opacity we must add to this the 
absorption due to bound-free transitions. From equations (1), by 


the same method used for x;, we obtain the total absorption opacity, 


Ka = 4.76 X 10% 3g(1 — X — V) — 


(11) 


where the summation is for all states with values of W.n (see eq. [12]) 
smaller than w. The most important term in x, is usually the term 
with largest W.,, and no large error is made in making g equal to 
the Gaunt factor for this term. Values can be obtained from the 
work of Menzel and Pekeris.2 In addition, a correction must be 
made at high densities for the fact that the integration for the transi- 
tion probabilities is made over a sphere of radius a, instead of over all 
space. This does not make so large a correction as might be expected, 
since the bound-free transitions involved are always those with wave 
functions that have small values beyond a.. 

Both these corrections can be combined into an effective Gaunt 
factor g. = g(w)ga, where g(w) is the Gaunt factor for the largest 
term in the summation of equation (11), and g, is the term involving 
the density correction. Table 2 gives approximate values of these g 


| | 
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factors. Instead of g(w), there is given the factor g,, the average of 
g(w) over the Planck distribution, as discussed in the next section. 

The factor [1 — P(w—W.,)] = [1 + Aexp(1/T.) (Win—w)]™ is 
nearly equal to unity whenever A is smaller than unity. For large 
densities and low temperatures, however, A is much larger than 
unity, and the expression above is very small for small values of w, 


TABLE 2* 


CORRECTION FACTORS FOR GAUNT FACTOR, g=g,g,; INTERPOLATION 
FACTORS TO OBTAIN MEAN OPACITY FROM ABSORPTION AND 
SCATTERING OPACITIES 


1.0 0.00 5.6 0.00 3.0 0.00 0.0 0.41 
.00 5.8 .00 0.50 
1.4 .O1 .4 .4 
.O1 6.0 .02 °.79 
1.8 .02 6.2 8 8 0.95 
6.4 .06 
2.0 .03 6.6 .08 02 1.0 
2.2 .O4 6.8 mere) .03 2 1.28 
258 .06 4 .O4 4 1.46 
2.6 .08 7:0 10 6 6 1.64 
2.8 .10 8 .O7 8 1.83 
12 
3.0 7.6 12 2:0 2.02 
$02 .16 7.8 2 2 2.21 
3-4 19 + 15 + 2.40 
3.6 23 8.0 12 6 20 6 2.60 
3.8 8.2 12 8 8 2.80 
8.4 II 
4.0 0.30 8.6 0.10 0.0 3.0 3.00 


* Logarithms are to base ro. 


rising rapidly to unity when w becomes larger than W., + 7.InA. 
Therefore, one further approximation has been made in equation 
(11). The factor [1 — P(w—W.,)] was omitted, and a new absorp- 
tion limit was used, whenever degeneracy is present (i.e., whenever 
InA is positive). The effective limit used is 


Z.\? (log A < 0) \ 
n ad log A (log A > 0) 


where E, is given in Figure 2, and log A in Table 1. 


| 

| 

| 
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To obtain the opacity in a form for easy comparison with Kra- 
mers’ formula, we multiply and divide by 7 = 11.14pa(1+X)T73”; 
and, to make calculation easier, we arrange the terms in the summa- 
tion in equation (11) in order of increasing value of W., rather than 
in order of increasing 2 or Z. The smallest value of W., will be re- 
labeled W,, and so on, so that the sth term in size is W,, where W,_,; 
< W, < W,,:. Then the total absorption opacity becomes 


Ka(w) 


5.30 X + X)(1 — X — v) 
X [F. + Fi(w)], Ws <w < Wess, 


I 
+4 (13) 
= 


s 
Lilo 
Monon 
c= 


F, P(W.e) . 


4 


The quantities Z,, ,, etc., are the values of Z, m, etc., corresponding 
to the oth absorption limit in the ordered sequence of W.n. The 
quantity F, is, therefore, a step function, changing value only when 
w passes through an absorption limit W,. 


THE AVERAGE OPACITY 


The important quantity in stellar calculations is the Rosseland 
mean of the opacity 


k= — |” (u = (14) 


Since the reciprocal of the opacity is averaged, it is not possible to 
obtain the average total opacity by adding the average absorption 
opacity and the average scattering opacity. For those temperatures 
where both opacities have the same magnitude, formulae must be 
devised to combine the two opacities, but for other temperature 
ranges the problem is fairly simple. 

In the low-temperature regions the scattering opacity x, is negli- 
gible compared with x, so that equation (13) is the expression to be 
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averaged. The method of calculation is that of Strémgren' and uses 
Stroémgren’s function 


S(u) = (vere — 1)3dx, (15) 


which is given in Table 1 for convenience. The function has been 
recalculated to obtain a somewhat greater accuracy than that of 
Strémgren’s original table. Intermediate values can be obtained by 
four-point interpolation. 

When x, predominates, the integrand of equation (14) has a fairly 
sharp maximum at u & 7, so that’ for slowly varying quantities we 
may replace w by 77, and place the quantities outside the integral. 
After the integration is performed, it is convenient to return to units 
of grams per cubic centimeter for density and degrees centigrade for 
temperature. The result is 


Ka © 7.23 X p(r + — X — cm’/gm, (16) 


where 
g = gag. (see Table 2) , I 
its 
— F, + F. > 100) | 


and where F, is given in equation (13). The function F, is given in 
Table 1 as a function of 7. 

The formula for x, without the factor (g/¢) is just Kramers’ formu- 
la, the correction factor ¢ being called the “guillotine” factor. The 
results below show that, in most cases, ¢ is larger than unity, indi- 
cating that the actual opacity is smaller than that predicted by the 
simple Kramers’ formula. Tables of (¢/Z) are given for four different 
mixtures of heavy elements in Table 3. Case 1 is for the Russell mix- 
ture. No values are given for high temperatures when x, becomes 
negligible with respect to k,. 


2 Tf the quantity Fs, in eq. (13) were independent of w, the maximum of the inte- 
grand would come at w = 77. However, the dependence of Fs on w makes this maxi- 
mum come nearer to 674. Therefore, for those parts of the opacity dependent on the 
free-free transitions we set wmax = 77a, and for those parts corresponding to bound- 
free transitions we set wmax = 607. 


TABLE 3 


VALUES OF THE LOGARITHM OF THE GUILLOTINE FACTOR LOGio (¢/g) 
(See Eq. [16]) 


CASE 1: CFhe=0.125, CK, Ca=CSi=0.0625, CNa, Mg =0.25, CO=0.5 


LOG + X) 

Loc 7 

—2.0 | —1.5| —1.0] —0.5| 0.0 0.5 1.0 ee 2.0 2.5 3.0 3.5 4.0 
5.6 0.10 | 0.28 | 0.56 | 0.88 | 1.23 | 1.62 | 2.02 | 2.46 | 2.92 | 3.43 | 3.08 | 4.55 | 5.14 
Ce | 0.09 | 0.17 | 0.44 | 0.76 | 1.09 | 1.48 | 1.90 | 2.34 | 2.80 | 3.30 | 3.83 | 4.41 | 5.00 
5.8 0.13 | 0.17 | 0.36 | 0.67 | 1.01 | 1.38 | 1.78 | 2.21 | 2.66 | 3.16 | 3.69 | 4.27 | 4.86 
5.9 0.23 | 0.27 | 0.38 | 0.65 | 0.98 | 2.30 | 2.68 | 2.09 | 2.54 | 3.03 | 3-50 | 4.33 | 4-72 
6.0 0.32 | 0.35 | 0.42 | 0.63 | 0.92 | 1.22 | 1.58 | 4.98 | 2.41 | 2.90 | 3.42 | 4.00 | 4.58 
6.4 0.2 0.30 | 0.35 | 0.55 | 0.81 | 1.12 | 1.47 | 1.85 | 2.27 | 2.76 | 3.28 | 3.86 | 4.44 
6.2 0.2r | 0.23 | 0.28 | 0.45 | 0.60 | 0.08 | 2.32 | 3.71 | 2-83 | 2.68 | 3-34 | 3.92 | 4.30 
6.3 0.20 | 0.22 | 0.26 | 0.37 | 0.56 | 0.84 | 1.18 | 1.57 | 2.00 | 2.47 | 3.00 | 3.58 | 4.16 
6.4 0.17 | 0.19 | 0.22 | 0.28 | 0.42 | 0.70 | 1.04 | 1.43 | 1.85 | 2.33 | 2.86 | 3.44 | 4.02 
6.5 0.11 | 0.13 | 0.15 | 0.20 | 0.32 | 0.58 | 0.92 | 1.28 | 1.71 | 2.19 | 2.72 | 3.30 | 3.88 
66 0.10 | 0.12 | 0.14 | 0.19 | 0.28 | 0.50 | 0.81 | 1.17 | 1.58 | 2.05 | 2.58 | 3.16 | 3.74 
6.7 0.13 | 0.15 | 0.17 | 0.21 | 0.28 | 0.45 | 0.72 | 2.07 | .47 | F.92 | 2-44 | 3.02 | 3.60 
68 0.22 | 0.23 | 0.25 | 0.27 | 0.33 | 0.47 | 0.70 | 1.01 | 1.38 | 4.82 | 2.32 | 2.88 | 3.46 
6.9 ©.3% | 0.32 | 0.34 | | 0.43 1 0. §3 | O.72 | 0.00. | 3-42 | 3.79 | 2-48 3-34 
7.0 0.41 | 0.42 | 0.45 | 0.49 | 0.54 | 0:62 | 0.77 | 3.00 | 2.30 | 1:67 | 2.12% | 2.63 | 3.26 
7.1 0.49 | 0.50] 052]0.54/0.5 | 0.77 | 0.908 | 3.25 | 3.6% | 2.02 | 2.52 | 3.07 
7.2 0.56 | 0.56 | 0.57 | 0.58 | 0.60 | 0.64 | 0.73 | 0.92 | 4.18 | 1.52 | 1.93 | 2.42 | 2.05 
7.3 0.63 | 0.63 | 0.64 | 0.65 | 0.67 | 0.69 | 0.74 | 0.87 | 1.06 | 1.41 | 1.83 | 2.32 | 2.84 
0.75 | | 6.95 0-7 o.77 | 0.79 | 0.83 | | 3.05 | 3.36 | | | 2-94 
778 © go | 0.90 | 0.g0 | 0.91 | 0.93 | 0.95 | 0.98 | 1.04 |] 1.16] 1.40] 1.75 |] 2.19 | 2.66 
1.33 | 1-33 | 2-33. | | | 3.396 | 2.38] B-42 | | 1.69 | 2.05 | 2.25 2-58 
7.8 1.59 | 1.50 2.590 | I-50 | 2.50 | 1-60 | | 3.66} 3.73 | 3.86 1 2:04 2.90) 2-55 
7.9 y.S2 | 1.8t. | | 1.82 | | 2.82 | 2.83 | | 2:00 ata 1 
80 2.00 | 2.00 | 2.00 | 2.00 | 2.00 | 2.00 | 2.01 | 2.03 | 2.06 | 2.12 | 2.24 | 2.41 | 2.62 
2.14 | 2-14 2.54 | 284] 2-84 | 2.26 | 1 2.40 2-65 
8 2 2.27 | 2.27 | 2.27 | 2.27 | 2.27 | 2.27 | 2.27 | 2.28 | 2.20 | 2.34 1 2.43 | 2.52 | 2.68 
8 3 2.36 | 2.36 | 2.36 | 2.36 | 2.36 | 2.36.1 2.36 12.36 | 2.37 | | 2.96% 2.32 
84 2.44 | 2.44 | 2.44 | 2.44 | 2.44 | 2 2.44 | 2.44 | 2.45 | 2.47 | 2.52 | 2.60 | 2.74 

CASE 2: CFe=0.3, CK, Ca=CSi=0.05, CNa, Mg =0.2, CO=0.4 
LOG p(1+X) 

Loc 

—20 | —1.5| —1.0| —0.5] 0.0 0.5 1.0 2.0 2.5 3.0 3.5 4.9 
5 6 $o.13} 0 31 | 0.59 | O.g1 | 1.25 | 2.64 | 2.03 | 2.46 | 2.92 | 3.43 | 3.07 | 4-55 | 5-14 
i +-o.09] 0.18 | 0.45 | 0.77 | £.18 | 1.50 | 2-98 | 2.33 | 2.79 | 3-20 | 3.83 | 4.42 | §.00 
5.8 +o og] 0.13 | 0 36 | 0.68 | 1.03 | 1.40 | 1.79 | 2.21 | 2.66 | 3.16 | 3.69 | 4.27 | 4.86 
5.9 +O.20] O. 2: 0.41 | 0.67 | 0.99 | 1.32 | 1.69 | 2.10 | 2.53 3.03 | 3.55 | 4. 4.72 
60 to 36] 0 40 | 0 48 | 0 69 | 0.97 | 1.26] 1.61 | t.99 | 2.41 | 2.90 | 3.42 | 4.00 |] 4.58 
61 +0. 35] 0.38 | 0.45 | 0.66 | 0.93 | 1.19 | 1.52 | 1.87 | 2.28 | 2.76 | 3.28 | 3.86 | 4.44 
6.2 to 26} 0 29 | 0. 36]055|081 1.09 | 1.41 2.76 | 2.35 2.62 3.14 3.72 | 4.30 
6.3 +o 16] 0.19 | 0.25 0.40] 0.63 0.94 1.28 1.64 2.03 2 48 3.00 3.58 4.16 
+007] 0.10 | 0.15 | 0.2 0.45 | 0.77 | | | 4.90 | 2.38 | 2-87 | $284 1 4-02 
65 —o ot] 0.02 | 0.05 | 0.14 | 0.30] 0.61 | 0.906 | 1.36 | 1.76 | 2.22 | 2.74 | 3.30 | 3.88 
6.6 +o 02] 0.02 | 0.04 | 0.10 | 0.23 | 0.48 | 0.82 | 1.21 | 1.63 | 2.08 | 2.59 | 3.16 | 3.74 
| 0.13 | 0.15 | 0.24 | 0.41 1.09 | 1.95 2.46 | 3.02 | 3.60 
68 +o. 24] 0.24 | ©. 2. 0.25 | 0.30 | 0.44 | 0.69 | 1.03 | 1.41 1.84 | 2.32 | 2.88 | 3.46 
69 +o 35] 0.36 | 0.38 | 0.40 | 0.43 | 0.52 | 0.71 1.o1 ¥.33 | 3.74 | 2.22 2.76 | 3.33 
7.0 +o. 44] 0.45 | 0.48 | 0.51 | 0.55 | 0.62 | 0.77 | 1.00 | 1.29 | 1.66 | 2.10 | 2.63 | 3.19 
Fae +0.50] 0.52 | 0.55 | 0.58 | 0.62 | 0.69 | 0.80 | 0.99 | 1.24 | 1.56 | 1.07 2.49 | 3.06 
7.2 $0.54] 0.56 | 0.59 | 0.63 | 0.66 | 0.70 | 0.77 | 0.93 | 1.15 | 1.44 | 1.84 | 2.36 | 2.93 
rae. +o. 61] 0.62 | 0.64 | 0.66 | 0.68 | 0.69 | 0.71 | 0 82 | 1.00 |] 1.30 | 1.72 | 2.23 | 2.79 
7.4 +0.741 0.73 | 0.72 | | 0.7 o.69: | 0.67 | 6.72 | 0.85 | | 1.66 | 4:32 | 4.6 
7.5 +o.90} 0.86 | 0.82 | 0.78 | 0.75 | 0.71 | 0.71 | 0.75 | 0.89 | 1.16 | 1.56 | 2.04 | 2.57 
7.6 +1 07] 1.03 | 0.98 | 0.92 | 0 88 | 0.86 | 0.88 | 0.93 | 1.05 | 1.30 | 1.62 | 2.03 | 2.50 
+1.30] 1.2 I. 2. ¥.20 | | 2.34 | 1.34 | 3.16 | ¥.20 | 2.50 | 3.74 | 2.07 | 2.40 
7.8 +1.52] 1.49 | 1.46 | 1.43 | 3.4% 3.40 3.5% 1.66 | 1.86 | 2.13 | 2.46 
7.9 +1.73| 1.71 | 1.71 | 1.70 | 1.69 | 1.67 | 1.67 | 1.69 | 1.73 | 1.82 | 1.98 | 2.20 | 2.48 
80 +1 go] 1.90 | 1.90 | 1.90 | 1.89 | 1.87 | 1.87 | 1.88 | r.90 | 1.96 | 2.09 | 2.28 | 2.51 
81 +2.04| 2.04 | 2.04 | 2.04 | 2.03 | 2.03 | 2.02 | 2.03 | 2.05 | 2.10 | 2.20 | 2.35 | 2-55 
+2.24] 2.2 2.24 | 2.24] 2.24.1 2.23) 9.23.) 2.24 | 2,25 | 2.28 2. 2203 
8.4 +2.32} 2.32 | 2.32 | 2.92 | 2.32 | 2.392 | 2.92 | 2.34 | 2.35 | 2:36 2.42) 2-66 
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TABLE 3—Continued 


CASE 3: CF. =0.2, CK, Ca=CSi=0.133, CNa, Mg =0.2, (9 =0.333 


LOG p(1+X) 

Loc T 

=2,0 | | —e.51 626 0.5 1.0 rs 2.0 205 3.0 3-5 4.0 
5.6.. 0.13 | 0.31 | 0.60 | 0.93 | 1.29 | 1.66 | 2.04 | 2.47 | 2 62 | 3.43.] 3:07 | 4.55 | s.3%4 
0.10 | 0.21 | 0.46 | 0.79 | 1.15 | 1.52 |] 1.92 | 2.34 | 2.79 | 3.20 1 
5.8. 0.13 | 0.17 | 0.37 | 0.69 | 1.04 | 1.41 1.82 | 2.22 | 2.67 | 3.16 | 3 69 | 4.27 | 4.86 
eae 0.22 | 0.27 | 0.39 | 0.66 | 0.98 | 1.32 1.70] 2.32 2.54 | 3.03 | 3-55 |] 4.13 | 4.72 
6.0: . 0.26 | 0.30 | 0.39 | 0.62 | 0.93 | 1 24 | 1.60 | 2.00 | 2.42 | 2.90 | 3.42 | 4.00] 4.5 
6.2.. 0.23 | 0.25 | 0.32 | 0.53 | 0.83 | 1.13 | 1.49 | 1.87 | 2.28 | 2.75 3.28.1 3.83 4-44 
6.2 0.19 | 0.20 | 0.25 | 0.42 | 0.68 | 0.99 | 1.35 | 1.73 | 2.13 | 2.60 | 3.14 | 3 92 |-4.30 
0.35: 0.16 | 0.17 | 0.19 | 0.31 | 0.52 | 0.85 | r.21 | 1 58 | 1.90 | 2.45 | 2.00 | 3.58 | 4.16 
6.4... O.T4 | 0.15 | 0.17 | 0.24 | 0.39 | 0.70 | 1.05 | 1.42 | 1.82 | 2.30 | 2.85 | 3 44 | 4.02 
6.5; O.12 | 0.13 | 0.15 | 0.20 | 0.31 | 0.56 | 0.91 | 1.28 | 1.67 | 2.15 | 2.72 | 3 30 | 3.88 
6.6 0.06 | 0.07 | 0.09] 0.14 |] 0.24 0.47 10.79 / 1.15 | 1.54] 2.01 2.560] 3.15 3-74 
6.7.. 0.06 | 0.07 | 0.09 |] 0.12 | 0.22 0.42 10.69] 1.04 1.41 I 2.41 3.01 3. 60 
6.8... 0.17 | 0.17 | 0.18 | 0.20] 0.26 | 0.41 | 0.64] 0 QO 2.33 2.86 | 3.45 
6.9... 0.20 | 0.26 | 0.28 | 0.32 | 0.36 | 0.43 | 0.64 ] o Qt | 5.24 | 2:44. 1 2:62 
7.0 0.34 | 0.33 | 0.35 | 0.37 | 0.41 | 0.48 | 0.64 | 0.88 | 1.18 | 54.12.02 |} 2.59 | 3.38 
ae 0.42 | 0.41 | 0.40 | 0.40 | 0 42 | 0.47 | 0 62 | 0.83 r.22 1.46 | 1.90 | 2.46 3.05 
7.2 O-5I | 0.47 | 0.43 | 0.41 | 0.41 | 0.46 | 0.57 | 0.77 | 1.03 | 1.37 | 1.80 | 2.34 | 2.03 j 
7.3 0.61 | 0.55 | 0.48 | 0.45 | 0.42 | 0.48 | o.s8 0.72 | 0.94 | 1.28 | 1.70 | 2.22 | 2.81 
7.4 0.73 | 0.67 | 0.60 | 0.56 | 0.54 | 0.58 | 0.67 | 0.77 | 0.92 | 1.22 | 1.62 | 2.12 | 2 69 

7.5 0.85 | 0.81 | 0.76 | 0.72 | 0.69 | 0.72 | o 78 | 0.90 | 1.01 | 1.26] 1.62 | 2.06 | 2 60 | 
7.6 1.00 | 0.98 | 0.95 | 0.89 | 0.87 | 0.90 | 0.97 | 1.06] 1.18] 1 41 ..3 2.07 | 2.53 
7.9 1.25 | 1.24 | 1.20 | 1.17 | 1.15 | 1.18 | £.23 | 2.29 | 2.41 | 1.59 | 1.82 | 2.12 | 2.50 | 
7.8 1.48 | 1.47 | 1.45 | 1.43 | 1.41 | 1.43 | 1.47 | 1.54 | 1.62 | 1.76 | 1.94 | 2 19 | 2.50 
7.9 1.70 | 1.70] 1.70 | 1.69 | 1.67 | 1.69 | 1.72 | 1.76 | 1.80 | 1.89 |] 2.05 | 2.26 | 2.61 
8.0 1.89 | 1.89 | 1.89 | 1.89 | 1. 88] 1.89 | 1 ot | 1 5:90 2.02 {.2.3§ | 2.33 
Ik a 2-05 | 2-05 | 2.05 | 2.05 | 2.04 | 2.04 | 2.05 | 2.07 | 2.10 | 2.14 | 2.25 | 2.40 | 2.50 
8.2 2:08 2:28 2.28 | 2.21 2.20 | 2.33 | 2.45 | 2.63 
Sy 3.; 2.26 2.26 2.26 2.26 2.26 2.20 2.26 2.28 2.30 2.34 2.40 2.50 | 2.67 
er a 2.38 | 2.38 | 2.38 | 2.38 | 2.38 | 2.38 | 2.38 | 2 39 | 2.30 [2.42 | 2.46 1) 2.55 1 3.70 
CASE 4: CFe=0.125, CK, Ca=CSj=0.125, CNa, Mg =0.125, CQ =0.5 
LOG p(1+X) 

Loc T 

=—2.0 | —0.5!1 0.0 0.5 Io 2.0 3.0 as 4.0 
5.6 O13 | 0.31 | 0.59 | 0.92 2.2 1 64 2.02 | 2.460 | 2.92 3.43 | 3.908 4.55 ee 
“ty ae O.1Ir | 0.20 | 0.45 | 0.78 | 1.11 1.49 | 1.90 | 2.34 | 2.80 |] 3.30 |] 3.83 4.41 5 00 
5.8 0.14 | 0.18 | 0.36 | 0 67 I.o1 r.38 | 1.78 | 2.2% 2.66 | 3.16 | 3.69 4.37 1 86 
5.9 0.28 | 0.31 | 0.38 | 0.63 | 0.94 | 1.29 | 1.68 | 2.09] 2.54 | 3.03 | 3.56 | 4.13 | 4-72 } 

6.0.. 0.32 | 0.34 | 0.40 | 0.60 | 0. go] 1.21 2.55 | 2.08 | 2.42 2.90 | 3.42 4.00 | 4.58 | 
6.1 ©.24 | 0-27 | 0.34 | | | 5.46 | 1.85 | 2.27 | 2.76 | 3.28 | 3.86 | 4.44 
6.2 0.16 | 0.19 | 0.24 | 0.40 | 0.65 | 0.96 3% 2.61 4.246.) 3.92 4.30 
6 3. O.1r | 0.14 | 0.19 | 0.31 | 0.52] 0 85 | 1.20] 1 58 | 2.01 2.48 | 3.00 | 3.58 | 4.16 
6.4 0.08 | o.11 | 0.15 | 0.24 | 0 41 | 0.72 | 1.07 2.47 | 2:88 | 2.35 | 2.87 3-44 | 4.02 
65 0.04 | 0 08 | 0.12 | 0.20 ©.34 | 0 62 | 0.96 7S 2.33 3.30 3.88 
6.6. 0.03 | 0.07 |} oO. 11 | 0.18 | 0.29 | 0.53 | 0 84 ] 1.21 1 61 2.08 | 2.59 | 3.16 | 3.74 
6.7 0.14 | 0.19 | 0.26 | 0.46 0.74 1 08 1.49 I. 2.44 3.02 3.60 
68 0.20 | 0.21 | 0.22 | 0 23 |] 0.2! ©.42 | 0.65 | 0.97 | 1.35 | 1.80] 2.31 2.87 3.46 j 
6.9 0.29 | 0.30 | 0.31 | 0.33 | 0.36 | 0.45 | 0.64 | 0.92 | 1.26 | 1.69 | 2.19 | 2.74 | 3.32 
7.0 0.35 | -0.36 | | 6.45 10-8 0.67 | 0.90] 1.22 | 1.59 | 2.08 | 2 61 3.19 
7.1 0.4% | 0.42 | 0.43 | 0.45 | 0.47 | Oo 5 0.65 | 0.384 | 3.17 | 1.51 | 1.97 | 2.40 | 3.07 
7.2 0.50 | 0 50 | 0.49 | 0.48 | 0.49 | 0 51 | 0 60] 0 80] 1.08 | 1.44 | 1.87 | 2.39 | 2 o4 
7.3 0.61 05 0.560 | 05610571059] 065 | 0 82] 1.01 1.38 | 1.81 2.31 | 2.83 
7-4 0-73 | 9.73 | 0.72 | 0.72 | 0.73 | 0.74 | 0.78 | 0.87 | 1.02 | 1.34 | 1.75 | 2.23 | 2.73 
0.99 | 0.90] 0 89 | 0 89 | 0 90] 0 92} 095] 1 01 | 1 38 | | 2.18 | 2.65 
£.35 I 31 I z.32 I.34 1 36 1.40 1.50 1 68 Ig! 4.392 
7.8 .57 | | 1.57 | 1.57 | 1.58 | 1.59 | | 1.65 | 1.72 | | 2.02 | 2.28 | 2.57 
8.0.. 2.00 | 2.00 | 2.00] 2.00 | 2.00] 2.00] 2.00] 2 o2 2.05 | 2.1% | 2.2: 2.40 | 2.61 
B.%..; 2.14 | 2.14 | 2.34 | 2.14 | 2.14 | 2.14 | 2.14 | 2.16 | 2.18 | 2.2 2.334 2:40 } 2.65 
8.2... 2.27 | 2.27 | 2.27 | 2.27 | 2.27 | 2.27 | 2.27 | 2.28 | 2.29 | 2.34 | 2.41 | 2.52 | 2.68 
8.3 2.36 | 2.36 | 2.36 | 2.36 | 2.36 | 2 36 | 2.36 | 2.36 | 2.37 2.41 2.47 2.56 | 2.72 
B48: 2.44 | 2.44 | 2.44 | 2.44 | 2.44 | 2.44 | 2.44 | 2.44 | 2.45 | 2.47 | 2.52 | 2.60 | 2.74 
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Values of (¢/g) for any temperature in the range 5 X 10° < T < 
10’ and for any density in the range 10°? < p < 104 can be obtained 
by graphical or four-point interpolation. Values for intermediate 
chemical mixtures can be obtained by making use of two properties 
of equation (16). In the first place, because of the maximum” in the 
integrand at w © 67, the most important term in ¢ is the one corre- 
sponding to the absorption limit W,, which is nearest in value to 
67. In the second place, other things being equal, the heaviest ele- 
ments have the greatest effect on the value of F,. These properties, 
plus a study of the computed results, enable one to set up an ap- 
proximate scheme for interpolation for any mixture of O, Na, Mg, 
Si, K, Ca, and Fe, diluted with H and He, subject to the following 
limitations: 


0.125 < < 0.25, 0.0625 < Csi = CKiCca < 0.133; 


0.125 < Cre < 0.3, 


where c is the proportion by weight of the element in the undiluted 
mixture. 

The results given in Table 3 for the four mixtures can be labeled 
t,, t., t;, and ¢,, respectively. The value of ¢ for an intermediate mix- 
ture is then a linear combination of the four tabulated ?’s: 


t= alT) (aiti + aot, + + ayts) 
+ B(T) (dit; + bat, + b,t,) 


Table 4 gives values of the coefficients a and 0 in terms of the con- 
centrations of the heavier elements, and of a and 6 in terms of the 
temperature. By this procedure it is possible to obtain the absorp- 
tion opacity for any specified case. 

In the high-temperature regions the absorption opacity is negli- 
gible, compared with the scattering opacity. From equation (2) 
we can find that the integrand in equation (14) has its maximum at 
w = 47',. Since the function V changes slowly, we obtain a formula 
for the scattering opacity. 


Ke © o.19(1 + X) N(4T,) cm?/gm, (18) 


44 PHILIP M. MORSE 


which is fairly accurate. The function N is given in Table 1, and 7, 
is the temperature in atomic units.’ 

In the temperature regions where x, and x, are the same order of 
magnitude, the two opacities should be added before averaging. To 
do this accurately would require graphical integration for each sepa- 


TABLE 4 


INTERPOLATION COEFFICIENTS FOR INTERMEDIATE CHEMICAL MIXTURES 
(See Eq. [17]) 


0.125 0.200 ©.250 ©. 300 
CFe 
|O.94/0.07/0. 10/0. 13/0.04/0.07/0. 10/0. 13/0.04/0.07/0. 10/0. 13/0.04]0 07/0. 10/0. 13 
a: |1.0 |o.4 fe) ° ° ° © |0.0 |0.0 |0.0 |0.0 fo) fe) fo) 
az |0.0 jo.0 | ° fo) ° © |0.6 jo.5 I I I I 
a; |0.0 10.0 jo.0 I I I I jo.4 |o.5 fe) ° ° 
a; |0.0 |o.1 I ° fo) © 10.0 {0.0 |0.0 |0.0 ° 
bh; fo) !0.9 |0.4 © |0.0 |0.0 © |0.0 |0.0 
I |0.0 I |0.0 |0.0 1 jo.8 |o.4 Jo.1 1 jo.8 |o.4 |o.1 
b, jo.1 jo.6 I |0.0 jo.5 © |0.0 jo.0 | fo 0.0 |o.0 

Usiks 650 6.2 6.4 6.4 ‘6235 6:6 Gig 

a 0760. 0.8 0.5 0.2 ©.0 073: 6.7 B20: 10:7 


rate case. However, by simplifying the form of «,(@) somewhat, it is 
possible to perform the joining integration in a manner that can be 
used in any case. We replace the step function /’, by a straight line, 
and thus make x,(w) proportional to w *. It is then possible to ob- 
tain the Rosseland average of (a/w)} + b in terms of the average 
values of (a/w)3 and b. The results in Table 2 give approximate 
values of the ratio (k/x,) as a function of the ratio (xk, x.). Once the 
values of x, and x, are determined from equations (16) and (18), the 
total average opacity x can be obtained from Table 2. 

Figure 4 gives the curves for average opacity for the Russell mix- 
ture (case 1) as functions of 7, for different densities and hydrogen 
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concentrations. The vertical scale includes the factor (1 — Y — Y) 
in order that the curves for low temperature be independent of 
and Y, the proportions of the light elements. Since the plots are on a 


) 


( 


log 


Fic. 4.— Average opacity « for the Russell mixture diluted with hydrogen, as a func- 
tion of temperature for different densities. The dotted lines are for the mixture diluted 
with 60 per cent hydrogen; the solid line, for Y = 0.30; the dashed line, for no hydro- 


gen. 


log-log scale, the slope of the curves gives directly the exponent of the 
dependence on 7. We note that at high densities the opacity is pro- 
portional to 7”? over most of the range, and that it is nearly inde- 
pendent of the density. At low densities the average dependence on 


4p 
| 
| 6 7 8 
| log T 
q 
| 
| 
| 
{ 
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T is as 77, but the exponent changes considerably with tempera- 
ture, being as large as --g/2 near a million degrees. At high tem- 
peratures the scattering predominates, and the dependence on tem- 
perature and density is less marked. 


DISCUSSION 


In the integral in equation (14) giving the mean absorption opac- 
ity, the integrand has a maximum for w somewhere between 67, 


\ 
2 

SS 

1 

0: 

8 


Fic. 5.—Values of the guillotine factor (¢ g) as function of 7 for three different values 
of p and for two different mixtures of heavy elements. The solid line is for the Russell 
mixture, case 1 of Table 3; and the dotted line is for case 2. Below the curves the differ- 
ent bound states are noted at the temperatures where the absorption due to this state is 
the most important part of the opacity. 


and 77.,." Therefore, the part of the opacity x,(w) that has the great- 
est weight in the Rosseland mean is for w s 67, & (7, 26,000). 
This behavior has been used earlier in the derivation of equation (17) 
for interpolation. It is illustrated by the behavior of the curves in 


Figure 5. 


| 
| 
) 
| 
| | 
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In Figure 5 are given curves of log (¢/g) against log (7) for three 
chemical mixtures and three densities. When log (t/z) = 0, the 
opacity follows Kramers’ law; when it is larger than zero, the opacity 
is smaller than Kramers’ value. A negative slope of the curve means 
that the opacity depends on 7 according to a negative power less 
than 7/2; a positive slope means that the opacity-curve goes down 
more steeply than a (7~*?)-curve. When the curves for different 
density overlie each other, the opacity is proportional to the density; 
when the curves for different p are spaced widely apart, the opacity 
is nearly independent of p. 

Below the curves are placed notations of bound states at points 
on the temperature scale corresponding to the absorption limit for 
the state (7 & 26,cooW.,.). According to the discussion above, the 
state noted should have a predominant effect on the opacity in the 
temperature range near the notation. That this is correct is shown in 
several places in the curves. For instance, for temperatures near 
25,000,000° (and also near 3,000,000°) the mixture having the most 
iron has the lowest value of (¢/g) and, therefore, the largest value 
of k,. In addition to the increase due to change in (¢/Z), the opacity 
is also increased because the value of T° for a mixture rich in iron is 
large. For certain temperatures an increase in the amount of iron 
by a factor of 2.5 at the expense of the other heavy elements will in- 
crease the opacity by a factor of 2. This increase does not involve 
an appreciable change in the mean molecular weight, and it would 
therefore produce an appreciable change in the predicted hydrogen 
abundances from those given by Strémgren' for the Russell mixture. 
The change in hydrogen abundance from Strémgren’s values would 
probably not be by a factor greater than 1.4, however. 

The results show that the hydrogen abundance computed for a 
star will depend on the particular mixture of heavy elements chosen. 
Therefore, the contours for constant hydrogen content on the Russell 
diagram, computed by Strémgren,' would be changed somewhat if a 
mixture of heavy elements differing from the Russell mixture were 
chosen for the stars. The change would be most marked if the as- 
sumed iron abundance were increased, and would perhaps be still 
more marked if the presence of appreciable amounts of still heavier 
atoms were assumed. 
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The general trend of the curves for (¢/g) indicates the various 
physical phenomena which are important. At the high-temperature 
end all bound states are poorly populated, and the opacity is small 
and proportional to the density, so that (¢/g) is large and independ- 
ent of density. As the temperature is reduced, the states become 
more populated, and (//2) decreases. Below T = 20,000,000" the 
K-shell absorption for iron loses its importance, and below 7 = 
4,000,000° the L-shell absorption preponderates. At about this point 
the K shells of the heavy elements begin to be filled, the effective 
nuclear charge for the L shells is thereby reduced, and (fg) increases 
again (the opacity falls again below Kramers’ value). As the temper- 
ature is further reduced, the opacity becomes less dependent upon 
density due to the filling of the L shells, so that (/ g) depends more 
strongly on the density, the curves for large density rising more 
steeply than those for small values of p. The curve for p = 0.01 even 
drops back down, owing to a bunching-together of absorption limits 
near the lower end of the frequency scale (and is therefore important 
at low temperatures). 

At temperatures below 400,000° the L shells are nearly filled, and 
the actual wave functions for the atoms depart so far from the hy- 
drogenic form that the use of an effective nuclear charge is not an 
adequate correction of equation (1) to give the correct form for the 
absorption cross-sections. Hartree wave functions, or variational 
forms, must be used for calculating the cross-sections. 

For densities less than 0.01 gm, cm‘, the quantity (/ g) will not 
differ very markedly from the curve for p(1+) = 0.01 shown, in 
the range of temperature considered. ‘The minimum value of (¢ ¢) 
will, in general, be unity; and the smaller the density, the lower must 
be the temperature before (/ g) increases again. 


SUMMARY 
The calculations outlined above make it possible to compute the 
mean opacity of a mixture of six heavy elements having’ proportions 
by weight equal to cz, diluted with proportions Y and V of hydrogen 
and helium, at a temperature 7° C and a density p gm/cm}. If the 
gas is a mixture of hydrogen and helium in any amount, with one of 
the four mixtures of heavy clements given in Table 3, we can com- 
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pute the absorption opacity x, by using equation (16) and Table 3. 
The scattering opacity x, can be computed by using equation (18) 
and Table 1. The total opacity « is then obtained by finding in Ta- 
ble 2 the value of (k/ x.) corresponding to the ratio (k./k;), which has 
been found, and multiplying (k/k.) by Ka. 

If the mixture of heavy elements does not correspond to one of 
those considered in Table 3 but has a composition intermediate 
between them, a value of x, can be obtained by the interpolation 
process given in equation (17) and Table 4. The total opacity can be 
obtained from Table 2, as before. Values of the opacity for densities 
larger than those considered in this paper can be obtained from the 
paper of Marshak." 
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A STUDY OF Het LINES IN THE BALMER 
CONTINUUM 


F. E. ROACH AND LEON BLITZER 


ABSTRACT 

Measurements of the equivalent widths of He 1 lines of the diffuse series have been 
made for 22 stars of classes Oe5 to B5. In all cases the higher members of the triplet 
diffuse series, located in the Balmer continuum, are abnormally weak. This weakening 
has been interpreted in terms of an increase in the continuous absorption coeflicient 
in the Balmer continuum. A comparison of our results with Unsild’s theoretical cal- 
culations suggests a high hydrogen content for the stars of this study. 

The two diffuse series of neutral helium offer a convenient test 
of the effect of the depression of the continuous background on the 
violet side of the head of the Balmer series (A 3647) upon the in- 
tensity of absorption lines. The first five members of the triplet dif- 
fuse series have wave lengths greater than \ 3647, while the higher 
members are in the Balmer continuum. The theoretical absolute in- 
tensities of these lines, as well as of the singlet diffuse lines, have been 
computed by Goldberg,’ who has shown that it should be possible 
to combine the two series into a single curve of growth,’ provided 
there is thermodynamic equilibrium. This is feasible, in spite of a 
rather marked difference in the damping constants of the singlet and 
triplet lines, because in neither case is the damping portion of the 
curve of growth reached in stellar spectra. 

However, when an attempt is made to combine all the singlet and 
triplet lines into a single curve of growth, it is found that the triplet 
lines in the Balmer continuum are systematically weakened, an 
effect which makes it necessary to construct two separate curves of 
growth. This was first observed by Struve and Roach} in the spec- 
trum of P Cygni. 

The spectrograms of O- and B-type stars for the present study 
were loaned to us by Dr. Struve from the collection of plates made at 
the McDonald Observatory with the two-prism quartz spectro- 
graph. Some of the plates have been examined in the near ultraviolet 


4 


t Ap. J., 90, 414, 1939. 2 [bid., 89, 623, 1930. 3 Ap. J., 90, 727, 1939. 
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by Struve, and a number were reproduced in his paper.‘ In all, we 
were able to measure spectrograms of 22 stars. 

The equivalent widths, with respect to the adjacent background, 
of all the available helium lines were determined in the usual man- 
ner. The lines in the portion of the spectrum available to glass spec- 
trographs have been previously measured by Williams* and by Rud- 
nick.® Our measures are in satisfactory agreement with both, al- 


TABLE 1 


EQUIVALENT WIDTHS OF HELIUM LINES 
(Unit =o.o1 A) 


| 
Star (HD) 3487/3470] | 4022] 4388/4144] 4009] 3926] 3872 
886 130} 120] 126) 38 | 30 | 16 123] 86 | 76 | 57 | 70 | 35 
2905 63} 58) 16 | 18] 11 7 ie 42 34 28 36 
3360 206} 226) 32 | 36 | 26.1. ... fe... 123] 89 | 80 | 68 | 68 | 27 
53904 121} 143] 104] 29 | 23 | 1§ | 12 | 8 5 
360822 10g] 115} 56!) 26] 14 9 6 | 7O 721] 35 Q | Io 
112] 78] 65] 26] 24 | 16 g 20] 14 84] so | 41 | 32 | 29 | 16 
149438 112] 108} gt} 38 | 23 84] 62 | 60 | 40 | 29 | IO 
157056 112) 113 S| 40 20 | 37 | 73 | 62 | 50! 29 9 
100702 123] 143] 116) 13 | 10 |.. 106} 82 | 83 | 62 | 52 | 25 
166182 120] 111] 34 | 29 | 27 | 19 | 11 6 
1903237° 82} 75] 64] 24 20 | 10 | 16 13 73 | 44 | 32 | 30 | 34 
198478 82} 76) 62) 21 | 23 14 | 12] 12] Of 14 85] 51 | 41 | 33 | 40 | 26 
205021 158] 119g] 103) 27 | 17 | 17 6 5 Oy Gee Se 103} 74 | 68 | 51 | 45 | 23 
200165 -| o7| 82 70 19 | 22] 14] 11 8 12 87] 54 | 40 | 37 | 44 | 32 
214080 82] go} 18 | 13 5 Gai 44-1 94 1.27 | 


* Data for 193237 are from Struve and Roach, A p. J., 90, 727, 1930. 


though we find that, for the stronger lines, Rudnick’s values are. 
slightly lower than ours. In Table 1 we record the values of the 
equivalent widths, 4, adopted in this study. We have averaged our 
measures with those of Rudnick and Williams whenever there were 
common values, and have used our measures alone for all other 
lines. The line \ 3705 is blended with 7/7 \ 3704 and was therefore 
omitted from Table 1. 

For cach star the observed values of log AX were plotted against 
the corresponding log Y{ as defined by Goldberg.? The values of 
log X{ in our Table 2 are from his Table 1, except for the higher 


4 Ap. J., 90, 699, 1930. 
SAp. J., 83, 279, 1936. 6 Ap. J., 83, 439, 1936. 
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members of the triplet series, for which they were determined from 
Goldberg’s asymptotic expression for the oscillator strengths. Fig- 
ures 1 and 2 show the resultant plots for some of the stars investi- 


TABLE 2 


LOG X$ FOR HELIUM LINES 
(From Goldberg) 


Transition x log Xf Transition log 
23Pp— gD..... 4472 +1.38 || 2P—4'D..... 4922 +0.7 
4026 +1.05 4388 +0.45 
3820 +o.81 oD 4144 +0.15 
3705 +0.59 7D 4009 —0.03 
8D 3634 | +0.43 83'D 3920 | —o.22 } 
3587 +0.27 3872 —0o.40 
10D) 3554 | +0.13 
123) 3512 —0.09 
1331) 34098 —o.18 
143D 3487 —o. 28 
153D 3479 —0o.39 
-30 T T T T -30 T T T T T T T ( 
e LEO v ORI | 
° y 
be 4 é 
00 10 20 00 10 20 
—30 T T T T T T T -30 T T T T T 
Cc MA p CEP 
e 4 
00 10 20 00 1.0 20 


Fic. 1.—Curves of growth. The ordinates are log /; the abscissae, log V/. Open 
circles designate triplet lines in the Balmer continuum; crosses, other triplet lines; closed 
circles, singlet lines. 


gated. It is apparent that the points representing triplet lines on 
the violet side of \ 3647 fail to lie on the main curve of growth. 
Goldberg's theoretical curve of growth for 20,000° K was next 
fitted to the observed points by moving it vertically and horizontal- 
ly, these shifts corresponding to the V and A, respectively, of his 
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paper. In the cases for which the data were sufficient, this procedure 
was followed independently for both sets of points. It was found in 
this way that for each star the same V will serve for both but that 
different A’s are required. In case there were too few points in the 
continuum to give a completely independent determination of both 
V and A, it was assumed that V has the same value for the points in 
the continuum as in the unaffected region, and A was then determined 
by a horizontal shift of the theoretical curve. The difference in A be- 


—3 Tr T T T T T -30 + T + 
9 CEP | 55 cyG 
= 1 it it 
-30 T T T T T T -30 T T + 
LYR HER 
4 — 
° 
0° 10 20 00 ike) 20 


Fic. 2. Curves of growth. The ordinates are log 4/\; the abscissae, log X$. Open 
circles designate triplet lines in the Balmer continuum; crosses, other triplet lines; 
closed circles, singlet lines. 


tween the two curves, which we shall call Z, may be interpreted as be- 
ing due to a smaller total number of helium atoms available for ab- 
sorbing the ultraviolet lines. This weakening of the ultraviolet lines 
with respect to those on the red side of \ 3647 is the principal obser- 
vational result of this study and is shown in the fourth column of 
Table 3. The sixth column lists the luminosity classifications which, 
with the exception of those in parentheses, are taken from Williams.’ 
In the seventh column are the measures of the so-called Balmer 
discontinuity, ), in the continuous spectrum,* where 


Int (3647 + ©) 
Int (3647 — 


D = logie 


7 Op. cit., p. 305. 
* Cf. the papers in the Annales d’ustrophysique by Arnulf, Barbier, Chalonge, et al. 
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That the cause of Z is the depression of the continuous background 
in the Balmer continuum becomes obvious upon examination of 
Figure 3, which correlates Z with the D of Chalonge and his collabo- 
rators.* The increase of Z with D indicates that this is the primary 
correlation. 

TABLE 3 
SUMMARY OF RESULTS 


STAR 
Vv Z SPEc. Lu. D* 
HD Name 
y Peg 0.38 0.08 B2 i +o0.13 
2905 eae x Cas 02 0.71 Bo (g) + .03 
33600 ¢ Cas 0.08 B3 i + .12 
5304 Cas °.68 Bop (g) — .19 
30512... v Ori 0.45 Bo (d) 
30822 Ori .4I Bo (d) 
30861. . Ori (br) 0.50 Oes i + 03 
44743... BC Ma 0.56 Br i + .o8 
52080. . eC Ma 28 Br (i) + 
o? C Ma 0.71 Bsp (g) + .o8 
p Leo 0.44 Bop g + .o4 
149438 t Sco .43 0.59 Bo d + (05 
157056. Oph 34 0.05 B3 (d) 
160762. Her 1.30 B3 i + .17 
164353. . 67 Oph 0.04 g + .13 
106182. 135 0.9! B3 (i) 
180163... n Lyr 1.26 B3 (i) 
193237... P Cyg .05 0.54 Bip (g) — 
198478. . 55 Cyg °.60 Be g + 
205021... B Cep .44 0.88 Bi i + 
g Cep .20 Bap g + .o6 
10 Lac 0.41 0.46 Oe5 i +o o4 


* From the papers by Arnulf, Barbier, Chalonge, ef al. in Annales d'astrophysique. 


In a study of lines of 77 1 in the Balmer continuum Struve and 
Sherman’? have observed a similar weakening. They have shown that 
the weakening, expressed as numbers of atoms, is to be interpreted 
as a change in the continuous absorption coefficient. This will be 
true in our case, provided there is no stratification of helium with re- 
spect to hydrogen. We have calculated the relative populations of the 
helium level, 2'P, as a function of the optical depth, following Swings 


9 Ap. J., 91, 428, 1940. We are grateful to the authors for the opportunity of exam- 
ining this paper in manuscript. 


| 
| 


Het LINES IN THE BALMER CONTINUUM 55 


and Chandrasckhar,'’’ on the assumption of an effective temperature 
of 20,000° K. Upon comparing these results with those of Swings 


HER 
% CAS YPEG 
9 CEP 
67 OPH 
x CE 
CAS OMA 
x 
55 CYG T Sco 
A ORI € 
¢ LEO 
| 
0.00 0.05 0.10 0.15 D 


Fic. 3. The observed weakening in the Balmer continuum plotted against D. The 
points for P Cyg and y Cas are not included. The symbol code is the same as for 
Vig. 4. 


TABLE 4 


Spectral Class Mean Z Spectral Class Mean Z 
0.48 B3. 1.03 
Bo. 59 Bs. 0.82 


* The value of Z for Ao is from Struve and Sherman (Ap. J., gi: 428, 
1910). Their measures for » Leo, a Cyg, and a Lyr are included in the 
mean. 


and Chandrasekhar for hydrogen, we obtain the entries of Table 5. 
It is immediately apparent that the relative stratification is negli- 
gible for normal and giant stars and is very small for dwarfs, so that 


10 97, 24, 1936. 
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we are justified in assuming that the ratio, 7, of the mass absorption 
coefficient within the line, «,-+ ¢,, to the mass absorption coefficient 
for continuous radiation, x, does not vary with optical depth. 

If we employ the nomenclature of Struve and Sherman,’ we have 


Ky + Oy NV (Kat + Oat) 


K K 


where V is the number of active atoms per gram of gas. Hence, while 
we have observed an effect which has been interpreted as A log .V, 


TABLE 5 


RATIO OF ABUNDANCE: HYDROGEN (BALMER) 
TO HELIUM (2'P) FOR Te= 20,000° K 


Tt (Optical Depth) net Giant Dwarf 
mediate 

| 5.400 5.380 7 325 
| 5-398 5.382 6.849 
| 5.360 6.159 
0.666. . 5-252 5.683 
| 5-373 5.6061 
5-302 5-541 
2.0... 5.462 5.553 


or Z, the cause of the measured apparent decrease in NV is an actual 
increase in the mass absorption coefficient, x, due to the photoioniza- 
tion of hydrogen. 

The observed Z of this paper may, therefore, be compared with the 
theoretical value of the increase in the mass absorption coefficient 
at the Balmer head. Unsold'' has computed this increase for pure 
hydrogen (his Table 30, p. 127) and for a ratio of hydrogen atoms to 
metallic atoms of 13.7:1 (his Table 28, p. 121). In Figure 4 are 
shown our observations, together with the two curves from Unsold, 
as a function of 6( = 5040/7). With the exception of two high points 


1 Physik der Sternatmosphiren, Berlin: Julius Springer, 1938. 
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(yn Lyrae and « Herculis) the observed data scatter about the curve 
for pure hydrogen. In principle, it should be possible to determine 
in this manner the percentage of hydrogen in stellar atmospheres, 
especially if the weakening of lines in the continuum is measured for 


Z T | T | T i T 
=o 

t M = 
7 LYR 
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° 
‘ 

@ 
I.C YPEGO OfCAS acyc@atyr = 
102 HER’ 67 OPH 
= 
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a x x H 4 
= 137 
@ OR! Oe oPH M 3 
T SCO X55 CYG = 
ac MA 
4 e cmaQ® a 
LAC v ori 

§ LEO = 

v SAG 
4 

B0BIB283B85 AO 
al ~ 
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Fic. 4. ~The observed weakening in the Balmer continuum as a function of @ 
(=soyo 7). Open circles designate dwarfs and intermediates; crosses, giants; closed 
circles, the means by spectral class. The data for a Cyg, a Lyr, » Leo, and v Sag are 
from Struve and Sherman.’ The two solid lines are from Unsild’s calculations of the 
increase of the continuous absorption coefficient in the Balmer continuum; the upper 
curve is for pure hydrogen; the lower, for 7/M = 13.7. 


stars of spectral classes later than Ao. From our data, together with 
the points from the paper by Struve and Sherman, omitting v Sagit- 
tarii, we have also plotted the averages by spectral class (solid 
circles). The results favor a much higher hydrogen content than 
H/M = 13.7. 

Those points which fall considerably above the mean curve, as 
those for » Lyrae and ¢ Herculis, are of particular interest. That 
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these points are associated with a Stark effect is ruled out, since 
there is, in general, no correlation of the weakening with absolute 
magnitude, a result which is impossible if Stark effect is involved.” 

On the other hand, the vertical displacements of our curves of 
growth, V, are definitely functions of luminosity. The mean value 
of V for the dwarfs and the intermediate stars is 0.35, and for the 
giants, 0.20. Goldberg's? Figure 5 indicates the same result, and he 
suggests that ‘“‘in the low-luminosity groups the value of V is a meas- 
ure of the average Stark broadening in the diffuse series, while in the 
giants V measures the degree of turbulence in the stellar atmos- 
phere.” 

Attention is called to the stars P Cygni and y Cassiopeiae. For 
these stars the Balmer continuum is found to be emission rather than 
absorption: D is negative. Nevertheless, we find appreciable values 
of Z, about the average for the respective spectral classes. This is 
not unexpected: the continuous absorption of hydrogen within the 
shells of these stars should produce a weakening of //e-lines if the 
latter are also produced in the shell."$ 


We are indebted to Miss Frances Sherman, who made a number 
of the microphotometer tracings for us at the Yerkes Observatory. 
We wish to express our appreciation also to the McDonald Observa- 
tory and to Dr. Struve for lending us the spectrograms for this study, 
and to acknowledge the value of discussions of the problem with 
Dr. Struve. 


STEWARD OBSERVATORY 
AND 
DEPARTMENT OF PuysIcs 
UNIVERSITY OF ARIZONA 
March 1940 


2 The role of Stark effect in the calculations of his theoretical curves of growth 
of B-type stars was discussed by Goldberg.? He assumed an average Stark broad- 
ening for all the diffuse series lines; whereas a rigorous solution should consider 
the increase of Stark broadening with serial number. It is apparent that if Stark broad- 
ening is significant it should enhance the equivalent widths of the higher members of 
the series and thus produce an effect the opposite of that which we have observed. 


13 Struve, Proc. Nat. Acad. Sci., 26, 117, 1940. 


HETEROCHROMATIC SURVEYS OF INFRARED 
VARIABLE STARS 


CHARLES HETZLER 


ABSTRACT 


Seven late M-type variables—T Cep, R Cas, R Leo, R Aql, o Cet, X Oph, and R Aqu— 
are included in this first of a series of papers on the heterochromatic photometry 
(particularly at 8500 A) of about forty infrared variable stars. The average relative 
magnitude ranges are 4.94, 4.29 (first two stars only), 4.73, 2.43, and 1.95 in the ultra- 
violet, photographic, visual, 8500 A, and g1oo-g250 A regions, respectively. 

The tri-color condition for a black star is shown to be that the plot (designated as 
“‘slope-index”’) of any two color indices gives a straight line whose slope depends only 
on the wave lengths and which passes through the origin. The indices involving the 
ultraviolet and infrared magnitudes come closest to meeting this condition, and temper- 
atures from these indices are therefore considered most reliable. Using these tempera- 
tures (averaging 2800° at maximum and 2100° at minimum), the differential absorption 
in the visual region between maxima and minima is calculated for each of the variables, 
being least for R Cas and R Leo, and around 2 mag. for each of the others, the average 
being 1.4 mag. 

It is shown that, for a black variable, star magnitudes at any wave length plotted 
against those at any other give a straight line of slope equal to the ratio of the wave 
lengths. Plots of visual against 8500 A magnitudes give straight lines (except for curva- 
ture near the maxima of stars having a delayed infrared drop and near the minima of 
stars having companions) of medial slopes 1.92, 1.87, 1.82, 1.94, 1.95, 1.80, and 1.92 for 
the seven stars, respectively, averaging 1.89, as compared to the theoretical value 1.54. 
The delayed infrared drop and the steeper slopes near maxima, particularly in o Ceti, 
are interpreted as an effect of band absorption which is not synchronous with, or pro- 
portional to, the temperature change; and the generally steeper slopes are due primarily 
to a superposition upon the black-body slope of an absorption which, on the average, 
increases linearly toward minimum; but changes of diameter would also contribute to 
this effect. The ultraviolet magnitudes plotted against the infrared magnitudes give 
straight lines of slopes near the theoretical value of 2.3; hence, band absorption is prob- 
ably the chief factor of departures in the visual region. The y-intercepts give the 
radii and possible changes in radii, if we know the individual absolute magnitudes. 

Extrapolation of the medial straight-line portions of the visual minus infrared and 
ultraviolet minus infrared slope-index curves yield the visual and ultraviolet magni- 
tudes, and hence the colors, of the companions to o Cet, X Oph, and R Aqu. 


In previous articles' the instruments and methods of heterochro- 
matic photometric surveys, particularly in the infrared, have been 
briefly described. The filter and emulsion combinations and the re- 
sulting effective wave lengths have already been indicated. This is 
the first of a series of articles which will give the magnitude tables 
and corresponding curves of about forty infrared (long-period) varia- 


"Ap. J., 82, 75, 1935; 83, 372, 1936; 86, 500, 1937. 
59 


60 CHARLES HETZLER 


bles, particularly at 8500 A, but also to some extent at gioo A and 
9250 A and at 3725 A, with visual-curves of the A.A.V.S.O. (kindly 
furnished by Professor Leon Campbell, of Harvard Observatory ) 
and, for a few stars, with photographic-curves by Miiller.? In this 
paper we present the data for the variables T Cep, R Cas, R Leo, 
R Aq], o Cet, X Oph, and R Aqu. Since all of these stars have spec- 
tra of late M type, it is of interest to compare them; and since the 
last three have known or suspected companions, certain deductions 
can be drawn from any departures from regularities found in the 
other four. 

Although it is generally recognized that stars do not radiate at 
their surfaces strictly as black bodies because of absorption bands, 
particularly of 77O in M stars and of CN in N stars, the amount and 
individual behavior of these deviations are little known, and compar- 
ison with some standard is fundamental in understanding them. The 
equation for the absolute magnitude of a black star at any wave 
length is (to an approximation accurate enough for our purpose) :* 


—— — 5logR. (1) 


The absolute visual magnitude of the sun, for which R is taken as 1, 
is known; hence we may readily compute C, for the effective wave 
length of the eye, about 5300 A for the sun. Having determined the 
visual constant, C,, the constants for other effective wave lengths 
can be computed by considering that for Ao stars, for which 7° is 
fairly well established at around 11,000°, all absolute magnitudes are 
equal. The values of C for the various wave lengths we are consider- 
ing (adopting 5500 A for the visual for red stars) are: C,,; = —1.23, 
Cy. = —0.72, C;; = —0.08, C;, = 0.02, Css = 0.93, Cor = 1.04, and 
C,; = 1.07, the subscripts defining the spectral region in hundreds of 
angstroms. The difference in magnitude at any two wave lengths is 
an index of the color and temperature. For various combinations 
of magnitudes at different effective wave lengths the tempera- 


2 Pub. Ap. Obs. Potsdam, 29, 2, 1938. 


3 Russel, Dugan, and Stewart, Astronomy, p. 733, 1927. 
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ture indices are readily transformed into the following tempera- 
tures: 


T= 72007 _ __18,350 10,000 

T 2-53 + 0.64 T + 1.65 T 55-85 + 0.91 (2) 
pa 3,800° 5,300° 

I + 2.16 | 57-93 + 2.30. 


These equations will be used in our considerations of the empirical in- 
dices from the data of the curves. 

Starting again with the general equation (1), it is readily seen that 
for the two wave lengths, A, and X.,, 


2 
My, = + Jo. — slog slog 8) (3) 


and 


(4) 


5(A, 


These equations are independent of the temperature. Thus, a 
plot of the absolute magnitudes of a “perfect variable” at two wave 
lengths is a straight line of slope equal to the ratio of the wave 
lengths. The y-intercept of this straight line depends upon the radius 
as well as upon the wave lengths. Such plots for actual stars will be 
designated as “‘bi-color,” or 2X, slo pe-indices. Changes of radius mean 
different intercepts at different points, amounting to an overall change 
of slope; but the slope-index curves would still be straight if the 
change in radius is uniform. The radii and changes of radii can be 
determined from the intercepts if the absolute magnitudes are 
known, as they are for a few long-period variables, such as o Ceti. 
This star is the only one whose radial velocity has been measured, 
showing a range from 64 to 52 km, sec, the greatest velocity of ap- 
proach being at minimum. Differential band absorption in actual 
variables would also alter the ideal slopes, since both absorption and 
radii increase toward minimum; and only if the amount of the one is 
known, as in o Ceti (if the change in radius is real), can the effect of 
the other be inferred. 
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It is evident that some criterion independent of the radius is need- 
ed to decide between changes of radius and of absorption, since dis- 
agreement of the slope with the ratio of the wave lengths does not 
necessarily mean deviation from a black body if the radius changes 
with the temperature, nor does agreement indicate radiation like a 
black body if absorption and changes in radius compensate one an- 
other. Considering again equation (1), this time at three wave 
lengths, we have 


C I 
My, Ih, = + 
and 
C I 


The left-hand members of these equations are the color indices, /,, 
and /,,. For an Ao star the color index is zero between any two 
wave lengths. Thus, at the temperature 7',,, 


C fx I C ys I 


Subtracting these equations from those preceding, we have 


I I I I e I I I I 
whence 


The tri-color condition for black-body radiation is, therefore, that a 
plot of any two color indices as ordinate and abscissa yield a straight 
line passing through the origin with a slope involving only the wave 
lengths, independently of RK and 7. We shall apply this ¢éri-color 
slope-index to our observational color indices. 
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Figures 1 3 show the heterochromatic-curves for the seven varia- 
bles. The corresponding data are given, together with the compari- 
son stars and their adopted magnitudes, in Tables 1-20. As indi- 
cated in the notes to the tables, the earlier portions of the curves at 
8500 A (P-emulsions) and the earlier points at 9250 A (Q-emulsions) 
are from data obtained with the 30-inch refractor of the Allegheny 
Observatory. All other observations were made with the 24-inch 
reflector of the Yerkes Observatory. With the Q and M emulsions 
(effectively at 9250 A and gt1oo A), because of their wide range of 
sensitivity, the bluer stars are of significantly shorter effective wave 
length, particularly with a reflector where there are no focal differ- 
ences. This makes the bluer comparison stars systematically too 
bright, for these wave-length regions,‘ relative to the variables. 
Most of the comparison star magnitudes at 8500 A were obtained by 
direct or indirect comparison with Vega, but at the other wave 
lengths they have been largely computed from visual magnitudes of 
the Henry Draper Catalogue. Supposing the bluer comparison-star 
magnitudes to be correctly computed for gtoo A and 9250 A, the 
variable-star magnitudes will be relatively too faint; but since the 
values adopted for the comparison stars are somewhat uncertain, no 
corrections for differences in effective wave length were attempted. 
The ranges of the variables are independent of the comparison stars 
whenever the same ones are always used. The streakiness and poor 
keeping qualities and possible differences in different batches of the 
M and Q emulsions are also to be considered as contributing to any, 
otherwise unexplainable, irregularities in the results. 

The magnitudes at g1oo A, at 8500 A, in the visual, and in the 
ultraviolet are shown as solid circles, while those at 9250 A and the 
photographic magnitudes of Miiller are shown as open circles. The 
weights of the points are indicated in the tables. The visual points 
are ten-day means which tend to smooth those observations. There 
are a number of instances where observations at 8500 A which de- 
part from the curves are associated with similar departures in the 
ultraviolet or at groo A and g250 A, and possibly many of the irregu- 
lar fluctuations are real. These fluctuations, however, are relatively 
insignificant. 

4 Ap. J., 86, 511, 1937. 
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Fic. 1.—Heterochromatic magnitude curves of T Cep (above) and R Cas. JD 
2427200-9200. Top to bottom: g200 A infrared (M plates solid and Q plates open 
circles), 8500 A infrared, visual (A.A.V.S.O. ten-day means), 3700 A uv (solid) and 
ptg. (open) by Miiller (see Tables 1-6). 
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lic. 2. Heterochromatic magnitude curves of R Leo (above) and R Agql. JD 
2427200-9200. Top to bottom: g200 A infrared (M plates solid and Q plates open 
circles), 8500 A infrared, visual (A.A.V.S.O. ten-day means), and 3700A uv (see 
Tables 10-15). 
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Fic. 3.—Heterochromatic magnitude curves of o Cet (above), R Aqu (left, below) 
and X Oph. JD 2427300-9200. Top to bottom: 9200 A infrared (M plates solid and () 
plates open circles), 8500 A infrared, visual (A.A.V.S.O. ten-day means), and 3700 A 
uv (see Tables 7-9 and 15-20). 
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TABLE 1 TABLE 2 
T Cephei at 8500 A T Cephei at 3700 A 
JD | Mag. |Quality] JD | Mag. |Quality|| JD | Mag. | Quality) JD Mag. | Quality 
242+ 242+ 242+ 242+ 
7475 1.58 6282 | 1.29 ft 8101] 10.91 6s502/ 11.99 fr 
490 | 1.41 vg 293 | 1.06 £ 114] 10.91 g 503] 11.97 fg 
498 | 1.29 g 294 | 0.97 g 131] 10.48 f 517| 11.97 fg 
0.65 fg 303 | 0.87 g 311 6.76 fg 10.61 
514] 1.06 g 310 | 0.96 ve 323] 6.90 f 715] 6.70 fg 
516 1.09 & 324 | 0.72 fg 335] 7.03 vg 730 7.29 f 
527 1.06 p 335 | 0.91 ve 351 7.78 g 738 7.26 vg 
530 | 0.99 336 | 0.91 g 391] 8.97 878 | 11.44 g 
537 | 0.87 g 344 | 0.93 pf 
549 | 0.93 g 351 | 0.90 4 The comparison stars with their computed 
558 | 1.29 g 362 | 0.93 g ultraviolet magnitudes were: HD 201731 A2 
565 1.48 f 363 1.29 3 6.76; 201684 A3 7.95; 202614 AO 8.28. 
570 1.38 377 
571 1.29 € 391 2. 4 
575 1.47 f 404 g 
581 | 1.58 f 443 | 3.15 g 
sea | 1.40 pf aes | 3.45 T Cephei at 9250 A and at 9100 A 
589 | 1.58 g so2 | 3.41 f 
602] 2.01 € 514] 3.39 vg JD Mag. | Quality} JD Mag. | Quality 
620 | 2.44 fg 552 | 2.63 f 242+ & Comp.| 242+ & Comp. 
639 | 2.73 g 614] 2.01 g 
654 3.00 & 619 | 2.01 vg 7537| 0.16 p 1 | 8502 3.42 fg 2 
679 3.16 vg 633 | 2.23 5 558; 0.45 es 523 2.64 gee 
680 3.15 & 637 1.99 vg 559] 0.45 ¢ 3 552 2.33 vg 3 
692 3.45 ft 647 2.02 4 717 1.27 : ae 709 | -0.09 fg 1 
698 | 3.32 4 661 | 2.03 f 718; 1.51 fg l 715] 0.46 fg l 
701 3.14 669 1.86 8304 0.33 730 0.57 
705 3.23 fg 682 1.44 vg 310} 0.32 fg 2 737 | O.31 pl 
709 3.31 fg 695 1.13 ve 0.47 fg 2 746 | 0.62 gel 
719 3.15 fg 697 1.13 Pp 324| 0.62 fg 2 761 0.76 el 
727 3.15 g 707 | 0.76 4 336] 0.66 fg 2 868 2.23 fg l 
39 | 3.11 pf 714 | 0.87 vg 341] 0.53 | fg 2 | 878] 2.62 
757 2.96 730 1.09 vg 0.72 | 
777 | 2.3 f 738 | 0.88 g 
799 | 1.80 € 746 | 1.29 g The Allegheny refractor was used through JD 
83) 1.71 6 754 | 1.64 g 2427718; thereafter the Yerkes reflector. Q 
861 1.2% p 761 1.30 € plates (9250A) were used through JD, 2428552; 
880 | 1.53 zg 813 | 2.73 vE thereafter M plates (9100A). The comparison 
888 1.11 VE 842 2.96 vg stars were: 1 = a Cep (2.30 mag. for the 
915 | 0.96 vE 858 | 3.32 ve refractor and 1.90 mag. for the reflector); 
915 0.88 fg 868 3.48 fg 2 =a Lyr (0.00 mag.) and 3 = a Aur (-0.74 
gel | 0.87 g 878 | 3.45 vE mag.). 
93 0.86 4 893 3.45 VE 
947 | 0.87 VE 902 | 3.47 € 
948 | 0.99 £ 918 | 3.39 zg TABLE 4 
964 1.25 & 926 3.24 4 
979 | 1.61 fg 979 | 2.38 fg R Cassiopeiae at 8500 A 
965 1.63 ve 991 2.31 fg 
goor | 2.11 ve 999 | 2.24 fg JD Mag. |Quality| JD Mag. |Quality 
O16 | 2.31 g 9011 | 2.2: fg 24e+ & Comp. | 242+ 
034 2.73 fg 02) 2.2 fg 
053 2.96 ve 041 2.15 fg 7227 2.57 pet T375 2.72 ve 
068 3.12 ve 2.01 fg ra 1.69 404 3.15 
090 3.15 vE 059 2.00 vE 245 1.57 ¢ 2 430 3.55 € 
1l4 | 2.69 vg 070 | 1.71 fg 255] 1.57 f 2] 453] 3.56 € 
138 ae. fe 077 1.48 fg 261 1.15 pf 3 457 3.55 vg 
134 | 2.24 109 | 1.19 272] 1.29 g 3 464 | 3.57 
197 | 2.02 p 112 | 0.75 fg 282 1.43 fg 3 477 3.95 pf 
199 2.01 Py 120 0.85 fg 290 1.44 fg 3 495 4.45 vp 
259 1.44 ? 125 | 0.87 fg 291 1.71 fe 581 3.03 ft 
e271 1.49 4 132 | 0.95 & 300 2.01 pf 614 2.72 f 
306 | 1.72 fg 637 | 2.31 4 
S22 1.73 ve 643 1.90 4 
The Allegheny 30-inch refractor was used 341 2.13 ve 651 1.59 ve 
through JD 2428068; thereafter the Yerkes 351 | 2.43 fg 3 667 | 0.86 f 
24-inch reflector. The comparison star was 352 | 2.30 g 669 | 0.66 é 
a Cep (2.01 mag.). 359 | 2.31 g 672 | 0.64 g 
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TABLE 4 (continued) TABLE 5 (continued) 
R Cassiopeiae at 8500 A R Cassiopeiae at 3700 A 
JD Mag. Quality | JD Mag. Quality JD | Mag. Quality | JD Mag. Quality 
242+ 242+ 242+ 242+ 
7677 0.77 g 8443 3.44 vg 8552 8.48 fg 6819 15.30 f 
680 0.65 vg 466 3.25 g 566 8.70 ¢ 820 15.25 ft 
687 0.57 €& 501 2.39 fg 761 13.65 £ 957 7.58 f 
690 | 0.76 g 514 | 1.89 g 
6398 0.87 fg 534 0.87 g [The comparison stars with their assumed or 
701 | 0.86 & 544 | 0.87 pf derived ultraviolet magnitudes were: BD +50° 
709 | 0.96 2 552 1.29 4 4208 = HD 224648 BO 7.17; BD +50° 4198 = HD 
727 | 1.16 f 556 | 1.30 g 224456 BY 8.82; BD +50° 4200 10.98 (faint 
741 | 1.34 fg 566 | 1.16 vg companion to R Cas); BD +50° 4202 9.42. 
758 | 1.59 g 580 | 1.29 vg 
777 | 1.89 g 594 | 1.59 g 
790 | 2.22 pf 602 | 1.60 g TABLE 6 
799 | 2.21 fg 609 | 1.58 & 
830 | 2.61 f 661 2.16 g R Cassiopeiae at 9250 A and at 9100 A 
832 | 2.54 vg 682 | 2.54 g 
919 | 3.72 vg 697 | 2.73 vE JD Mag {Quality | JD | Mag. | Quality 
931 | 4.13 g 715 | 2.92 g 242+ & Comp. | 242+ & Comp. 
934 3.96 ve 740 3.49 ve 
964 | 3.78 vg 72 | 3.56 g 7380 1.93 | fg 1 | 8545 1.14 fg 4 
966 3.95 4 761 3.56 vg 8317 3.24 os 547 0.72 g 32 
985 | 3.39 pf 813 | 4.09 g 324 3.53 | pf 3 552 | 1.11 ge 4 
8002 | 3.55 fg 814 | 3.87 vg 332 3.26 | fg 3 556 | 0.97 fe 4 
016 3.28 € 820 3.68 pf 352 3.29 f 3 559 1.64 fg 4 
034 3.05 g 840 3.03 € 353 2.68 f 3 566 1.54 ge 4 
044 3.04 i§ 847 3.44 vg 378 3.94 fg 3 740 2.83 es 4 
062 | 2.50 vg 858 | 2.83 vg 361 3.44 | pf 4 752 | 3.59 ve 1 
068 2.12 fg 870 2.85 vg 384 4.26 fg 4 761 3.07 ma 3% 
0388 1.13 ve 879 2.44 & 391 4.71 f 4 813 3.22 | a | 
090 1.10 f 893 2.54 fg 405 3.58 f 4 814 3.22 fg 1 
094 1.02 ve 902 2.54 ft 443 3.00 fg 3 954 0.70 a | 
113 0.99 ve 911 2.31 fg 534 1.14 fu 957 0.52 fe 2 
128 | 0.97 ve 925 | 1.75 fg 
140 | 1.06 p 954 | 1.06 fg [he Allegheny refractor was used for JD 
156 1.16 vg 957 1.06 fg R427380. Q emulsion was used through JD 
186 | 1.60 Pp 964 | 0.87 r 428566; thereafter M emulsion. The compari- 
187 1.59 g 969 0.88 fg son stars were: 1 = B Cas (1.93 mag.); 2 = 
197 | 1.89 g 975 | 0.76 fg bh Aql (0.49 mag.); 3 = a Lyr (0.00 mag.) and 
213 | 2.12 pf 981 | 0.86 g | =a Aur (-0.74 mag.). 
303 3.99 g 994 0.87 pf 
316 | 3.90 g 998 | 1.21 pr 
319 | 4.15 vp 999 | 0.87 fg TABLE 7 
324 | 4.14 vg 9019 | 0.865 f 
332 | 4.14 vg o29 | 1.01 fg o Ceti at 8500 A 
345 | 4.39 ve 041 | 1.40 fg 
362 | 4.15 vp 055 | 1.65 fg JD | Mag. |Quality | Jb | Mag. |Quality 
363 | 4.56 4 059 | 1.64 f 242+ e4e+ 
379 | 4.64 pf 069 | 1.73 fg 
381 | 4.44 g 078 1.86 fg 7315 1.94 4 7485 | -0.90 g 
391 | 4.38 g 116 | 2.22 fg See 1.94 vE 489 | -0.95 & 
404 | 4.38 g 133 | 2.45 fg SSS 2.35 pf 496 | -0.69 é 
335 2.38 pr 498 -0.77 g. 
The Allegheny refractor was used through 345 2.19 ‘g 503 | -0.75 f 
JD 2428068; thereafter the Yerkes reflector. 346 1.95 f 511 | -0.92 & 
The comparison star was B Cas (1.59 mag.) 349 1.95 ve 657 1.46 fg 
except: 1 =a Lyr (0.00 mag.); 2 =a Aql 357 2.17 € 677 1.61 pf 
(0.43 mag.); and 3 =a And %.61 mag.). 358 1.96 pf 679 1.75 fg 
367 2.37 vp 680 1.65 ve 
374 1.94 £ 687 1.75 f 
TABLE 377} 1.81 g | 698| 1.63 
404 1.37 g 701 1.67 g 
R Cassiopeiae at 3700 A 405 1.65 709 1.60 
JD Mag. | Quality! JD Mag. | Quality|} 430 0.52 f 13 1.46 £ 
242+ 242+ 454 | -0.43 4 719 1.37 fg 
458 -0.49 4 731 1.15 f 
8092 8.11 f 18345 | 14.65 rg 464 | -0.80 g 740 0.68 f 
105 8.14 fg 361 | 15.68 f 467 | -1.03 f 757 0.23 ¥ 
131 7.62 f | 534 | 8.69 fg 477 | -0.92 vp 769 | -0.20 g 
325 | 14.00 | 550 8.53 fg 483 | -0.92 782 | -0.86 
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TABLE 7 (continued) 
o Ceti at 8500 A 


Notes to Table 8: 

The comparison star for all observations ex- 
cept those of 1933 was HD 18884 for which 
the magnitude 5.65 was adopted, based upon 
an adopted ultraviolet magnitude of 4.93 for 


HD 18883. The 1933 magnitude was based upon 
Mag. Quality || observations at Allegheny (Ap. J. 82, 84, 
1935) plus an adopted ultraviolet magnitude 
7789 | -0.96 pf e460 | -1.08 fg of 9.90 for BD -3°357 A3. 
799 | -0.77 fg 466 | -1.18 fg 
807 | -0.92 g 468 | -1.03 fg 
820 | -0.75 fg 471 | -1.15 ve TABLE 9 
-0. 
o Ceti at 9250 A and at 9100 A 
846 | -0.76 491 | -1.15 
| | sol| -1.11 Jp [| Mag. [Quality] JD | Mag. [quality 
865 | -0.58 ve 506 | -1.15 ve |) 242* 242+ 
ae || 7372] 0.32 g | 8477 | -1.43 ve 
034 | 2.08 552 | -0.39 pf 492 | -1-33 
0329 1.94 vg 557 -0.40 g 430 -1.12 vp 501 -1.23 
040 2.04 g 590 0.52 g 458 | -1.54 ve 506 | -1.35 fg 
043 1.89 vg 595 0.50 g 467 -1.55 fg 524 -1.25 
os2 | 1:90 aiz | -0.20 a 709 | 0.23 g 814 | -0.76 fe 
o70| 1.74| ve | 814] -0.20] vg || 850) -1-44 
o72| 1.76] pf | 820] -0.18 g 646 -2.65) pf fe 
094 | 0.94 840 | -0.43| pr -1-12 | Si 
113 | -0.46| fg | 842 | -0.63 | ve i 
128 | -0.83 846 | -0.43 g | 
140 | -0.75 849 | -0.57 g 471 -1.45 975 | 0.45 ve 
The Allegheny refractor was used through JD 
is6 | -0.63 p g6e7 | -0.45 g 2427646; thereafter the Yerkes reflector. 
157 -0.86 fg 868 -0.42 ve Q plates (9250 A) were used through JD 
162 | -0.82 7 877 | -0.50 fe 248524; thereafter M plates. A value of 
177 | -0.63 vg 878 | -0.42 zg -0.40 mag. was adopted for a Cet, the com- 
183 | -0.60 g 893 | -0.20 vg parison star, at this wave length. 
18? -0.65 g 902 | 0.03 g 
-0.73 fg 907 | 0.00 g 
199 | -0.62 f gia | 0.17 g TABLE 10 
208 | -0.63 € 926 | 0.33 fg 
209 | -0.63 zg 936 | 0.68 pf R Leonis at 8500 A 
220 | -0.65 954 | 0.75 f 
230 | 987 | 0:94 & || sp | Mag. [Quality] JD | Mag. |Quality 
232 | -0.20 g 964 | 0.95 f 242+ 242+ 
237 | 0.10 970 | 1.05 f 
363 | 1.24 975 | 1.46 te 7349 | 0.52 f | 7827 | 2.25 f 
444 | -0.63 | | 9133 | -0.57 fg 358 | 0.48 fe 8268 | 2.14 vE 
S61 | 0.48 841 | 2.25 vE 
The Allegheny refractor was used through JD pias oo § 658 2.25 & 
2428072; thereafter the Yerkes reflector. 408 0.28 t 871 2.36 € 
The comparison star was a Cet (-0.20 mag.). . & 888 1.95 g 
420 0.47 914] 1.38 VE 
431 0.57 f 924 1.49 ve 
TABLE 8 442 | 0.70 ft 940 | 1.57 pf 
o Ceti at 3700 A 468 | 1.14 pf 1s5 | 2.14 f 
482] 1.46 fe 190 | 1.43 pf 
gp | Mag. |quality | JD | Mag. |Quality || #64) 2-4 fe 198 | 1.72 fg 
242+ 242+ 513 | 2.14 217 | 1.55 fe 
525 2.15 i 218 1.42 ve 
7300 8199 fg 550 | 2.44 ve 231 | 1.72 f 
to 9.8 g 208 | 6.43 fg S42} 2.43 & 245 | 1.70 fr 
7350 209 | 6.50 fg 556] 2.57 fe 252] 1.57 g 
6115 6.92 g 220 6.56 fg 561 2.46 pf 275 1.30 fe 
128 | 6.29 vg 232 | 7.50 f 562 | 2.90 f 284] 0.57 ve 
129] 6.37 f 233 | 7.27 fe 566 | 2.23 pf 301 | 0.05 f 
138} 6.11 g 237 | 7.66 Pp 573 | 2.26 fg 3oz | 0.21 g 
147 | 6.08 f 460 | 5.65 fg 581] 2.14 fg 303 | 0.16 
157 | 6.07 f 468 | 5.07 fg 715] 0.28 vE 307 | 0.14 
162 6.10 fg 477 4.89 fg 727 0.47 C4 314 | -0.14 & 
6.08 g 501 | 4.70 fg 735 | 0.56 fg 317 | 0.01 g 
183 | 6.08 fg 524 | 5.35 g 758 | 0.99 fg 321] 0.05 pf 
187 6.03 fg 552 6.46 i 783 1.43 f 323 0.05 € 
194 6.05 fg 557 6.41 fg 806 1.95 f 331 0.47 € 
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TABLE 10 (continued) 
R Leonis at 8500 A 


TABLE 12 (continued) 
R Leonis at 9250 A and at 9100 A 


JD Mag. |Quality | JD Mag. Quality JD Mag. |Quality] JD Mag. | Quality 
242+ 242+ 242+ 
8334 0.47 g 8682 1.18 fe 8335 -0.71 fg 8618 -0.44 vE 
335 0.37 fg 691 1.69 é 337 -0.37 fe 614 -0.57 fg 
341 0.45 613 2.14 543 -0.15 fg 619 -0.41 
343 0.56 ve 820 2.08 4 352 -0.02 i 620 -0.56 fg 
350 0.85 g 845 1.91 ‘g 46C 1.53 pf 813 1.00 
352 0.87 g 848 1.69 g 471 1.54 vg 820 1.22 VE 
466 2.51 860 1.36 482 1.53 gl2 -0.66 vE 
471 2.63 ve 867 1.19 g 502 1.43 fe 31] -0.41 fe 
481 2.52 ve 878 1.00 fe 510 1.19 fg 959 | -0.42 
489 2.50 g 691 0.57 fe 518 1.15 + 9054 1.19 fe 
502 2.45 4 894 0.47 ve 614 -0.50 fg 
510 2.44 ? 904 0.37 & 
518 2.31 "4 gl2 0.15 fg The Allegheny refractor was used through JD 
560 1.42 fg 925 0.17 i 2427727; thereafter the Yerkes reflector. 
568 i637 ve 927 | -0.14 ? Q plates (9250 A) were used through JD 
592 0.63 fg 931 0.15 vg 2428616 except that on JD 2428614 both Q and 
597 0.37 & 955 0.47 vg M plates were used; thereafter M plates were 
609 0.15 g 959 0.57 f used. The value 1.00 was adopted for the 
618 0.15 g 965 0.70 vg magnitude of the comparison star, a Leo, at 
619 0.19 4 975 0.70 fg these wave lengths. 
630 0.26 992 0.98 
635 0.43 fg 9007 1.16 4 
642 0.50 é 010 1.30 fg TABLE 13 
660 0.63 g 030 1.72 fg 
670 1.00 f 054 1.95 fg R Aquilae at 8500 A 
677 0.99 vg 057 1.95 fg 
JD Mag. Quality JD Mag. Quality 
The Allegheny refractor was used through JD 242+ 242+ 
2426252; thereafter the Yerkes reflector. 
The comparison star was a Leo (1.00 mag. at 7238 | 2.09 f 7691] 4.42 p 
8500 A). 245| 2.17 vg 704| 4.62 Pp 
258 2.28 ve 714 4.52 € 
261] 2.26 vg 4.40 
TABLE 11 272] 2.99 pf 723 | 4.23 
283 2.57 644 1.72 
R Leonis at 3700 A 288} 2.55 ve 847] 1.57 p 
290] 2.70 ve 55 1.85 € 
JD Mag. |Quality | JD Mag. /Quality 291}; 2.71 ve 861] 1.95 f 
242+ 242+ 299 3.29 ft 871 1.95 g 
3.39 £ 680 | 2.08 fé 
6141 | 11.28 8330 7.07 fg 312] 3.42 f 888 | 2.33 fg 
166 | 12.60 fg 500 | 13.25 fg 325] 3.52 ve 915] 2.86 ve 
231 | 11.22 $12 | 12.37 fg 345] 3.85 fg 941] 3.72 
301 6.92 g 614 6.70 f 352] 3.66 fg 948 | 3.80 € 
307 7.00 fg 625 6.69 f 359] 3.84 f 964] 4.02 VE 
317 7.00 3.94 pf 980] 4.28 
516 1.87 f 8002 4.57 pf 
The fundamental comparison stars and their 517] 2.15 f 003} 4.55 pf 
adopted ultraviolet magnitudes were: HD S22} 1.86 pf 014] 4.60 pf 
64748 FO 6.33 and HD 64561 K5 8.33. 527] 1.58 Pp 034] 4.68 ve 
530] 1.69 O69 | 3.66 ve 
538 1.61 f leg 2.09 ve 
TABLE 12 548] 1.69 g 140] 2.09 g 
556; 1.87 f 146] 1.74 g 
R Leonis at 9250 A and at SIOO A 557 2.09 fg 157 1.73 ve 
560 2.07 fe 210 2.48 pf 
JD Mag. |Quality | JD Mag. 568} 1.94 2.49 
242+ 242+ S75 2.09 4 271 3.72 f 
2.16 fg 276) 4.63 p 
7372 | -1.07 6303 | -0.35 589} 2.44 277| 3.90 
376 | -1.38 305 | -0.44 ve 609; 2.40 293] 3.70 
430 | -0.93 g 314 | -0.30 g 628] 3.29 pf 294] 3.73 £ 
535 0.28 316 | -0.48 fg 3.52 298] 4.04 
558 1.43 Pp 317 | -0.24 fg 654| 4.01 pf 303] 4.99 fé 
561 0.28 vp 323 | -0.85 658} 4.00 308 | 4.10 
718 | -0.85 f 331 | -0.55 fe 679] 4.23 fg 315| 4.44 £ 
727 -1.01 f 334 -0.71 fg 680 4.20 324 4.15 pf 
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TABLE 13 (continued) 
R Aquilae at 8500 A 


TABLE 15 (continued) 
R Aquilae at 9250 A and at 9100 A 


JD Mag. |Quality | JD Mag. [Quality JD Mag. | Quality] JD Mag. | Quality 
242+ 242+ 242+ 242+ 
8331 3.93 f 8737 1.53 g 6441 1.63 f 8745 1.53 4 
351 3.01 f 740 1.59 g 442] 1.40 fg 750} 1.63 vE 
353 3.15 g 745 1.60 fg 444| 1.51 fg 761] 1.89 fg 
377 2.44 750 1.57 453 | 1.63 fg 618} 1.96 fz 
383 2.28 vE 761 1.66 g 459} 1.40 g 841] 2.34 fz 
395 2.27 € 813 2.94 g 460; 1.52 VE 863 | 3.20 fg 
409 2.25 veg 818 3.09 fg 471] 1.82 vE 864] 2.56 fe 
441 1.87 4 841 3.45 f 481 1.51 & 877 2.85 fg 
442 1.87 & 858 3.76 f 620 3.42 Pp 878 2.85 f 
444 1.85 fg B63 2.88 6 621} 3.35 fa 902] 3.53 f 
453 1.71 f 864 3.86 fg 660} 2.77 fg 903 | 5.16 fe 
457 1.78 vg 877 3.82 vE 737| 1.81 £ 9053 | 1.21 pf 
460 1.57 878 3.95 740) 1.63 055] 1.36 fe 
465 1.55 fé 902 4.21 g 
471 1.76 g 903 4.12 fg The Allegheny refractor was used through JD 
481 1.87 vg 976 3.37 vg 2427718; thereafter the Yerkes reflector. 
509 1.84 ve 986 3.33 g Q plates were used through JD 2428471, and 
610 4.77 Fg 998 3.10 f M plates thereafter. The magnitude adopted 
614 4.10 g 999 3.29 fg for the comparison star, a Aql, was 0.49. 
633 4.05 9011 2.75 fg 
643 3.90 g 029 2.15 vg 
660 3.40 ve 035 1.81 fg TABLE 16 
669 3.10 g 041 1.60 pf 
678 2.67 f 043 1.67 fg X Ophiuchi at 8500 A 
683 2.57 4 053 1.65 fg 
695 2.38 f 055 1.81  g JD Mag. Quality | JD Mag Quality 
709 2.38 vE 063 1.96 fg 242+ 242+ 
715 2.25 vg 069 1.84 fg 
077 2.06 fg 8465] 1.94 fe 8737 | 1.96 
471} 2.11 2.23 
The Allegheny refractor was used through JD 481] 1.63 € 761) 2.04 g 
2428069; thereafter the Yerkes reflector. 483 | 1.76 f% 813] 1.66 ve 
The comparison stars were: a Aql (mag. 0.43,|| 490] 2.09 g 816} 1.79 & 
JD 7238 through 7740 and 8069 through 8610); 502| 1.87 € 841] 1.88 & 
a Cyg (mag. 1.00, JD 7844 through 8043) and 509} 1.54 pf 858 | 2.33 & 
HD 164406 (mag. 3.10, JD 8614 and after). 610] 3.52 g 864] 2.57 ve 
614] 3.53 g 877 | 2.33 ve 
633 | 3.47 fg 3.73 vE 
TABLE 14 643 | 3.53 986] 3.63 fg 
660} 3.40 998] 3.40 pf 
R Aquilae at 3700 A 669} 3.51 4 999] 3.53 fg 
678 | 3.33 fe 
JD Mag. alit JD Mag. alit 683 | 3.26 & 29] 2.8 fe 
242+ 242+ 695 2.94 035 2.91 fg 
709| 2.62 4 041} 2.80 pf 
8303 | 12.76 fg |8444 | 8.00 fg 715] 2.67 6 063 | 2.75 fe 
319 | 12.79 fg 622 | 12.81 fg 069] 2.73 fg 
334 | 12.86 fg 636 | 12.71 fé 
441 8.28 fg 647 | 12.67 fg The Yerkes reflector was used throughout. 
740 8.28 fg The comparison stars were: a Aql (0.43 mag.) 
through JD 2428610; thereafter HD 164406 


The fundamental comparison stars and their 
adopted ultraviolet magnitudes were: 
178162 B9 9.00 and HD 178216 GO 10.10 


HD 


(3.10 mag.) 


TABLE 17 
TABLE 15 X Ophiuchi at 3700 A 
R Aquilae at 9250 A and at 9100 A | 
JD | Mag. |Quality | JD | Mag. (Quality |) e481] 7.62 fg | 8647] 9.91 g 
242+ 242+ 490 | 7.98 665 | 9.61 
636 | 9.90 g 695 | 9.81 g 
e fundamental comparison stars and their 
| 3.42 353 | 2.98 KO 8.44; HD 172 AL 7.117 HD 172588 
317 | 4.00 pf 409 | 2. fe ia. 
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TABLE 18 TABLE 19 (continued) 
X Ophiuchi at 9250 A and at 9100 A R Aquarii at 8500 A 
JD Mag. |Quality | JD Mag. |Quality JD Mag. | Quality; JD Mag. | Quality 
242+ 242+ 242+ 242+ 
8471 1.63 6812 | 1.44 pf 9055} 4.69 9078 | 4.40 f 
481 1.53 ve 813 | 1.32 fg O72) 4.75 fg 1l2 4.33 fg 
4681 1.37 fg 618 | 1.15 fg 073 | 4.97 4.06 & 
483 1.68 fg 841 | 1.25 g 
660 S77 f 986 2.63 g The Yerkes reflector was used throughout. A 
678 2.03 fg 993 | 2.42 fg field comparison star, HD 222643 K2, was 


Q plates were used through JD 2428481; there- 
after M plates —_ A). The comparison 
star was a Aql 


mag. 0.49). 


TABLE 19 
R Aquarii at 8500 A 


assumed to be of magnitude 3.80 at 8500 A. 


JD Mag. [Quality | JD | Mag. |Quality 

242+ 242+ 

8464 | 2.46 fg  |6849 | 1.76 g 
465 | 2.50 fg 861 | 2.08 fg 
707 | 4.43 f 877 | 2.42 g 
708 | 4.78 vg 878 | 2.51 ve 
715 | 4.94 f 893 | 2.28 ve 
735 | 3.80 fe 902 | 2.66 
751 | 3.61 fg 907 | 2.78 f 
fe 925 | 2.97 vg 
814 | 2.30 ve |9035 | 4.33 f 
844 | 1.87 f 041 | 4.35 g 


TABLE 20 

R Aquarii at 9250 A and at 9100 A 
JD Mag. Quality | JD Mag. Quality 
242+ 242+ 
6465 | 1.37 8861 1.32 
712 | 3.60 877 | 1.59 pf 
714| 3.17 878 | 1.98 f 
715 | 3.47 fg 902 | 1.80 
735 | 2.63 f 9041 | 3.37 f 
751] 2.81 vg 078 | 3.41 fg 
763 | 1.80 


The Yerkes reflector was used throughout. WM 
plates were used except on JD 2428465. The 
comparison star, HD 222643 K2, was assumed 
to have a magnitude of 3.60 at 9100 A. 


| 
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T Cephei.—The curve (Fig. 1) at 8500 A is very similar to the 
visual-curve, but the average magnitude ranges from five maxima 
and four minima are 2.56 and 4.63, respectively. All infrared max- 
ima attain about the same value, while the one visual maximum 
which is conspicuously lower is rather broader, so that the area under 
the curve from a magnitude below maximum is of the same order as 
for other maxima. Two minima are definitely lower visually and cor- 
respondingly at 8500 A. Over two maxima and three minima the 
ultraviolet-curve practically parallels the visual, the ranges being 
4.75 and 4.70. The photographic range by Miiller is less than the 
visual in the ratio 3.22: 3.67 over two cycles, probably because band 
absorption in the visual more than compensates for the temperature 
change. The ranges at g100 A and 9250 A are 2.0 and 1.9, as com- 
pared with 2.6 at 8500 A. 

T Cep obviously does not everywhere radiate like a black star. 
Since the ultraviolet-, photographic-,and visual-curves are practically 
parallel, some of the color indices will differ only by a constant giving 
a 45° line not passing through the origin. The ultraviolet-minus-8500A 
and the ultraviolet-minus-gioo A indices, however, give a slope of 
index 0.92, passing through the origin. The theoretical value of the 
slope is 0.94. The temperatures from the various indices derived 
from formulae (2) are as shown in the accompanying table. If 2975° 


INDICES TEMPERATURES 
WAVE 
LENGTHS 

Max. Min. Max. Min. 
5500-8500. ... 4.098 (5 7.15 (4) 1700° 1240° 
3700-8500. . 5.87 (2) | 8.12 (3) 2975 2320 
4200-8500......| 6.85 (2) | 8.02 (2) 2160 1900 
3700-QI00..... 6.50(1) | 8.95 (1) 2860 2235 
3700-9300..... 6.40 (1) | 9.10 (1) 2920 2230 


at maximum and 2320° at minimum are chosen as most probable, 
the visual magnitude range for a black body is 2.72, corresponding to 
4.87 observed and indicating a differential absorption between maxi- 
mum and minimum of 2.15 mag., provided there are no changes in 
radius and no band absorption in the ultraviolet and at 8500 A. 

In the plot of visual and 8500 A magnitudes (Fig. 4) the points for 


| | 
| | 
| 
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the ascending branches of the curves lie near the slope of index 2.00, 
while those of the descending branches indicate 1.83, averaging 1.92. 
Since the infrared maxima are flatter, particularly following the 
visual maxima, both branches are perhaps best represented by the 
slope-indices 2.33 and 1.78 toward maximum and minimun, re- 
spectively, for the ascending branches, and 2.44 and 1.60 for the 
descending branches. Thus, toward the minima the slope-indices 
are nearest the theoretical value of 1.55, indicating little change in 
absorption, but a nearly uniform graying; but toward the maxima 
the absorption lessens more or less linearly, but not always in phase 
with the temperature change. The ultraviolet-minus-8500 A slope- 
index of 2.1 is near the theoretical value of 2.3. The photographic- 
minus-8500 A index of 1.75, compared to the theoretical value of 2.0, 
might indicate an increase of diameter toward maximum, but prob- 
ably differential band absorption also obtains in the photographic 
region. 

R Cassiopeiae. The infrared maxima are rather less variable than 
the visual maxima, and the average magnitude ranges from five 
maxima and four minima are 3.07 and 5.65, respectively (Fig. 1). 
As in T Cep, the photographic range (Miiller) is less than the visual 
in the ratio 4.80: 5.40, while the ultraviolet is greater than the visual 
in the ratio 7.22:5.80. The ranges at g1oo A and at 9250 A are 2.6 
and 2.9, to 3.6 at 8500 A. The plot of the ultraviolet-minus-8500 A 
and visual-minus-8500 A indices does not pass through the origin, 
and the slope-index is 1.5, compared to the theoretical value of 2.4; 
but the line formed by the ultraviolet-minus-8500 A and the ultra- 
violet-minus-g1oo A indices passes through the origin with a slope 


INDICES TEMPERATURES 
WAVE 
LENGTHS 

Max. Min. Max. Min. 
5500-8500......| 5.49(5) | 8.07 (4) 1560° I110° 
3700-8500......| 6.90(3) | 11.2 (2) 2050 1800 
4200-8500......}| 6.60(1) | 9.10 (1) 2240 1715 
2900-OL00......| 7-084) x20 (x) 2700 1750 
3700-9300...... 7 | 12.05 (x) 2600 1800 


equal to the theoretical value of 0.94. Indices from these spectral 
regions should therefore give the best temperature estimates. Choos- 


sapnztusew jensta juareddy o nby y pur (a2ogn) [by y 


‘(mojaq) oa] pur X (aa0gn “74324 07 dad J, pur sey xaput adojs [ensta 


I— ° I £ 
see 
8 — ‘ 
£— 
Guan ‘ 
° 
38 
° 
aee 


—z 6 


gio— zo 


<4 


+ 


| 
| 
4 
| 
| 
| 
| 
{ 
| | 
a j 
} = 
| 
i 
N 
| 
| 
| — 4 
| 
in 
Got 8, —## 
&°¢ 
re) 
eo. 
38 
eo S20 
e 
nm 
° — 
| 
| 
| 


76 CHARLES HETZLER 


ing 2650°-1800° as most probable, the black-body visual range would 
be 5.10 mag., which corresponds to 5.90 observed, showing a difler- 
ential absorption of 0.80 mag. 

The visual-minus-8500 A plot (Fig. 4) gives 1.95 as the slope- 
index of the ascending branches and 1.80 as the slope of the descend- 
ing branches of the curves, the average of 1.87 being close to that for 
T Cep. But there is no indication of curvature, as with that star. 
The photographic-minus-8500 A slope-index of 1.85 and the ultra- 
violet-minus-8500 A index of 2.1 (in exact agreement with T 
Cep) are fairly near the theoretical values of 2.0 and 2.3; and possi- 
bly only slight, if any, radial changes are indicated. 

R Leonis.— The average magnitude ranges from five maxima and 
five minima at 8500 A and in the visual region are 2.16 and 3.83, the 
curves being very similar (Fig. 2). Over two cycles the ultraviolet 
and visual ranges are in the ratio 6.05:4.22, while the ranges at 
gtoo A and at g250 A average 1.70 and 1.52, to 2.00 at 8500 A. The 
ultraviolet-minus-8500 A index plotted against the visual-minus- 
8500 A index intercepts below the origin and has a slope of 1.9, com- 
pared to the theoretical value of 2.4; the line formed from the ultra- 
violet-minus-8500 A and the ultraviolet-minus-gt1oo A indices passes 
through the origin with a slope of 0.93, in agreement with T Cep, 
R Cas, and the theoretical value. The temperatures from the various 
indices are given in the accompanying table. The temperatures of 


INDICES TEMPERATURES 
LENGTHS 
Max. Min. Max. Min. 
5500-8500......|  5.79(5) | 7.30 (5) 1490" 1220" 
3700-8500..... 6.75 (2) | 10.60 (2) 2070 1860 
3700-9300......| 7.30(2) | 11.45 (1) 2050 1850 


2670° at maximum and 1860° at minimum would give a black-body 
visual range of 4.18, in practically exact agreement with the corre- 
sponding value of 4.20 observed over the same cycles as for the tem- 
peratures, the differential absorption being zero. 

The plot of visual and 8500 A magnitudes (Fig. 4) for both ascend- 
ing and descending branches gives the slope 1.82 without much scat- 
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ter or definite indication of curvature: the ultraviolet-minus-8500 A 
slope is 2.71; and the ultraviolet-minus-9250 A slope is 3.2. These 
are all greater than the theoretical values of 1.55, 2.3, and 2.5; and 
our data would indicate that, since there is no differential absorption, 
but merely a uniform graying in the visual region, the star expands 
toward minimum, the amount being readily derived by equation (4), 
if absolute magnitudes are determined. 

R Aquilae-The average magnitude ranges at or near maxima 
and minima over seven cycles are 4.76 visually and 2.31 at 8500 A 
(Fig. 2). Over two cycles the ultraviolet range is less than the visual 
in the ratio 4.55:4.81, while for three maxima and two minima at 
g1oo A and for one cycle at 9250 A the ranges are 1.92 and 1.86, re- 
spectively, to 2.40 at 8500 A. The visual curve practically par- 
allels the ultraviolet curve, but the ultraviolet-minus-8500 A against 
the ultraviolet-minus-gioo A slope passing close to the points and 
through the origin is again the theoretical o.94. The temperatures 
from the various indices are as shown in the accompanying table. 


INDICES TEMPERATURES 
WAVE 
LENGTHS 

Max. Min. Max. Min. 
5500-8500..... 4.59 (7) 7.00 (7) 1820° 1265° 
3700-8500......] 6.45 (2) | 8.70(1) 2775 2200 
3700-QI00..... 6.60 (2) 9.40 (1) 2830 2150 
3700-9300......| 6.50 (1) g.20 (1) 2890 2200 


Adopting 2800° at maximum and 2200° at minimum the black- 
body visual range would be 2.76 mag., as against 4.81 observed over 
the same cycles as for the temperatures, the differential absorption 
being 2.05 mag. in the visual region, if the others have none. 

The visual-minus-8500 A line (Fig. 4) has an average slope of 
1.94; but toward maximum the slope is more nearly 2.4 for both 
branches, while toward minimum it is about 1.3 for the ascending 
branches and 2.0 for the descending. The ultraviolet-minus-8500 A 
slope is 2.0, and the ultraviolet-minus-g1oo A slope is 2.3. Probably 
there are radial changes as well as absorption changes; and since the 
latter two slopes are less than the theoretical values of 2.3 and 2.5, 
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respectively, there are indications of an increase in diameter toward 
maximum, as was already suggested for T Cep and R Cas, while R 
Leo apparently expands toward minimum. Probably the present 
data are insufficient for a final decision, but they suggest that there 
are striking differences in the infrared variables. 

The four variables thus far considered——T Cep, R Cas, R Leo, and 
R Aql—are all of similar spectra, M5~M8. Plots of visual and ultra- 
violet magnitudes against infrared magnitudes give straight lines, 
particularly toward minima. This feature can be used in considering 
three further stars of similar type-—-o Cet, X Oph, and R Aqu—of 
which the visual and ultraviolet magnitudes toward minima are af- 
fected by the presence of known or suspected companions. Extra- 
polation of the straight lines beyond the points where they become 
curved because of the companions gives the true visual and ultra- 
violet magnitudes of the variables for any particular total magnitude; 
hence, it also gives the visual and ultraviolet magnitudes, or the 
colors, of the companions. 

o Ceti.— Many of the maxima (Fig. 3) show a delayed drop of the 
infrared points, a feature already noticeable, though less conspicuous, 
in T Cep. The infrared maxima overlie the projection of the visual- 
curve a magnitude or so below maxima, over which the visual peak is 
usually, but not always, displaced to the left.’ The visual peak 
shows more fluctuation in range, as well as in position, than the in- 
frared. The effect upon the visual-infrared plot (Fig. 4) is an up- 
bending toward maximum. There is also a decided upbending to- 
ward minima, which is interpreted as being caused by the known 
companion, since it does not occur in the stars already considered. 
Extending the line toward minimum of slope 1.60, which presumably 
is affected neither by the differential maxima nor by the companion, 
the visual magnitude of the variable for total visual magnitude 9.2 is 
9.7, giving 10.3 as the visual magnitude of the companion. This is in 
agreement with the average estimates of 9.8 by Joy,° of 9.6 by Pettit 
and Nicholson,’ and of 11.1 by Van Biesbroeck.* 

In the ultraviolet we use our observations obtained at the 1933 

5 Pub. A.A.S., 9, 124, 1937. 7Ap. J., 78, 327, 1933. 

6 Ap. J., 63, 281, 1926. 8 Ap. J., 82, 84, 1935. 
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minimum at the Allegheny Observatory’ with the same filter as that 
used over two maxima at the Yerkes Observatory. Extending the 
ultraviolet-minus-8500 A line of slope 2.0, we obtain 11.3 as the 
ultraviolet magnitude of Mira, corresponding to a total ultraviolet 
magnitude of 9.8. Hence, the ultraviolet magnitude of the compan- 
ion, based upon an adopted ultraviolet magnitude of 9.90 for the 
A3star BD —3° 357, is 10.1. Adopting the values in the following table 
for the magnitudes of the two comparison stars used at Yerkes and 
at Allegheny, the results for o Ceti A and B are: 


MAGNITUDE 
STAR 

Ultraviolet |Photographic Visual 
BD—3°357 (A3)...... 9.90 9.85 9.77 
BD—3°355 (Ma). 11.34 9.65* 9.08 
o Ceti (A+B) 9.8 9.75 9.2 
o CetiA.. 48.3 10.7 6:7 


* Derived from the preceding star and differing from the HD esti- 
mate of 1o.4 assumed at Allegheny. 


These magnitudes are all interrelated; and the values are fairly con- 
sistent, except possibly in the photographic region, where the mag- 
nitudes have been rederived from the few observations at Allegheny.’ 
In this region the A and B components were nearly equal, so that 
10.3 and 10.7 are about the limits corresponding to total magnitude 
9-75. 

The temperatures from the various indices, after correction for the 
effect of the companion, are as shown in the accompanying table. 


INDICES TEMPERATURES 
WAVE 
LENGTHS 
Max. Min. Max. Min. 
5500-8500. ..... 4.70 (6) 7.6 (3) 1800° 1175° 
3700-8500...... 6.38 (2) 9.3 (1) 2800 2035 
3700-9300 6.00 (1) 10.9 (1) 3000 1915 
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Temperatures of 2800°-2035° would give a visual-magnitude range 
of 3.3. The corrected visual range from six maxima and five minima 
is 5.66, as compared to 2.52 at 8500 A. The visual range correspond- 
ing to these temperatures is about 6.1, but the various maxima are 
extremely variable and the indicated differential absorption of 2.8 
mag. is doubtless also variable. The very steep slope of 3.7 of the 
visual-minus-8500 A line toward maxima probably indicates a rapid- 
ly changing absorption; but the medial slope of 1.95 is similar to 
that of the other stars. The star o Cet most resembles T Cep; and, 
except near maxima, the various slopes of the two stars are practi- 
cally identical. 

NX Ophiuchi.—This variable has a companion of K-type spectrum 
for which Pettit and Nicholson deduced a visual magnitude of 8.9.7 
The points at 8500 A (Fig. 3) have almost as great a range between 
maxima and minima as the visual, and the infrared-curve may rise 
while the visual-curve maintains a constant minimum. Extending 
the line of slope-index 1.80 of the visual-minus-8500 A plot (Fig. 4), 
the visual magnitude of the variable is 1o.o (in general agreement 
with Van Biesbroeck, as against Pettit and Nicholson’) for total 
apparent magnitude 8.6, giving 8.9 as the visual magnitude of the 
companion. The ultraviolet-minus-8500 A line of slope 2.1 gives 
11.2 for the variable, corresponding to a total apparent magnitude 
of g.6, or to 9.9 for the ultraviolet magnitude of the companion. The 
temperatures from the various corrected indices are given in the ac- 
companying table. Possibly 2g00° 2400” is fairly close to the actual 


INDICES TEMPERATURES 
Wave 
LENGTH 
Max. Min. Max. Min. 
5500-8500... | 1666° 1370 
3700-8500... 6.0 (1) awn) 2920 2410 
3700-QI00..... 6.1 (1) 8.8 (1) 3000 2250 


temperature range of the variable. The corrected visual range corre- 
sponding to the infrared range (8500 A) of 1.7 mag. is 3.0 mag., while 
for a black body it would be 2.1 mag., indicating a differential ab- 
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sorption of 0.9 mag. The ranges at 8500 A and gtoo A are 1.8 and 
I.1, respectively, from two cycles. 

R Aquarii.— Our observations definitely indicate the presence of a 
companion by the strong curvature of the visual-minus-8500 A plot 
(Fig. 4) toward minimum. Extending the line of slope 1.92, the vis- 
ual magnitude of the variable corresponding to a total visual mag- 
nitude of 11.2 is 12.8, making the visual magnitude of the companion 
11.5. This cannot be far off, although Pettit and Nicholson ob- 
tained 10.0.7. Unfortunately, few ultraviolet observations of this 
star were obtained; hence the color remains unknown. The magni- 
tude range at gioo A is 2.2 mag., as compared to 2.5 mag. at 8500 A. 


YERKES OBSERVATORY 
September 1939 


CONTRIBUTIONS TOWARD A PHYSICAL MODEL 
OF y CASSIOPEIAE 


RALPH B. BALDWIN 


ABSTRACT 

The light-curve and the temperature-curve of y Cassiopeiae are correlated in a 
peculiar manner. An explanation is that the photosphere of the star is variable in size 
and that a reduced color temperature depends on an increase in the extent of the ra- 
diating surface relative to the available energy. 

The law of conservation of angular momentum may be used to explain these corre- 
lations. The broad absorption lines are observed to change in width; and if this is 
interpreted as a change in the velocity of rotation, we immediately have relative 
radii of the star for all dates covered by measures of the absorption width. The original 
absolute magnitude and the light-curve permit computation of a temperature-curve, 
which is found to be in agreement with the observed temperature-curve. 

This principle allows a determination of the radii of the layers o1 emission. It is found 
that the atmosphere is stratified and that it pulsates in a manner similar to the photo- 
sphere but with a moderate lag. 

Equivalent widths of certain sharp absorption lines of helium and hydrogen have 
been measured, and the method of Struve and Wurm has been applied to give the 
distance from the photosphere to the “shell” in which the lines were produced. The 
results are consistent with the radii determined from the broad absorption and emission 
lines; the sharp absorption lines are produced in the outermost layers. 

Variations in the emission intensity ratios of the two components of the lines of ten 
elements are shown to exist and are discussed. 


The 216 spectrograms of y Cassiopeiae on which this study is 
based cover the period 1926-1938. Most of the plates are on Process 
emulsion and were taken with the Michigan one-prism spectrograph ; 
the dispersion is 40 A/mm at Hy. Because of the unusual trans- 
parency of this glass prism the ultraviolet spectrum was photo- 
graphed beyond the limit of the Balmer series. 


THE LUMINOSITY 
Many observers have discussed the problem of the luminosities of 
stars of classes B and Be, and they generally agree that the Be stars 
are slightly brighter. Mean absolute magnitudes of Bo stars range 
from —2™5 to —3™4, according to the work consulted. However, 
space absorption is not considered in the earlier estimates. Including 
the effect of space absorption, Stebbins, Huffer, and Whitford" tind 
a mean absolute magnitude of Bo stars of —3™9. KE. G. Williams” 
1 Ap. J., 90, 209, 1939. 2A p. J., 83, 305, 1936; Mt. W. Contr., 541, 29, 1936. 
82 
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value of the absolute magnitude of y Cassiopeiae, derived from inter- 
stellar line intensities, is —2™5. The original color excess' was 
+orr, at least part of which is an intrinsic property of the star; the 
remainder is due to space reddening. The dimming of any star by con- 
tinuous interstellar absorption (A ~~) is 7£ on the visual-magnitude 
scale." Consequently, the uncorrected absolute magnitudes found for 
y Cassiopeiae must be decreased by a quantity between o™oo and 
o"77. A direct trigonometric parallax’ yields —o™03, which seems 
too faint. In the following discussion the absolute magnitude will 
be considered as — 2's. 


THE BROAD ABSORPTION LINES 


HISTORICAL 

In 1932 Struve! showed that the broad absorption wings which 
flank the emission lines of hydrogen in Be stars have the contours of 
normal absorption lines which have been abnormally broadened by 
rapid axial rotation. Struve did not specify what actual velocities 
were present for each star, but showed only that the wing contours 
were of the form to be expected if rotational velocities were large. 
If we make any reasonable assumption for the mass and density of 
a Boe star, the maximum velocity of rotation which could obtain 
is of the order of 300 km sec, derived from Jeans’s criterion for 
liquid bodies.® 

It is not established that the laws found for liquid bodies are ap- 
plicable to the gaseous and compressible materials of stars. It may 
be that even higher rotational velocities would not disturb the sta- 
bility of the average Boe star where centrifugal force does not become 
equal to gravity until the velocity is greater than 500 km/sec. 

The first measurable part of the 178 wing of y Cassiopeiae is dis- 
placed 16 A from the center on microphotometer tracings;* and if 
this were to be interpreted as a Doppler effect, it would correspond 
to nearly 1000 km sec. Obviously then, a measure of the displace- 
ment of any definite point in the wing of a line will not, by itself, 
give the rotational velocity but will give a velocity which is propor- 
tional to, and larger than, the rotational velocity. 

3 Schlesinger, General Catalogue of Stellar Parallaxes, 6, 1935. 


4 Ap. J., 76, 300, 1932. 5 Astronomy and Cosmogony, 198, 1928. 
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THE MEASURES 

In general, the structure at any strong line is that of a double re- 
versal. The continuous spectrum is weakened by many very broad 
absorption lines. An emission line approximately bisects each of the 
the broad absorption lines and is in turn divided into two compo- 
nents by a sharp absorption line (Fig. 1). 

As described by several authors,’ the spectrum of y Cassiopeiae 
was relatively constant from 1gtr to 1928, when slight changes be- 
gan. The star has passed through three increasingly violent cycles 


(5) 


Fic. 1.-Ideal cross-section through a complex line, showing a broad absorption 
line bisected by an emission structure, and the latter, in turn, divided by a central 
“absorption” dip. The indicated points are the places at which the seven velocity- 
curves of Fig. 2 were measured. 


of change since 1928'??4 and is at present (1939) beginning a fourth. 
In each cycle the two emission components varied oppositely in in- 
tensity (V/R variation), while the emission and absorption lines 
grew narrower; and the components of the former even blended in 
1934 and 1937 (single-line stages). During the latter cycle the color 
temperature first dropped 5500° K from 13,500 K and then returned 


6 Curtiss, Pub. Obs. U. of Michigan, 2, 1, 1916. 
7 Ibid., 3, 1, 1923. 

8 Tbid., p. 10. 

9 Tbid., p. 256. 

1 Lockyer, M.N., 93, 362, 1933. 

™ Cleminshaw, Ap. J., 83, 495, 1936. 


6 Ibid., p. 670. 

'7 McLaughlin, Ap. J., 84, 
8 Thid., 85, 181, 1937. 

99 Heard, J.R.A.S. Canada, 32, 


235, 1930. 


353, 1935. 


2 Baldwin, Ap. J., 87, 573, 1938. 
3 [bid., 89, 235, 1939. 

"4 Tbid., go, 204, 1939. 

5 Lockyer, M.N., 96, 2, 1935. 


24 Cherrington, Perkins Obs. Contr., No. 


20 Heard, Ap. J., 81, 341, 1935 


2" Cherrington, Pop. Astr., 46, 183, 1938. 


22 Cherrington, Pub. A.S.P., 49, 91, 1037. 
23 Perkins Obs. Contr., No. 11, 1938. 


12, 1939. 
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nearly to normal." The light fluctuated in an irregular manner, 
although minimum brightness occurred at minimum color tempera- 
ture. ‘The number of visible emission lines increased markedly dur- 
ing this cycle. 

In none of the earlier investigations had anyone studied the be- 
havior of the broad underlying absorption lines. Inspection of the 
plates disclosed that throughout the entire series there were broad 
lines of hydrogen, helium, and ionized iron, while ionized magnesium, 
ionized silicon, and, surprisingly, ionized nickel appeared in absorp- 
tion near the single-line phase of the 1937 cycle. 

Complete measures were made for width of Hé, He1 4026, Fe 
4233, Siu 3863, Mg ut 4481, and Ni 11 4067. The point of setting 
was, of course, rather arbitrary; but an attempt was made to measure 
the first point where the wing could be distinguished from the con- 
tinuous spectrum except at //6. At this line there seemed to be a 
more or less definite break in the contour of each wing, a sudden 
dip; and the settings were made on these points and not in the ex- 
treme wings. As a result, //6 is usually given as narrower than Fe 1 
4233, although it is really broader. These features are actually parts 
of the 176 wings and are not to be attributed to V m1 4097 and NV in 
4103. If y Cassiopeiae is in class Bo, these lines will be present, but 
they should not show the measured changes of radial velocity. 

The individual measures for each line of all elements scattered 
considerably from plate to plate, but this is only what would be ex- 
pected in view of the diffuseness of the absorption. Nevertheless, 
when the measures from the two hundred odd plates are compared, 
it is seen immediately that definite trends are present. Prior to JD 
2420900 in 1932, the absorption widths at //6 were apparently con- 
stant. There are not enough plates covering the 1928 cycle to show 
possible small changes in absorption width. 

25 Greaves, Observatory, 62, 252, 1930. 

26 MLN., 98, 435, 1938; Observatory, 61, 118, 1938. Temperatures by Greenwich 
observers. I have adjusted their scale. 

2? Hunter, Martin, Observatory, 62, 237, 1930. 


C.R., 203, 1320, 1930. 
29 Observatory, 60, 320, 1937. 3! Williams, University of Michigan, unpublished. 


3° Observatory, 62, 76, 1939. 3? Huffer, Ap. J., 89, 139, 1939. 


86 RALPH B. BALDWIN 


Plates taken after the middle of 1932 tell an entirely different 
story. The broad absorption line, 17/6, gradually sharpened until in 
the latter part of 1934 it was more than 200 km, sec narrower than 
formerly. This narrowing was essentially in phase with the narrow- 
ing of the emission lines in the 1934 cycle. The absorption then wid- 
ened at a more rapid rate than that of the narrowing and came ap- 
proximately to the original width at the beginning of 1936. This 
was followed by a rapid decline in absorption width to a value 400 
km sec less than normal during the 1937 cycle. 


+ +600 4 
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+200 4 
L-200 4 
+ - 300 a 
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Fic. 2.—Velocity-curves for absorption and emission structures at //6. (1) and (7) 
are the curves for the red and violet absorption wings, respectively; (2) and (6) are for 
the red and violet emission edges; (3) and (5) are for the red and violet emission maxi 
ma; (4) is the central ‘‘absorption’”’ velocity-curve. 


In an earlier paper'’ it was pointed out that the red emission 
velocities were slightly out of phase with the violet emission veloci- 
ties. The absorption edges were also out of phase, but cach edge was 
in phase with the emission-curve for the corresponding side (Fig. 2). 
The narrowest absorption came at the time of the single-line emis- 
sion state; and the absorption widened as the emission widened, un- 
til in 1938 the absorption was even broader than it was originally. 

Complete measures of the lines He 1 4026 and Fe 11 4233 indicate 
a slight, but real, decrease in velocity width at the 1934 cycle and a 
much more prominent narrow stage during the 1937 cycle. The 
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Meg u, Siu, and Ni 11 lines were not visible until 1936 but thereafter 
were unanimous in showing a narrowing and subsequent widening. 
The mean curves of Figure 3 are simply plots of the average absolute 
values of the velocity displacements of each edge of each line for each 
plate. Even though the edge curves of //6 are slightly out of phase, 
they have been combined in this manner for comparison purposes. 
The //6 curve in the 1937 phase passes through its changes slightly 
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lic. 3.-- Mean velocity half-widths of six broad absorption lines 


earlier than do the curves of the other elements, and the He 1 maxi- 
mum width in 1936 comes a bit later. Otherwise, the six mean curves 
are different only in the amount of narrowing. 

A complication is found for //6 (and other 7 lines), namely, that 
there is a V/R variation of relative absorption-wing intensities 
which is essentially in phase with the V,/R curves for the 16 emis- 
sion (Fig. 4). The absorption is strengthened on the side adjacent 
to the stronger emission component, an effect similar to one discov- 
ered in the spectrum of 25 Orionis by Miss Dodson.*} The absorp- 
tion edge on the side of the weaker emission component apparently is 
not affected until the emission components reverse intensities. This 


33 Ap. J., 84, 180, 1936. 


4233 Fe = 
100 
po 
40% He, 100 
38635 
. M 
= 200 
44BI 100 
| 200 
100 
| 
| 
| 


88 RALPH B. BALDWIN 


additional absorption may be at least a partial explanation of the 
lack of agreement in phase of the absorption edges of //6. ‘The nor- 
mal state at this line is one in which the violet absorption component 
is slightly stronger than the red. 

A general weakening of the absorption lines could produce an 
apparent narrowing. Such a diminution of strength is present during 
1937 but not in 1934. It is suggestive of general absorption in the 
extensive atmosphere but does not appear large enough to account 
for the observed changes of line widths. 


3900 3500 (000 7000 7500 8000 800 


Fic. 4.Emission V/R curve (top) and absorption-wing V/R curve (bottom). 
Note subcycle in latter curve near JD 2429100. The emission V/R curve is not well 
defined near this time. 


If the widths of the diffuse absorption lines are caused by Doppler 
effect, they will vary linearly with wave length; whereas Stark effect 
and other widening influences are not linear with wave length. 
Stark effect is symmetrical for hydrogen only and scarcely affects 
the lines of the ionized metals. As these lines are also very broad and 
symmetrical, it follows that the widening process is not Stark effect. 

On two plates chosen at random near the single-line stage, meas- 
ures were made on the Balmer lines as far as Hu. The results indi- 
cate a constant velocity width. The essential similarity of phenom- 
ena in the broad metallic absorptions and in the hydrogen wings sug- 
gests a rapid rotation as the broadening agent. 


INTERPRETATIONS 


In earlier paragraphs it was shown that the great widths of the 
metallic absorption lines, and probably those of hydrogen and heli- 
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um, are due to rotation of the star. A change of line width implies 
either a change in the rotational velocity or an increased general ab- 
sorption in the shell. The latter mechanism was undoubtedly pres- 
ent, particularly in 1937. In 1934 a sharp decrease in line width oc- 
curred practically independent of an increase in the general absorp- 
tion. In addition, the absorption and emission widths are closely re- 
lated. It seems probable that the rotational interpretation is the 
more logical hypothesis. In this view, a decrease in the rate of rota- 
tion is dependent upon an increase in the radius. This is a necessary 
result of the principle of conservation of angular momentum in a 
closed rotating system. An increase in the line-width, then, is indica- 
tive of a contraction of the photosphere. As these fluctuations in 
line-width do occur, we infer a correspondingly variable apparent 
surface. 

Fortunately, we can check this line of reasoning observationally. 
A brightness range of about one magnitude was measured photo- 
clectrically..? Correlated with the light-curve is a temperature- 
curves 5! which at one point was 5500° K lower than normal. An 
explanation of this temperature drop in terms of superimposed 
Paschen continuous emission*! has been shown to be untenable."’ 
The only theory which seems capable of explaining simultaneous 
changes of light and heat is one in which the photosphere itself is al- 
ternately large and small. The necessary range of photospheric radii 
is less than 2. 

Kosirev* offers a logical explanation of the anomalous conditions 
obtaining in stars with extended atmospheres where the color tem- 
perature is very much lower than the ionization temperature. He 
shows that the energy distribution from a star with a deep, tenuous 
atmosphere will depart greatly from that of a black body, in the 
sense that the star will be rich in ultraviolet radiations. We can 
think of Kosirev’s energy distribution as coming from a raised 
photosphere, although the usual concept of a photosphere is not real- 
ly applicable. The ‘‘photosphere”’ in the terminology of this paper 
simply refers to some average layer within the greatly extended lay- 
ers of the star which produce continuous absorption. 

If we disregard such effects as general absorption would produce 
on the broad lines, and interpret the measured widths as pure Dop- 

34 Woolley, A.NV., 95, 115, 1934. 35 M.N., 94, 430, 1934. 
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pler effect of rotation, we may then apply the law of conservation 
of angular momentum. 

The true velocity of rotation (v sin 7) cannot be determined by 
measures of apparent width. However, the theoretical contours of 
absorption lines for all rapid rotational velocities are similar; thus, 
the measures of apparent widths must be nearly proportional to the 
actual rotational velocities, and in a computation involving relative 
velocities we may substitute the measured line-widths. 

In any closed rotating system angular momentum must be con- 
served according to the formula. 


RV = Constant. (1) 


Thus, if the effective radius is increased, the velocity of rotation must 
decrease proportionally. Chandrasekhar* has suggested that angu- 
lar momentum is conserved in the atmospheric layers of an individ- 
ual, nonpulsating star; i.e., the outer layers rotate more slowly. 

The measured line-widths may now be substituted into equation 
(1), and relative radii may be determined for every date from 1926 
to 1938. 

If, for any time, the actual photospheric radius can be determined, 
the photospheric radii at all other dates are known. To find the 
initial value, we go to the light-curve for the apparent magnitude 
and the temperature-curve for a temperature of corresponding date. 
At JD 2428800, m = +2.03 and T = 8400 K. The original appar- 
ent magnitude was +2.25, and an original absolute magnitude of 
—2.5 has been measured by Williams.? Consequently, the absolute 
magnitude at JD 2428800 was — 2.72. These quantities may be in- 
troduced into the well-known formula 
5900 


log R = — 0.20 M, — 0.02, (2) 
and the radius for this date found in terms of the radius of the sun 
as aunit. This value is 10.3 ©. Greaves’ states that the radius of 
a star determined from its color temperature and absolute magnitude 
may be in error by as much as a factor of 2 or 3 because the stellar 
radiation may differ decidedly from black-body distribution in the 


36 Ap. J., 88, 106, 1938. 
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inaccessible parts of the spectrum. However, this should not serious- 
ly affect the relative radii determined for y Cassiopeiae at different 
times. The constants in equation (2) are adjusted by observation on 
the sun, and several observers have shown that y Cassiopeiae radi- 
ated energy with black-body distribution in the spectral range of 


TABLE 1 
COMPUTED RADII AND TEMPERATURES 


JD Velocity 
km sec 
454 + 2™25 —2™50 10.3 © 11,050° K 
6500.. 454 2.57 10.3 11,350 
7000 42 +2.07 —2.68 11.0 11,250 
7100 a 405 +2.04 —2.71 11.6 10,900 
7200 297 +2.02 12.4 10,400 
7300 ed 304 +1.99 —2.76 12.8 10,250 
7.400 359 +1.97 13.0 10, 200 
7500 Her, 354 +1.95 —2.80 3.2 10,150 
7000 342 +1.92 —2.83 10,000 
7700 5 329 +1.90 —2.85 14.2 9g, 800 
7800 302 +1 .87 —2.88 9,300 
7QOo 330 +1.84 —2.91 10,150 
8000 308 +1.82 10,950 
8100 402 +1.79 —2.96 11.6 11,950 
8200 i 441 +1.77 —2.98 10.6 13,100 
8300 és 377 2.75 — 3.00 12.4 11,500 
8400 342 3.72 13.7 10,750 
302 +1.70 — 3.05 9,850 
273 +1.63 —3.12 9,35° 
8660 ae 260 +1.57 —3.18 18.0 9,250 
8700 . 260 +1.62 —3.13 18.0 9,100 
8775 270 +1.92 —2.83 3 8,500 
S800 278 +2.03 —2.72 16.8 8,400 
8860 311 +2.46 —2.20 15.0 8,000 
Sgoco 338 +2.33 —2.42 13.8 8,700 
goo. +2.30 —2.45 12.0 9,700 
Q100 433 +2.17 —2.58 10.8 10,950 
g200 477 + 2.25 —2.50 9.8 11,500 


4000 6300 A during all stages." °* From the above value of the 
radius the radii at all dates were computed and listed in Table rt. 
Then, with the computed radii and the absolute magnitudes derived 
from the observed light-curve, we may re-enter equation (2) and 
solve for the temperature at each date. 

The computed temperatures agree within the errors of measurement 
with the observed temperature-curve for the 1936-1938 cycle, although 
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there is a possible slight phase lag of computed curve over the meas- 
ured temperature-curve. If real, this means that the photosphere 
went through its changes ahead of the reversing layer. For 1938 the 
computed temperature is very slightly too low; and the original tem- 
perature, 1926, is found to be 11,050° K against a measured value of 
13,500 K.*° These discrepancies are in the right direction, however. 
If the reversing layers were larger in proportion to the photosphere 
in those times than in the period of great change, the difference is 
easily explained. It requires a relative decrease of only 12 per cent 
in the size of the photosphere to account for the 2500° K original dif- 
ference. The reversing layer would have to be largest relatively 
when the photosphere was smallest absolutely. 

A similar result could probably have been obtained on the as- 
sumption of an increase in general absorption in the atmosphere, but 
the fact that good agreement is found by the former method is con- 
sidered significant. 

For the numerical work the velocities used were those of //6, but 
with this limitation: because of the excess absorption strength on 
the side of the stronger emission component, it was felt that the 
velocities there might be affected. Therefore, the velocities were 
read from the smooth curve drawn through the points on the weak 
side only, thus shifting from side to side with the V, Rk curve. This 
precaution was hardly necessary, as a temperature-curve derived 
from the mean widths of the absorption velocity-curve was in excel- 
lent agreement with the adopted curve. 

Huffer’s photoelectric light-curve* and the observed and com- 
puted temperatures are shown in Figure 5. The observed tempera- 
tures are from three observatories and are reduced to the scale of 
Rk. C. Williams,*? of the University of Michigan. We have added 
0.15 and 0.12 to the Greenwich gradients’**? and to those of Cha- 
longe and Satir,?* respectively. Observations have been made by 
Vanderkerkhove’® *’ on the Greenwich scale but are not shown in 
the figure because of their relatively large scatter. The general agree- 
ment is good. The measured temperatures could have been reduced 
to the Greenwich scale equally well for such a computation, 

The temperatures, of course, can be computed accurately for only 


37 Pub. Obs. U. of Michigan, 7, 159, 19309. 


y CASSIOPEIAE 93 


those times when a light-curve is available. Thus, we are confined to 
dates after JD 2428400. Fortunately, the largest temperature varia- 
tions must have occurred after this date. Regardless of this, we may 
compute approximate temperatures for earlier dates. The apparent 
magnitude of y Cassiopeiae during the period 1931-1936 is not 
known; had it been of outstanding brilliance as in 1936, it would al- 
most certainly have been noticed. Likewise, it could not have been 
very subnormal, or this too would have been seen. The star was 
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Fic. 5. Photoelectric light-curve by Huffer (¢op) and theoretical temperature- 
curve (bottom). Published measured temperatures are shown. Crosses (X) are by 
Chalonge and Safir; open circles, by R. C. Williams; filled circles, by Greenwich ob- 
servers. 


often used as a standard of comparison for the variable a Cassiopeiae. 
Therefore, it has been assumed in the calculation that the brightness 
rose gradually, but continually, from 1931 to 1936, and a tempera- 
ture-curve was computed on that basis. This curve shows a second- 
ary dip in temperature at the time of the 1934 cycle, an effect which 
might have been suspected after the observed drop in the corre- 
sponding 1937 cycle but which certainly was not measured. If we 
assume that the apparent magnitude was normal in the interim, we 
also get a temperature dip, but deeper than on the above assump- 
tion. In any event, a mild drop in temperature for the weak 1934 


| 

\ 
| 


904 RALPH B. BALDWIN 


spectroscopic cycle is indicated, and a larger decrease for the radical 
1937 phase. 

A change in the original absolute magnitude changes the absolute 
values of the computed radii but not the relative values, and conse- 
quently the temperatures are independent of the absolute magni- 


tude. 
COMPUTED RADII OF THE EMITTING LAYERS 


Possibly the principle may be carried one step farther. On this 
rotational model the width of a spectral feature determines the radius 
or height at which it is formed, unless the lines are affected by radial 
motion. In a later section of this paper it will be shown that no sig- 
nificant pulsational velocities (greater than 60 km, sec) were present 
during a period of nearly two years. Consequently, we may deduce 
the general extent of the emitting layers by comparing their velocity 
widths with any reasonable assumed original velocity of rotation of 
the photosphere. This has been done for an original velocity of 300 
km_ sec, with the results listed in Table 2. If the original velocity 
were not 300 km/sec, the radii of the emitting layers will be in error 
by a constant proportionality factor. 

In Table 2 the radii for the different elements are derived from 
the measured edge displacements of each emission component.'? The 
hydrogen components are decidedly asymmetric, and the degree of 
asymmetry may depend on the depth in the atmosphere from which 
the Balmer emission can escape. If we can see down only to the 
photosphere in the Balmer wave lengths, the displacement of the 
red component will be measured too small in the case of an actual at- 
mospheric expansion. It may be necessary to invoke some such 
mechanism to explain the observed velocities and thus the emission 
asymmetry, and therefore the unequal extent of hydrogen emission 
on the approaching and receding sides of the star may be illusory. 

The hydrogen radii of Table 2 are shown in Figure 6, along with 
the V, R curve of emission intensity ratio for //6. It is found that the 
violet emission component is stronger than the red when the emit- 
ting layers are increasing in size, and that the red component is 
stronger during a decrease in the volume effective in emission. It is 
probably an apparent change of volume rather than a bodily shift 
of material, as the velocities are not observed to change sign. ‘The 
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atoms are probably continually streaming through the emitting lay- 
ers, and the apparent size of the envelope changes as the density and 


radiative flux vary. 
The only apparent exception to this rule comes about JD 2429100; 
but then the emission structure is complex, and neither the emission 


TABLE 2 
RADII OF THE EMITTING LAYERS 


HyDROGEN 
Rp Het Siu | Mgu | Niu 
R J 
97200... -| 12.4 1.93 18 21 ca pent 
79000...) 1.§R| 2.13R| 24 27 SG: SRO 
8000...| 12.7 1.7R| 1.68R| 22 27 
8100 11.6 1.8K) 1.72R)| 21 22 Tile’ 22 
8200. 10.6 2% 17 22 23 
12.4 1.60R| 1.18 27 19 23 18 22 220] 21 
8400. | 33 21 22 2 26 23 22 
8500 ES: 5 1.4 56 26 27 20 30 25 25 
8600...| 17.1 ‘7 2.70 52 37 34 28 Xe) 30 38 
8700.. 18.0 64 690 70 84 66 55 
8800...) 16.8 nad 88 110 58 84 107 52 81 
13.8 1.4R| 1.40R]| 70 48 53 51 62 50 53 
9006... .| (12.0 1.4R| 1.90R]} 32 30 26 29 30 41 40 
. 10.8 ? 24 26 23 2 22 20 20 
9200... | 9.8 1.4R| 2.57R| 27 30 22 27 22 23 2 


widths nor the V/R ratio estimates are reliable. There is an indica- 
tion of a subcycle in the proper direction in the absorption V/R 
curve, and the emission curve of V/R is usually in phase with that 
of the absorption. 

On McLaughlin’s hypothesis** we should find the opposite condi- 


38 Proc. Nat. Acad. Sci., 19, 44, 1933. 
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Fic. 6. Computed changes of photospheric radius and extent of emitting layers of 
hydrogen in the atmosphere (top). Hé emission V/R curve (bottom). Note that, when 
the violet emission component is stronger than the red, the hydrogen atmosphere is 
expanding, and vice versa. The hydrogen asymmetry may not be real. 
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Fic. 7.-Extents of emitting layers of five elements and variations in radius of the 


photosphere. Each of these layers appears to be symmetrical. 
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tions to obtain, the red component theoretically being the stronger 
one during an actual expansion of the atmosphere, while the violet 
is strong during a collapse. 

Figure 7 shows the changes in the radii of the emitting layers for 
the remaining elements of Table 2. The velocity-curves'’ of the edges 
of the emission components of these elements are nearly symmetrical 
so that mean radii only are given. 

If average values of the width, or height, of the emitting layers 
are taken for the period JD 2428300-2429200, the elements appear to 
be stratified as follows: Si highest, then H, with Niu, He 1, and 
Mg u grouped, and Fe 11 lowest. Similarly, the lines of the elements 
present in diffuse absorption suggest stratification. The S71 lines 
are narrowest, followed by He 1, Niu, and Mg u, which are grouped, 
with Ve 1 last. Hydrogen cannot be listed with the latter because of 
the slightly different manner in which the lines were measured. It 
should be borne in mind that a weak absorption line, or a line on an 
overexposed or underexposed plate, will be measured too narrow, 
while an emission line in the violet is apparently narrower than one 
in the blue because of slit effects. Attempts have been made to 
eliminate these effects, but both may still be operative, and thus the 
stratification may be somewhat other than that given. 


THE SHARP ABSORPTION LINES 


An outstanding and entirely unexpected set of spectral changes 
occurred between JD 2427880 and 2428500. These changes have 
been described by Lockyer,'® '° McLaughlin,"’ and in the second pa- 
per of this series.’ Briefly, the //e 1 lines, which are normally broad 
and dish-shaped, developed strong, sharp absorption cores. The 
first line to appear was \ 3889 (2°S—3'P), which arises from a meta- 
stable lower level. Soon other //e 1 lines appeared, including \ 5015 
(2'S—3'P), X 3965 (2'S—4'P), and Xd 3614 (2'S—5'P), which also 
possess a common metastable lower level. The last line was the only 
measurable feature in the region of the Balmer emission continuum. 
There were six other //e1 lines visible in more or less strength, as 
well as the Balmer lines of hydrogen. Several of the stronger lines 
of Fe 11 and \ 4067 of Niu were also present, although they were 
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extremely weak. The 2S level of hydrogen is metastable, as are the 
lower levels of the Fe 11 and Ni 1 lines. 


DATES OF APPEARANCE 

These sharp absorption lines appeared just after a red maximum 
of the emission V’/ Rk curve prior to the 1937 cycle. A search was 
then instituted on plates taken at the corresponding phase before the 
previous cycle. Suggestions of absorption cores at AX 3889, 3965, 
and 4026 of He 1, and at 176 were found. Thus, the absorption before 
the 1934 phase was weak, and this set of spectroscopic changes was 
not of great significance. The absorption lines suddenly became 
strong just before the 1937 cycle, and this was the most radical ever 
seen for any Be star. At present, in the fall of 1939, the emission 
VR curve is just past a red maximum, and the sharp absorption 
lines have returned in greater strength and numbers than ever be- 
fore. Many lines never until now seen in a “shell” spectrum are 
present. 

Two vears after the 1937 cycle came the beginning of a fourth 
cycle of spectroscopic changes. Thus far the shorter the period be- 
tween cycles the stronger have been the changes; the strength of the 
spectrum of the “‘shell” also seems to be an indication of the nature 
of the subsequent variations. Therefore, it seems likely that the 
1937 phase is about to be repeated, possibly in more cataclysmic 
form than before, i.e., in an increase in light and a corresponding de- 
crease in color temperature. 

THE VELOCITIES 

All of the “shell” lines which appeared before the 1937 cycle were 
extremely sharp, and even on microphotometer tracings showed 
practically no wings. Because of their signal sharpness, velocity 
measures were made on the two edges and centers of all lines. The 
results are given in Table 3. The helium and hydrogen velocities 
were similar but not identical. Irregular, though real, changes were 
found from plate to plate, and the He1 lines varied more than did 
the hydrogen lines. They are shown in Figure 8. Emission structures 
at the He 1 lines were present, but not measurable, on early plates of 
this period. Emission velocity measures for /7e 1 were made on some 
of the later plates, and for hydrogen they were made on all plates. 
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When the hydrogen absorption-line velocities are plotted over the 
emission structure, it is seen that they lie in the normal central dip 
of the emission line. However, when the helium absorption lines are 
plotted over the helium emission, it is seen that the absorption 


TABLE 3 


HYDROGEN AND HELIUM VELOCITIES* 


HyDROGEN HELIUM 

JD 

a EE | EM | VAE| CA | RAE| VAE] CA | RAE| VAE] CA | RAE 

(1) (2) (3) (4) (5) (6) (7) (8) (9) | (10) | (rt) | (12) 
807. ° + 2 — 64| + 
8182 ..| —22] —27| —35| 83] —4o] + 
8206 ..| —16} 66) —12] 
8256 —32| —35| —83} —42] + 2/— 86] — 2! 
8315 —32| —32| —32| - 3i— 75) — 
8342 —35| —40]| —97| —45| — —16 
8301 —27| —25|..... — 26 — 95} —16}| —21| +16 
83092 ..| —29]| —30] —7o| —20] +30/—103] —26)..... —33]..... 
$420 —~ 35) 32 — 28! —106} —67| —31)..... 7| 

| 


* All velocities are in kilometers per second. EE stands for mean emission edge velocity; EM, mean 
emission maxima velocity; VAE, violet absorption edge velocity; CA, central absorption velocity; RAE, 
red i edge velocity. 

Values in columns 7, 8, and 9 prior to JD 2428342 are mean velocities from all helium lines. After 
this date these columns give the velocities found for the lines \ 3889 and A 3965. The last three columns 
give the measured velocities for AX 3820, 4026, and 4472. 


gradually shifts until it partly overlies the violet emission component. 
Because there is no diminution of absorption intensity on the violet 
side, i.e., the emission is almost completely absorbed, we have addi- 
tional proof that the sharp lines of helium, and consequently of hy- 
drogen, originate in layers of atmosphere more remote from the 
photosphere than do the emission lines. The He 1 absorption veloc- 
ities show secondary changes super:mposed on a constant accelera- 
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tion by which the velocity changes from +10 to —6o0 km/sec in the 
six hundred days prior to JD 2428500. 

Peculiarly enough, the helium lines arising from metastable levels, 
d 3889 and d 3965, show this continuous acceleration throughout 
the entire period. The remaining //e 1 lines do not possess metastable 
lower levels, and they all show the acceleration until JD 2428342, 
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Fic. 8. .f: R lation between sharp absorption (dotted lines) and the two emission 
components of hydrogen. The two edges and center of the absorption line are given. 
Shaded regions show portions of emission which have been absorbed by the outer 
“shell” gases. B: Similar structure in lines of helium, showing an even greater over- 
lapping of absorption over emission. C: Mean velocities of the hydrogen emission 
edges. D: Mean velocities of the hydrogen emission maxima. /:: Mean velocities of 
the hydrogen sharp absorption lines. F: Mean velocities of the helium sharp absorp- 
tion lines. From JD 2427885 to 2428340 the velocities are the means of all lines. After 
the latter date the solid lines represent A 3889 and \ 3965, which arise from metastable 
lower levels. The dotted line is for the normal I!nes AX 3820, 4026, and 4472. 


when they break off sharply to lower velocities. The hydrogen lines 
also veer off, but not as far as the helium lines. This may possibly be 
a consequence of the partial metastability of the 2S level of hydrogen 
under the conditions of electron density to be found in the upper at- 
mosphere. 

An odd fact is noticed in this and other Be variables. The mean 


| 
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velocities of the emission edges or emission maxima practically du- 
plicate the velocities of the central dip in the emission. Here the hy- 
drogen emission edge velocities were almost identical with those 
measured for the sharp absorption, which arose in another layer, 
even to the secondary changes. This arises through a decided emis- 
sion asymmetry. The mean emission velocity curves are also shown 
in Figure 8. 
THE EQUIVALENT WIDTHS 

Through the kindness of Dr. J. F. Heard, of the David Dunlap 
Observatory, who made microphotometer tracings of plates of y 
Cassiopeiae taken during the second absorption-line stage, it has 
been possible to measure the equivalent widths of the sharp absorp- 
tion lines of helium and hydrogen over the entire sequence. 

The measured equivalent widths of these lines are listed in Ta- 
ble 4. The unit is one-thousandth of an angstrom. The column 
headed ‘3889 H7¢ marks the intensity of the blend of the helium 
and hydrogen lines. If, as seems probable, the hydrogen line is pro- 
duced lower in the atmosphere than \ 3889, the effective continuous 
spectrum reaching the helium atoms is reduced, and consequently 
the helium line is stronger than a simple subtraction of intensities of 
blend minus H¢ would indicate. The probable intensity of 1¢ may 
be interpolated from the intensities of Hy, H6, He, and Hn. A pre- 
liminary computation suggests that /7e 1 3889 is nearly of the same 
equivalent width as the measured width of the blend. Therefore, 
the intensity of the blend is a first approximation to the intensity of 
d 3889. The ratio of the strength of the blend to X 3965 is given in 
column 11. The measured mean ratio is 4.76. The theoretical ratio 
of strengths of AX 388g 3965 is 5. Consequently, the lines are pro- 
duced in a region so close to the photosphere that the dilution factor 
is not extremely small, and the difference in metastability of the 
2°S and 2'S levels is not yet apparent. With dilution factors such as 
are found in the Orion nebula the ratio of strengths may be very 
large. 

Unexpectedly, it was found that the X 3889-//¢ blend appeared to 
be weakest, relative to \ 3965, when the remaining hydrogen lines 


were strongest. 
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DISCUSSION 
Because the 2'S and 23S levels of helium are metastable and the 
2'P and 23P levels are not, it is possible to compute the dilution fac- 
tor by using relative intensities of absorption lines arising from these } 
levels in “shell” stars and in normal stars. The method has been | 
completely described by Struve and Wurm.*? 


TABLE 4 
EQUIVALENT WIDTHS OF THE SHARP ABSORPTION LINES IN 


THE SPECTRUM OF y CASSIOPEIAE f 
HELIUM HyDROGEN 
JD 
242+ 
88 388 
(1) 3820 3965| 4026] 4144] 4388] 4472] 5016] 5876 Hy | Hi | He | Hn 
(2) | (3) | (a) (5) | | (7) (8) | Gro) | | | (3) | (rg) | (5) 
8182.. 260) $85) 133) 154|....| 26 | 235]... ...| 4.40] 202) 166] 98} 149 
8206 625] 2591 200) ..... 3571... 2443] 473) 2751... 
8250.. 04} 923] 230] 252) 56 | 53 |] 3.86) 514) 322) 278 
. 125] 495} 80] 68]....] 30] 120)... + of 202) 177) 
So] 106)....]....] 3-64) 145] 92] 125) 55 
8461... 74 
|| 


In order to maintain homogeneity of the measures, a plate of y 
Pegasi, taken with the Michigan spectrograph, was run on the David | 
Dunlap Observatory microphotometer. The resultant equivalent 
widths were assumed to be those of a normal star and were used in | 
the computation. They are comparable to the values determined by l 
Williams.*° Rudnick" has also measured equivalent widths of lines | 
in the spectrum of y Pegasi; they show fair agreement with measures 

39 Ap. J., 88, 84, 1938. 

4° Ap. J., 83, 279, 1936; Mt. W. Contr., No. 540, 1936. 41 Ap. J., 83, 439, 1930. 
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Fic. 9. Equivalent widths of the sharp hydrogen and helium absorption lines 
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by Williams and the author for the weak lines but are systematically 
smaller for strong lines. Table 5 lists twelve lines whose equivalent 
widths have been measured. The unit of equivalent width is one- 
hundredth of an angstrom. 

From the dilution factor we may compute the radius of the helium 
shell relative to that of the photosphere, and if the radius of the 
photosphere is known, we shall have the absolute radius of the shell. 
This procedure is subject to rather large experimental error: never- 
theless, it should give an approximate value of the radius of the 


shell. 
TABLE 5 
EQUIVALENT WIDTHS OF STANDARD LINES IN THE 
SPECTRUM OF 7 PEGASI 


Line Baldwin | Williams} Rudnick Line Baldwin Williams] Rudnick 

8008.35 18 18 26 449 400 182 

40290. 143 98 1 ..| 409 510 

38 36 34 H&-3880 30% 529 

4388 98 80 Q2 


It should be recalled here that this process depends on the dilu- 
tion of black-body photospheric radiation at a particular color tem- 
perature. For the stars discussed by Struve and Wurm*? this method 
should give fair approximations to the real radii. In the atmosphere 
of y Cassiopeiae the case is different; there is detinite emission in 
helium lines whose lower levels have n = 2. Consequently, there 
will be strong emission at the transition 2'P—1'S. This line is emit- 
ted below the shell and thus will disturb both the effective color tem- 
perature of the radiation reaching the shell and the dilution factor. 
Therefore, the computed radius of the shell will be too small, and we 
shall determine only a lower limit of size. 

The radii of the photosphere are derived earlier in this paper; the 
values of the minimum relative radii of the shell, computed from the 
measured dilution factors, are given in Table 6, as are the minimum 
absolute radii. When the changes of radius of the shell are plotted 
on the same time scale as the mean intensity-curve of the sharp heli- 
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um absorption lines, they are seen to give nearly a reflected curve; 
i.c., the stronger the line, the closer to the photosphere it is produced. 

In making the computations the intensity ratios of AX 3965-4388, 
3905-4144, 3889-4026, and 3889-4472 were used. The assigned 
weights were 4, 2, 1, and 1, as \ 4144 was weaker than \ 4388, and 
since \ 3889 could not be observed independently. The value of the 
equivalent width used for \ 3889 was five times the observed width 
of X 3965; so we are really comparing \ 3965 with A 4026 and X 4472. 


TABLE 6 

CHANGES OF THE RADIUS OF THE SHELL 

JD Radius of ng Radius of 
242+ Photosphere Egy Shell 

Shell 

BOON 11.8 226 42 
. 11.6 38 
8128. Re: 11.4 3.0 34 
$152... sis 10.8 3.0 32 
8206 10.7 26 
82506. Pa 11.6 2.4 28 
$315. , 12.6 2.4 30 
8315. 12.6 2.9 37 
8342... 12.9 | 35 
8301. 13.6 30 
8392. 13.6 33 
14.3 3.1 44 


In each case the value was close to the measured equivalent width 
of the //¢-3889 blend. 

Earlier in the paper it was shown that changes in the effective 
volume of the atmosphere in which emission occurred could be found 
by comparing emission velocities at different times with an original 
photospheric velocity of rotation. Part of Figure 6 is repeated in 
Figure 10, with the minimum values of the radii of the shell super- 
imposed. It is seen that the minimum distance of the shell from the 
photosphere comes at the time when the emitting layer has reached 
minimum size. The shell is shown to be definitely outside, and its 
radius changes in a manner which is roughly proportional to the 
changes in size of the layers of emission. 
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The displacements of the He 1 lines show a continuous increase in 
negative, or expansional, velocity. In spite of this, the shell appar- 
ently approached the photosphere and then receded. Probably the 
shell is merely an optical effect produced by atoms streaming through 
a region where conditions are right for the production of sharp ab- 
sorption lines. The irregular nature of the intensity curves near JD 
2428340 is strong evidence against a discrete shell hypothesis. 

The shell layers are certainly 


+ 1000 7 and yet the » = 2 levels of he- 
L750 | lium are highly populated. This 
Boa i | means that radiation of wave 
length shorter than 584.4, 
1'S—2'P, is penetrating the lay- 
( (_) 4 ers where emission is produced. 
4 Thus, the Lyman continuous ab- 
oe \ | sorption is not complete. 
The sharp absorption lines 
L750 were present about a year anda 
is | half. In that time the maximum 
L 1250 4 velocity of expansion was ap- 


Fic. 10.—Computed radii of the helium about 30 km 
“shell,” compared with the radii of the During this period, then, it 
hydrogen emitting layers. The heliumlayer is clear that there were no ex- 
of emission is very similar to the hydrogen pansional components of motion 
layer, except that it appears to be sym- 
pa comparable to that of the ro- 

tation. The emission and ab- 
sorption V, R curves meanwhile changed sign. In both cases the red 
component was originally stronger, while the violet components 
gained as the red lost intensity, and were the more prominent fea- 
tures at the end of the shell stage. This seems to indicate that the ex- 
pansional component of the motion of the atmosphere is not the 
factor which controls the VR changes and that future models of the 
Be stars must be built on a modification of Struve’s rotational hy- 
pothesis.” 


42 1031: 
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THE V/R CHANGES OF THE EMISSION LINES 


B stars with emission lines are relatively rare; yet they actually 
total over four hundred established cases. In a large percentage of 
those which have been observed intensively the hydrogen emission 
lines are double and undergo cyclical V/R changes. However, in all 
the V/R variables only five stars are known to the author to show 
similar changes in the lines of any other element. They are 8 Mono- 
cerotis,*} where the variation is undoubted, Persei,**# Capricorni,*3 
the Be component of VV Cephei*® where the variation is suspected, 
and vy Cassiopeiae. In the first three cases it is only the strongest 
line or two of Fe 1 which shows V/R changes similar to those of 
hydrogen. For the companion of VV Cephei the H and K lines may 
vary, although the origin of these emission lines is not definitely 
established as in the atmosphere of the Be star. In the spectrum of 
y Cassiopeiae lines of at least ten different elements show the V/R 
transpositions. They are: H, Feu, Siu, Mg u, Cau, Niu, Cru, 
Nat, Het, a line tentatively assigned to Mg 1, and the unidentified 
line \ 6444. 

The VR curves for Hy and H6, three multiplets of Fe um, Sim, 
and Mg u, have been given for the period March 23, 1935—October 
24, 1938, in the second paper of this series." Estimates of the V/R 
ratios were made later for the K line of Ca 1, Ni 11 4067, and occa- 
sionally Cr 1 4558, when it could be distinguished from Fe 1 4556. 
Dr. Cherrington of the Perkins Observatory kindly sent several 
spectrograms of the visual region covering the time when the emis- 
sion lines were narrowing, and on these it was found that the D lines 
of Vat also showed the expected variations and that the range was 
larger than for any other element except hydrogen. The number of 
plates was not sufficient to give complete curves. Estimates were 
also made of \ 6318, which is possibly produced by Mg1, the D3 
line of He1, and \ 6444 of unknown origin. The latter line appears 
to be blended with another emission feature, as the V/R estimates 
were not dependable. 

Estimates of the ratios of intensities of the two emission com- 

43 McLaughlin, Pub. Obs. U. of Michigan, 4, 175, 1932. 

44 Dustheimer, Pub. Obs. U. of Michigan, 7, 171, 1930. 

4s Goedicke, Pub. Obs. U. of Michigan, 8, 1, 1930. 
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ponents have been made for all Balmer lines from H6 to HX. From 
these, two general conclusions may be reached. The lines of shorter 
wave length show a smaller range of V/R, and the changes at any 
line precede those of the line of next longer wave length by an inter- 
val which increases toward the long-wave end of the series. The 
phase differences for lines of shorter wave length than Hy are not 
noticeable. They are, roughly, 30, 60, and 80 days between He and 
H6, Hé and Hy, and Hy and H8, respectively, while Cherrington’s 
estimates* indicate a lag of the order of 300 days for Ha relative to 
Hp. 

In the rising part of the V/R curve for the 1934 and the 1937 
cycles, Hé preceded Hy; but in the descending part of the 1934 
curve, Hy preceded H6 and therefore passed through its changes 
more rapidly. An effect of the same nature is indicated for the de- 
scending branch of the latter cycle, but this cannot be proved, as 
Hy was not completely resolvable until near maximum width was 
reached. 

Numerous emission lines of Het appeared after the 1934 cycle 
and possessed both components. With the single exception of the 
D3 line in the visual spectrum, these components remained almost 
exactly equal in intensity. The D3 line showed changes similar in 
nature and direction to those of the other elements. 

The amplitudes of the VR curves for all elements were much less 
than those of the hydrogen lines. In spite of the relatively small vari- 
ations, real changes could be detected and estimated. These changes 
for all lines are shown graphically in Figure 11. For most elements 
the curves continue to rise until approximately the beginning of the 
single-line stage. The exact date of maximum is usually hard to de- 
termine; yet the Sz 11 lines and the lines of the b‘F —z'D° multiplet 
of Fe 11 reached maxima of V/ RK during the narrowing stage and then 
approached equality before the components blended. 


CORRELATIONS 
All the curves for the ionized metals are related to each other. 
Correlations may be made which suggest that the process producing 


46 Pub. A.A.S., 9, 2, 34, 1937- 
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the V/R changes is the same for all. The hydrogen, helium, and 
neutral metals give curves which are not as accordant. 

Only two elements are represented in the latter group—sodium 
and probably magnesium. Both the D lines of sodium are weak, and 
the blends caused by atmospheric absorption lines make estimates of 
VR difficult. 
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Fic. Emission V/R curves 


If we take the maximum values of the V/R ratios for the period 
JD 2428300 700 as determined from smooth curves drawn through 
the estimated points, it is to be expected that the accidental errors 
of estimation will largely cancel out and that an index of the true 
changes will be formed. This period of time corresponds to the nar- 
rowing phases when V was greater than R. The second period shown 
in Figure 11 is after the single-line stage, when R was greater than 
V. The latter values are not used in the discussion, as a second pair 
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of emission lines blended with the original components of several 
lines, perhaps with all, and seriously affected the accuracy of the 
estimates. The above-mentioned maximum V/R values are given 
in column 2 of Table 7; columns 3, 4, and 6 are self-explanatory; 
while column 5 lists the date when the two emission components be- 
came equal after R had been stronger than V’. 

It is seen that when the range of V/ R is plotted against the ioniza- 
tion potential there is a weak relationship. The elements which have 


TABLE 7 
DATA CONCERNING THE EMISSION LINES 


Maximum Ionization Date of 
Element V/R Potential Equality 
Range Volts Weight JD 242+ V/R 
JD 242+ 
Heit (D,) 24.47 4.000 8400 8625 
Mg il 1.55 14.97 24:32 8070? 8620 
S21. 1.43 16.27 28.06 8110 8475 
1.50 11.82 40.07 8150 8635 
Fe 11* t.3% 16.16 55.84 8175 8625 
Fe tt I.10 16.16 55.84 8325 8600 
16.16 55.84 8200 8435 
P25 18.4 58.60 8215 8625 
Mog 1? (6318) 1.20 7.61 24.32 8300? 8625 
*bsP—z1D>. bsF —z'Fe. 


the higher ionization potentials exhibit, in general, smaller VR 
ranges. On this chart the hydrogen values are not accordant (Fig. 
12, A). 

The atomic weights are also related to the range of V/ RX in that, 
excepting helium, the heavier elements are least variant. Hydrogen 
again shows too great a range to fit the curve, while the other resid- 
uals are within the errors of estimation (Fig. 12, B). 

An additional correlation may be noted. If the range of V, R is 
plotted against the Julian day when the emission components be- 
came equal after having had the red component stronger, there is a 
rclationship for all elements, again excepting hydrogen. The ele- 
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ments which showed the greatest range of V/R reached equality 
first. The dates for Mg 1 \ 4481 and Mg 1 6318 are extrapolations. 

The value for the line at \ 6318 falls beautifully on the last corre- 
lation curve; but if the other two are used to predict the atomic 
weight and ionization potential, we arrive at extremely improbable 
results. It seems that the first two curves are valid only for the ion- 
ized metals and that the placing of the values for the D lines is for- 
tuitous. If the relationships shown in the charts are due to physical 
causes and are not purely chance coincidences, they lend added sup- 


T T T T T 
30 45 60 /JD.2428000 ¢ 4 8400 
+18 = 


Fic. 12.-Correlations with emission V//R ratios. 4: V/R — Ionization potential. 
B: V/R — Atomic weight. C: V/R — Time when emission components became 
equal. Filled circles are for ionized metals; open circles are for neutral elements. 
Hydrogen values are indicated by arrows. The correlations appear valid for only the 
ionized metals. 


port to Cherrington’s identification’’ of \ 6318 as a line of a neutral 
metal. 
DISCUSSION 

In attempting to get a better picture of the structure of such a 
star as y Cassiopeiae from the above data, let us recall that Goe- 
dicke*> has shown for the Be companion to VV Cephei that the hy- 
drogen emission is stratitied. //8 is emitted with measurable inten- 
sity over a much larger volume than is Hy, and Hy over a larger 
volume than 1/6. As the V/R lags suggest a stratification effect in 
y Cassiopeiae, it is probably valid to assume that such is the case, 
and thus for hydrogen the range of V/R and the phases of the V/R 
curves are functions of the effective level of emission. 

The existence of the strong and broad //e 1 absorption lines in the 
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earlier stages suggests an original high excitation temperature, prob- 
ably higher than the observed original color temperature of 13,500° 
kK. As the color temperature dropped to 8000° K, the ionization 
temperature rose considerably. The recent presence and increasing 
strength of the helium emission lines indicate that the atoms were 
ionized and then underwent recombination. The ionization of the 
neutral helium atom requires nearly 25 volts, or a temperature of the 
order of 20,000° K. The increase in the emission strength"! of the 
hydrogen lines and continuum appears to be correlated with the 
change at the helium lines. Swings, Edlén, and Grandjean’ have 
identified four of the emission lines measured by the author during 
the single-line stage as due to Fe in. If their identifications are cor- 
rect, the iron atoms must first be doubly ionized and then excited to 
produce the observed emission. The weakness of the lines suggests 
that they are produced by excitation rather than by recombination. 
The presence of the Fe 11 lines probably is indicative of a high ion- 
ization temperature of the star during the period when the color tem- 
perature was very low. 

The production of emission lines of Fe 11 and //e 1 may be recon- 
ciled with the presence of the emission lines of neutral metals, such 
as Vatand Mg1, by the theory of Kosirev.** 

To say that the atmosphere of y Cassiopeiae appears to be strati- 
fied does not mean that the atoms themselves are confined to any 
given region. They may, and probably do, stream through these re- 
gions all together, as is shown by the similarities of their velocity- 
curves. In this view they would be stratified in the sense that they 
produced emission from definite layers in a manner similar to the 
Balmer lines. 

SUMMARY 

In conclusion it may be said that, although a complete physical 
model of y Cassiopeiae cannot now be constructed, the data present- 
ed point toward the fact that such a model must be based on rota- 
tional motions with pulsational velocities superimposed. The prin- 
ciple of conservation of angular momentum gives a satisfactory ex- 
planation of the correlations and dissimilarities between the observed 
light- and temperature-curves. 


47 Ap. J., 90, 378, 1939. 
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The size of the detectable atmosphere is greatly variable, and the 
emitting and absorbing layers are probably stratified. Certain ab- 
sorption lines characteristic of a ‘‘shell” spectrum are produced in 
regions outside the emitting layers. The conception of a discrete 
shell of hydrogen and helium must be discarded, as the atoms are 
observed to stream through a nearly stationary region in which the 
absorption lines are produced. 

The VR changes of the emission lines of ten elements are similar 
but show real differences which seem to be related to the atomic 
weights and the ionization potentials of the elements. Any future 
model must explain, for instance, why the curves for the lines of the 
three multiplets of /e 1 are not identical, and why elements whose 
velocity-curves are similar show VR curves going in opposite di- 
rections at the same time. A typical example of the latter point is 
found before the single-line stage, when the S7 1 lines were approach- 
ing equality after reaching a maximum of VR, while lines of Ca 1, 
Mg u, ete., were still showing an increasing ratio of intensities. 


I wish to express my appreciation to Dr. D. B. McLaughlin, of the 
University of Michigan Observatory, for lending me the long series of 
spectrograms of y Cassiopeiae and to Dr. Ernest Cherrington, Jr., 
of the Perkins Observatory, for lending me several spectrograms of 
the visual region. Dr. J. F. Heard, of the David Dunlap Observatory, 
rendered invaluable aid by making microphotometer tracings of the 
spectrograms taken during the “shell” stage. The page copy has 
been made possible by a subsidy from the Research Council of the 
Graduate School of Northwestern University. 
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THE KINETIC TEMPERATURE OF THE 
SOLAR REVERSING LAYER 


PEARL J. RUBENSTEIN 


ABSTRACT 


A redetermination of the curve of growth for the sun by means of least squares has 
resulted in a value of the temperature of the solar reversing layer for iron atoms 


T = 5400° + 1300° (m.e.) 
and a value of the damping constant 


1079 + X 107° . 

In a series of recent papers” 3 the theoretical expression for 
the curve of growth has been derived and compared with Allen’s* 
observations of the equivalent widths of lines in the solar spectrum. 
The present note concerns a further investigation of the same ma- 
terial and of the new tables by Allen,’ which extend the observations 
to 3924 A. A least-squares solution gave the theoretical curve which 
best fits the observational data. 

Theory shows that HW” can be given as a function of Y,, the opti- 
cal depth at the center of the line. For small values of X,, 


W 


log + log + log + log (1 Be: 4 


V2 
For intermediate values of X,4, 

2108 


W v 
log = log Xo log (1 + 


and for large X,, 


log = log — log 2 + } log + log + log , 


t Ap. J., 84, 462, 1936. 


2Ap. J., 84, 474, 1936. 4 Mem. Comm. Solar Obs., 1, No. 4, 1934. 
3 Ap. J., 87, 81, 1938. 5 Ibid., 2, No. 6, 1938. 
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KINETIC TEMPERATURE OF SOLAR LAYER II5 


where 
| 5.510 x 10 35° 


Since XV), the excitation potential, is essentially constant for lines 
of a given multiplet, all the lines have been plotted on a single dia- 
gram with an arbitrary zero point for log X,. Various multiplets are 
combined on the diagram by sliding the log XY, zero point until the 


TABLE 1 

log Xo log N O-C 
—0.92 —5.98 4 0.101 +o.10 
—5.82 14 . 242 .O4 
—0.53 —5.74 2 .121 — 
—0.32.... —5.5§2 17 .142 + 
—0.12.. —5.48 17 .148 — .14 
+0.68.. —5.02 30 -¥32 — .06 
+0.87.. —4.95 36 .149 — .06 
+1.08.. —4.84 25 .00 
+1.26.. —4.81 28 .O71 .02 
+1.44.. —4.75 21 + .o1 
+1 65. —4.71 .00 
+1.84.. —4.65 20 .078 + .02 
+2.05.. —4.58 20 .o81 + .03 
+2.43.. —4.46 25 .097 + .08 
+2.64.. —4.37 8 .136 + .07 
#3. 75: —4.47 4 O71 
+3.00. —4.18 .147 + .08 
+3.13. —4.20 2 .119 — 
| —4.03 ©.050 +0.02 


best fit with an assumed master-curve is obtained. The master- 
curve used in this study was that given by Menzel, Baker, and 
Goldberg.’ 

The resulting graph was divided into equal parts along the log X, 
axis, and all the points in each strip were averaged to give the nor- 
mal points for the least-squares solution. The first and second col- 
umns of Table 1 give the co-ordinates of the normal points; the 
third column, the number of observations: and the fourth, the mean 
error obtained from the combined dispersions in the two co-ordi- 
nates. The last column gives the difference between the observed 
and computed values of log (IW, \) for given log Y,, computed with 


116 PEARL J. RUBENSTEIN 


the least-squares values of the constants. Figure 1 shows the nor- 
mal points with their mean errors in the log (IV’/) co-ordinate and 
the final theoretical curve. 

A least-squares solution was carried out to determine the best 
values of the constants v/c and I'/y. An attempt to solve for both 
and separately resulted in values which differed insigniti- 
cantly from each other, in view of the large mean errors. In the 
final solution I'/y was therefore set equal to I)». The curve was 


WwW 
LOG 5 


——THEORETICAL CURVE 
@ NORMAL POINTS wl 
MEAN ERROR 


Qo +10 +20 +30 


Fic. 1 


divided into three parts, corresponding to the three theoretical ex- 
pressions; and points in the transition regions were omitted. 
Using the relation 


v = 1.289 X 10! Je j 


where yp is the atomic weight, one arrives at the following value of 
the temperature for iron atoms: 


= 5400° + 1300° (m.e.) . 


The solution also gives 


X 10°* + 0.40 X 10°° (m.e.). 


| 
i 
-50 
LOG Xx, 

| 
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From Figure 1 it is seen that the normal points fit very well on 
the theoretical curve determined. with two constants. The tempera- 
ture derived in the solution is not the same as that previously found, 
but the mean error is so large that one can draw no certain conclu- 
sions from the differences. One can merely state that the tempera- 
ture falls within the range of values determined by various other 
methods. The value of I'/y agrees well, within the limit of error, 
with that found by Menzel." 


I wish to express my appreciation to Dr. D. H. Menzel and Dr. 
T. E. Sterne for their kind assistance and many suggestions in con- 
nection with this paper. 


HARVARD COLLEGE OBSERVATORY 
February 12, 1940 
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RECENT EVIDENCE ON THE NUCLEAR REACTIONS 
IN THE CARBON CYCLE 
H. A. BETHE 


ABSTRACT 

New experimental results provide better estimates of the probabilities of nuclear 
reactions involved in the carbon cycle than those available previously. The lifetimes of 
C and N' are found to be almost equal, namely, 2.5 10° years and 4X 10° years, 
respectively. The abundances of C?? and V's should be similar in main-sequence stars. 
The combined abundance of C and V should be about 1 per cent by weight—in suf- 
ficient agreement with observational evidence. The ratio V'+:.V'5 in the sun should be 
2X10°:1, considerably greater than on earth (260: 1). 

In the last months several experimental investigations have been 
published which permit more accurate estimates of the probability 
of the nuclear reactions involved in the carbon cycle.' Curran and 
Strothers’ have observed the radiative capture of protons by NV" and 
have measured its cross-section. Burcham and Smith’ have found 
the reaction + H' = C” + and Holloway and Bethe? have 
measured its cross-section. The other reactions involved, viz., the 
radiative capture of protons by C’ and C%3, had already been 
measured** when the theory of the carbon cycle was put forward. 

The reaction 


NS + H' = C” + He! (1) 


has been found to be exceedingly probable. The measured cross- 
section is compatible with the assumption that all protons hitting 
the N'’ nucleus produce reaction (1). The effective width I used 
in “I’’ as a measure of the probability of nuclear reactions (cf. I, eqs. 
[4], [7], and [11] and Table V) comes out! to be 1.0X 10% e.v. (electron 
volts), while only 10° e.v. was assumed in I, Table V. The reaction 


™H. A. Bethe, Phys. Rev., 55, 434, 1939 (quoted as “I”’). 

2 Nature, 145, 224, 1940. 

3 Nature, 143, 795, 1930. 4 Phys. Rev., §7, 747, 1940. 

5 R. B. Roberts and N. P. Heydenburg, Phys. Rev., 53, 374, 1038. 

6 P.T. Dee, S. C. Curran, and V. Petrizilka, Nature, 141, 642, 1938. 
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is therefore one hundred times more probable than was assumed 
in I, which reduces the mean life of a N'> nucleus at the center of 
the sun to 20 years. At the same time it means that the reaction 
N® + H' = C’ + He*is presumably about 10° times more probable 
than the radiative capture of the proton, NV’ + H' = O' + ¥, be- 
cause the “width,” I’, is of the order of 10 e.v. for most radiative 
capture reactions. Of a million C’* nuclei undergoing the transfor- 
mations of the carbon cycle, on the average only one will be trans- 
formed into O", while all the others will return to C’?. In I, the ratio 
was assumed to be only 104, which was already quite sufficient to in- 
sure the conservation of practically all the carbon and nitrogen in the 
star. 
The measurement of the reaction 


is even more important. Curran and Strothers? find that about 
1.510" nuclei O'' are formed in a thick Va,N target by protons 
of 960 kv. (kilovolts) initial energy. The excitation curve shows indi- 
cations of resonances, so that conclusions from the high-energy 
(500 goo kv.) cross-section on the probability of thermonuclear re- 
actions (proton energy, a few kilovolts) are not very certain. From 
the curve given in footnote 2, the cross-section comes out to 
be 1.8X10°3° cm? between 890 and 940 kv. proton energy, and 
0.33 107-3" cm? between 550 and 650 kv. From either measurement 
we obtain I’ =60 e.v., using equation (11) of I. This is a very high 
probability for a radiative capture reaction which is equaled only by 
the reaction:’ Li? + H' = Be’ + y. If we assume the same I for 
low proton energy (which may not be justified because of the reso- 
nances), the probability of the thermonuclear reaction (2) will be 
twelve times greater than was assumed in I, and therefore the mean 
life of a N'4 nucleus at the center of the sun will be reduced to 
4& 10° years (the previous value was 5X 10’). 

Including the carbon isotopes for which we take over the values 
obtained in I, we have for the lifetimes of the nuclei involved in the 


7W. A. Fowler, E. R. Gaerttner, and C. C. Lauritsen, Phys. Rev., 53, 628, 1938. 
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carbon cycle at the center of the sun (assumed density, p=8o; 
hydrogen content, 35 per cent; temperature, 2X 10’ degrees): 


2.5 X 10° years | 

(3) 

4 X 10° | 

20 | 


All these values may easily be wrong by about a factor of 3 either 
way. The main conclusions to be drawn from the figures are: 

1. The lifetimes of C’? and NV" are now almost equal, whereas 
previously the life of V'* was assumed to be twenty times longer than 
that of C’. This means that C’ and N™ should be about equally 
abundant in main-sequence stars, which seems to be in agreement 
with the astrophysical evidence, in contrast to the previous assump- 
tion (V'+ twenty times as abundant as C”). 

2. The total time required for a complete carbon cycle is now 
only 6X10° years, as against 5X10’ years previously assumed. 
Therefore, in order to obtain the same energy evolution as before, 
we can now assume the combined abundance of carbon and nitrogen 
to be only one-tenth of the value assumed in I. This means a com- 
bined abundance of about 1 per cent by weight, or about 0.5 per 
cent for each of the two elements, C and V. This again seems to be 
in sufficient agreement with astrophysical data, whereas the abun- 
dance of 10 per cent for nitrogen, assumed in I, appeared to be much 
too high. 

3. The ratio of the abundances of the two nitrogen isotopes in the 
sun should be V'4: NV = 2 X 10°:1. This is very much greater than 
the abundance ratio on earth,* viz., 260:1. It is, of course, not vet 
quite certain whether our values for the lifetimes of N‘4 and N* in 
the sun are correct, but the value for V'’ seems rather well estab- 
lished from the nuclear observations, while that for V‘4 seems well 
substantiated by the astrophysical evidence (points 1 and 2 above). 
If we assume the lifetimes to be correct, there is certainly a very 
great discrepancy between the abundance ratios V‘':V" at the 


8 A. L. Vaughan, J. H. Williams, and J. T. Tate, Phys. Rev., 46, 327, 1934. 
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center of the sun and on the earth. This makes it likely that the 
earth was separated from the sun when the equilibrium between the 
elements in the carbon cycle was not yet established and the sun was 
still at a lower temperature. 

Similar discrepancies have previously been found? in the abundance 
ratio of the hydrogen isotopes and in the existence of lithium on the 
earth. The ratio H':H? is about 5000:1 on the earth, while it must 
be at least 1o'*:1 at the center of the sun. Li is found in some 
quantity on the earth and must be exceedingly rare at the center of 
the sun of the same order as deuterium. In all cases mentioned 
(NS, H?, and Li*?), nuclei which should be very rare in the sun, 
because they undergo rapid nuclear reactions, are found in some 
abundance on the earth. The simplest explanations would be that 
cither (a) the earth was separated from the sun when the sun was 
much cooler than now and the respective nuclear reactions had not 
yet destroyed the H?, Li, and N's nuclei originally present or (0) 
the earth originated from some region of the sun near the surface 
which was cool enough to retain approximately its original chemical 
composition. Assumption (>) requires that there be very little mix- 
ing of the various parts of the sun, so that concentration gradients 
can persist. Even though this may be true ordinarily, it seems not 
very likely at the time of a catastrophe such as the formation of the 
planetary system. Assumption (a) may seem more attractive; more- 
over, an early separation of the earth from the sun is also indicated 
by the approximate equality of the age of the earth (from radio- 
active measurements) and of the universe (from the red shift). More- 
over, it seems that such a freak event as the formation of the plane- 
tary system could have occurred more easily at a stage when the 
density of stars in the universe was larger and the probability of 
collisions correspondingly greater than at the present time. 

DEPARTMENT OF PHysIcs 
CORNELL UNIVERSITY 
April 1940 


9. Teller and G. Gamow, private communication. 
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SPECTROPHOTOMETRIC OBSERVATIONS OF THE 
MINIMUM OF ¢ AURIGAE 


C. A. WIRTANEN 


ABSTRACT 

The variation in brightness of ¢ Aurigae during the recent eclipse has been observed 
photographically at eight different wave lengths, using the 1o-inch Cooke camera 
equipped with a coarse grating and a low-dispersion prism. Plates were measured with 
the Schilt microphotometer. Seven early-type stars in the same field, ranging from 3.3 
to 7.3 mag., served to establish the zero point and the scale at each wave length. 

The light-variation of ¢ Aurigae was observed spectrophotometri- 
cally by crossing the low-dispersion prism on the 1o-inch Cooke 
triplet with a coarse objective grating. The aperture was reduced to 
25 inches to improve the color-curve of the objective. With this 
aperture, measurable spectra of seventh-magnitude stars can be ob- 
tained with a 1-hour exposure. The spectra were widened by means 
of a fine screw which moves the plateholder carriage in right ascen- 
sion, the carriage being shifted by 0.01 mm once every 6 minutes 
during the exposure. On the original plates the width of the spectra 
is 0.12 mm, while the interval from 6 to He is 2.7 mm. A part of 
the field is reproduced in Plate III, in which ¢ Aurigae is shown in 
three phases with one of the comparison stars. Eastman Super 
Panchro-Press plates were used; they were developed 15 minutes in 
Rodinal 1:32 at 68° F., with potassium bromide added to reduce 
background fog. 

The plates were measured with the Schilt microphotometer, using 
a 1-mm diaphragm, which produces an illuminated area 0.2 mm in 
diameter on the plate. Two regions in good focus were selected for 
measurement, the yellow part of the spectrum being measured at 
two points on either side of \ 5500, and the violet part at six points 
from Hy to Hn. Settings on these eight points were made with an 
accuracy of about 0.03 mm. Figure 1 shows the general appearance 
of the “tracings,” the upper curve being that of a Bog star. 

The magnitudes of ¢ Aurigae were read from characteristic curves 
determined by the central images and the grating images of the 
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PLATE III 


a 
Cc 


¢ AURIGAE AND THE REFERENCE STAR HD 32188 AT THE RIGHT 


(a) November 13, 1939; maximum. (/) January 28, 1940; partial 
eclipse. (c) January 6, 1940; total eclipse. 
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seven comparison stars. Three plates in different regions of the sky 
yielded the following values for the grating constants: 


First Order | Second Order 


1.90+0.04 3.04 40.05 


For the two points near \ 5500, Harvard visual magnitudes reduced 
to the I.A.U. photovisual system' were adopted for the seven com- 
parison stars. Photographic magnitudes were obtained by applying 


20 


30 


40 | ZETA AURIGAE MAX | | 
| ” MIN | 


Fic. 1. ~Tracings of ¢ Aurigae at maximum and minimum referred to the same zero 
point. In the topmost tracing of HD 31069, Bo, the zero point is abritrary. The vertical 
lines indicate the points measured. 

I 


King’s color indices.’ These were regarded as suitable for the points 
in the violet, since the average spectrum of the comparison stars is 
close to Ao and since there is practically no progression of spectral 
type with magnitude. The comparison stars are listed in Table 1. 
Eight reduction curves, one for each wave length, were used for each 
plate. In every case of normal exposure the magnitudes of ¢ Aurigae 
in the yellow were determined from its grating images; magnitudes 
in the blue and violet were determined from its central image. 

The magnitudes, as derived from each plate, are given in Table 2. 

* Harvard Obs. Mimeograms (3d ser.), No. 1. 

2 Russell, Dugan, and Stewart, Astronomy, p. 734, 1927. 
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TABLE 1 


Star Sp. mMpv Mpg Star Sp. Mpv Mpg 
42080: B3 | 3.30 | 3.06 || HD 32206... AS 
33462... B8 | 6.94 | 6.87 
31069... Bg | 6.02 | 5.98 31866. . | 
AO: ||) to 
TABLE 2 


1940 Jan. 


Probable error of one entry 


5.16. 


Wave Lenctus 
5660 5470 4260 4180 4040 3940 3910 3870 
3.61 | 3.82 | 5.20 | 5.34 | 5.40 | 5.64 | 5.72 | 5.65 
3-05 | 3-77 | 5-23 | 5-37 | 5-52 | 5-090 | 5-75 | 5-00 
3.56 | 3.87 | 5.26 | 5.52 | 5.46 | 5.63 | 5.74 | 5.77 
3-73 | 3-81 | 5-33 | 5-49 | 5.52 | 5.78 | 5.82 | 6.09 
3-65 | 3.81 | 5.23 | 5.36} 5.53 | 5.72 | 5-79 | 5.78 
3-59 | 3-79 | 5.34 | 5-42 | 5.58 | 5.74 | 5.71 | 5.86 
3-54 | 3-79 | 5-27 | 5-47 | 5-40 | 5.74 | 5.82] 5.71 
3-71 | 3-93 | 5.40 | 5-52 | 5.06 | 6.02 | 6.23 | 6.29 
3.78 | 3.95 | 6.10 | 6.35 | 6.52 | 7.53 | 7.50 | 7.35 
3257 | A503] 5.621 6.13 | 6.16 
3-77 | 3-83 | 5.37 | 5.58 | 5.67 | 6.15 | 6.33 | 6.05 
0.04 | 0.05 | 0.04 | 0.04 | 0.05 | 0.07 | 0.08 | 0.09 


* Twelve-minute exposure. 


t Plates fogged, poor seeing. 


|| 
— 
G.M.T. 
1939 Dec. 12.19.... 
| 
17.09 
18.09.. | 
18.96. | 
19.00. 
Ig.10 | | 
EO. | | 
19.35 | | 
| 
} 
21.04. | 
21 . | 
27: 
26.13T 
26.10T. } 
29.05... 
27:10... 
28.04.. 
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The bottom line gives the average internal probable error at each 
point. This was derived from the interagreement of the plates ob- 
tained during one night, excepting the two plates taken on January 
26. On these two plates the fogging due to the full moon and the 
poor definition due to bad seeing combined to make the settings with 
the microphotometer uncertain. A graphical representation of the 
observations is given in Figure 2. It will be seen that the points in 
the violet do not present a steady progression with wave length. 


DECI7Z 18 19 20 JAN6 27 28 
Oo 
x 
05 —  \ 
\ \ 
\ + + 
¢ 25660 \ 
O 5470 \ 
Q 4260 \ & 
20— + 4 — 20 
3940 a 
& 39/0 
3870 


Fic. 2. Graphical representation of Table 2. Mean values of each night are sub- 
tracted from the means of December 12 and 15. 


This is due to the fact that the points at \ 3940 and at \ 3910 are 
within the heavy wings of the H and K lines in ¢ Aurigae. The meas- 
ures are in good agreement with those of other investigators at previ- 
ous eclipses, as may be seen by a comparison with the summary pub- 


lished by Hall.s 


I wish to acknowledge a number of suggestions made by Dr. 
EK. T. R. Williams, and particularly to thank Professor Vyssotsky 
for his continued interest after having suggested the problem. 

LEANDER MCCORMICK OBSERVATORY 


UNIVERSITY OF VIRGINIA 


3 Ap. J., 90, 317, 1930. 


NOTES 


NOTES ON PROPER-MOTION STARS* 


This note is a continuation of two earlier ones.’ Table 1 consists 
of two parts, the first of which contains 38 subdwarts found recently 
by the writer. This list brings the total number of known subdwarfs 
to 140. The second part of Table 1 gives eight of the ten stars listed 
by Luyten’ as having a low luminosity. The remaining two stars 
could not be identified from the rough co-ordinates given by Luyten. 
It is seen that the white dwarf is confirmed, being of the 40 Eridani 
type; but some of the other stars are in no way exceptional, and only 
two (or possibly three) of them appear to be subdwarfs. 

Six binaries were discovered. They are Wolf 531 (AB=12.5 Mo; 
15", Am = 1.2); Wolf 676 (13.0 M3+, 3”, A = 23); Wolf 1495 
(11.3 M2, 3”, A = 3); Ross 534 (A = 11.5 K6, B = 14.3 M3, 7”); 
Ross 486 (10’"); and Wolf 1561 (8”). The last two pairs occur in 
Table 2. 

The white dwarf Wolf 672 noted earlier’ was found to be of the 
40 Eridani type; and its parallax, from the companion, 0704. 

Over 1500 proper motion stars have now been observed by the 
writer for spectral class, including most of the stars with large trigo- 
nometric parallax accessible from this latitude. Using these stars for 
purposes of calibration, spectral parallaxes were derived for all stars 
observed; stars with spectral parallax not less than o’09 for which 
no trigonometric parallax is as yet available are given in Table 2. It 
is estimated that the probable error of the spectral parallaxes is 
about 20 per cent, apart from possible errors in the apparent magni- 
tudes; the latter are mostly estimates by the writer. 

The planet Pluto was observed with the same spectrograph, on 
panchromatic film. No absorptions such as shown by the outer plan- 
ets could be seen. 

* Contributions from the McDonald Observatory, University of Texas, No. 23. 

Ap. J., 89, 548, 1939; 91, 269 and 366, 1940. 

2 Harv. Ann. Cards, Nos. 521 and 525, 1940. 
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TABLE 1 

Star a 6 Myis Sp m (Sp.) 
Karpov Anon.. .... 0:04.1 + 8:29 12.9 1710 K4: 0.007 
Ross 570..... E36 — 7 30 10.7 0.34 A6 002 
—11 18 12.0 0.36 Go: 003 
Ross 797... . 5 49.4 —14 24 0.53 Fo 005 
Ross 612... 6 56.3 + 6 34 12.0 ©.50 G6 005 
Hertzsp. 11. 8 37.7 +36 37 11.5 0.43 Go 004 
Furuhj. 24. 8 58.9 +39 13 11.8 0.47 F7 003 
Wolf 334... 9 51.9 +33 06 12,2 ©.49 G2 004 
Ross 892... . IO 22.3 + 1 55 0.38 F8 004 
Ross 106... . IO 44.2 +56 59 12.6 0.51 Kr 005 
Furuhj. 31. Ir 05.8 +45 58 E225 0.74 Kr 005 
+36° 2165 07.3 +36 17 9.9 0.52 F4 005 
+14°2481... I2 07.0 +13 50 10.2 0.40 F6 005 
Wolf 420*..... I2 27.4 +12 31 14.3 0.37 K3: 003 
Wolf 1447... I2 29.4 +15 50 $827 0.35 A6 ool 
Wolf 449.... 1250.7 +13 06 Lies 0.38 F6 003 
Ross 979. . 13 05.7 — 8 20 13.6 0.48 Ko: 003 
L 1338-3 .. 13 27.9 +24 44 12.6 0.96 K3 007 
Wolf 508 . 13 43.8 +16 00 0.39 Ko 004 
Ross 1025... S204 +25 58 0.42 G5 004 
Ross 1027.... 13 56.4 +48 09 12.5 1.45 K4 009 
Ross 993... +41 59 0.68 K3 006 
Ross 1038... sal E5-6%-60 + 6 25 12.0 0.96 G2: .004 
Ross 510. ES +17 04 0.97 .009 
Ross 805. . A 2606 —II 02 ©. 46 F8 .004 
Wolf 624. 15 56.4 + 5 41 12.8 0.51 F8 002 
Ross 527... .:|. £6.24.6 +48 23 13.0 G8 003 
Furuhj. 45. | +44 55 0.75 G5 006 
+17°3154 +17 20 8.8 2.2 G1 O15 
Wolf 754... + 1 06 13.4 0.47 K3 005 
—8°4501. — 8 44 10.8 0.46 F5 004 
Ross 135. —24 37 11.8 0.50 G1 004 
Wolf 816. . ..| 18 05.0 +20 o1 12.5 0.53 Kr 005 
Wolf 834. . +20 29 43,2 1.00 K5 00g 
Wolf 1465. 18 30.4 — 8 20 14.1 1.28 K4: 004 
Ross 711. . 18 31.4 +28 38 II.4 0.36 A8& 002 
Wolf 851. . +0 52 12.3 1.90 K4 
Ross: 767. +34 36 F4 002 
Kpv 1126-68 cal + 7 59 15.0 0.94 
Wolf 536. . 14 33.1 02 12.8 0.56 K3 006 
—15°1776 oy 7 16.6 —15 10 9.7 0.54 G5 O13 
12 +16 17 1.3 1.37 K4 006 
Ross 484..... 10.8 0.65 K4 o18 
Kpv 1126-25 + 9 05 11.6 0.50 Ko 007 
L 1064-75... 18 11.3 + 1 29 12.5 0.73 Ms 09 

| 


* Subdwarf if motion is real; stellar, not nebulous. 


TABLE 2 


Star 
a 6 Myis Sp x (Sp.) 
(= \). +60: 50 10.8 M3 10 
Wolf.a7(=B)._. 48 | 9.5 | 0.74 | 
Wolf 66. . 13.7 0.82 | M7 14 
Ross 10....... 12.0 0.46 | M5 
Ross 555. +56 45 12.2 | Ms 
Ross 15. . +58 17 0.79 | M4 
Ross 550 0.30 | M5 13 
Wolf 1322... 2 aha | | 30:7 | 
Ross 504. . 6 49 o 8 | Ms og 
Wolf 1530 pee 39 39 13.3 | 0.53 | M6+ 14 
ark 4 40.7 + 6 19 11.8 0.22 Ms 13 
52 -< . 
Ross 42 153 03 O20 1.96 | Mi 09 
Wolf 237... + 9 45 11.4 0.30 | M6 
Ross 79 +82 06 1.30 M2 10 
Ross 417... 4 21 0.93 My’ 18 
— 44°3045. 0.67 | Ms+ 10 
+ 36°1638 (=A) 43 12.0 1.12 | Ms 13 
Ross 882.... o.g1 | 10 
L. 829-19 AB.. 8 +348 0.64 | M6 22 
Ross 686 (=A). 8 6 0.62 Ms 
Ross 687 (=B) 8 6 5 1 0.47 Ms 
| 9 +35 5! 13.0 0.33 | Ms 60 
R | I2.2 0.08 17 
OSS 104 | 10 54.8 Aas 7 
SS | I 22 
Ross 689... | 12 00.6 | 3 11.0 1.40 Ms5 23 
0. 12 18.4 +64 35 0.73 M4 
oss 095. 312 19.6 —17 28 
Wolf 414. | 12 23.9 3 | 2.53 M4+ 13 
Wolf 30 IO 23 0.50 M4 12 
Ross 845. . | 8 I 5° 14 3 0.51 M5 Os 
Ross 848 IT 33 13-5 0.79 Mo oy 
| 14 14.6 | = O06 Ms 
Ross | 11.4 M4 10 
1449.5 | +23 58 o7 
Ross 802...... 5 74 M5 13 
Ross 868 (=A) | 16 M5+ og 
Ross 867 (=B) 3 0.47 Ms 
. 4 | 42 3 5 fo) 7 I 
Furuhj. 48. | 17 16.7 +41 50 M 
L 1064-75 | 18 11.3 ) 
Wolf 108 “9 2.5 0.73 Ms on 
2040.5 | +5457 | 15.3 | 1.87 | M7+ 
OSS 205. 57-0 +16 07 
+0°4810 | 2 ©. 30 M3 10 
Wolf 1561B | 13-5 0.08 Mo og 
5 ‘| 223201 9 18 0.68 M7 
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